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Guide to the reader 



When writing this text, I had a two-fold aim and audience in mind which is roughly re- 
flected in the division of the text into a main part and the appendices. 

First of all, I wanted to give a gentle introduction to (elements of) quantum probability 
theory using quantization of the damped harmonic oscillator as leitmotiv. This is the focus 
in the main text. As audience I had (graduate) students in physics in mind, people like 
myself when I started the research of this thesis, although also others with similar level or 
interest are welcome to have a look. 

My aim was to explain some of the physical intuition and motivation behind the, sometimes 
overwhelming, math machinery of quantum probability theory in order to convey some 
of its intrinsic appeal. To this end, we start in the main text from a piece of quantum 
mechanics that is treated in almost all basic textbooks: the quantization of the (undamped) 
harmonic oscillator from the Heisenberg and Schrodinger point of view. We show how 
both treatments are special instances of a quantum probabilistic quantization procedure: the 
second quantization functor. As a corrollary, quantum probability theory is a generalization 
of the quantum mechanics of Heisenberg, Dirac, Schrodinger etc. as we know it. We then 
proceed to build up the quantum probability machinery needed to quantize the damped 
harmonic oscillator. I have tried to stress the ideas and intuition in the main text and defer 
the details of the proofs and reasoning to the appendices, but I have not been entirely 
consistent in that. 

Secondly, this thesis also reflects the research I have done. Actually, parts of this research 
were used and extended in lectures that my thesis supervisor, dr Hans Maassen, gave in a 
summer school in Grenoble in 1998. In particular, Proposition 4.1 and Proposition 6.3 of 
his lecture notes [11] correspond to my theorem 2.5 and lemma 1.30. The details of their 
proofs are given in appendix B and A. 18, respectively. 

It goes without saying that my own writing has benefited greatly from the theoretical expose 
that dr Maassen provided in his lecture notes [11]. 

Steven Teerenstra 
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Chapter 1 

Quantum Mechanical Harmonic 
Oscillator 

1.1 Dirac's Quantization procedure 

Following Dirac's quantization rules [5], we translate the (1 -dimensional) classical har- 
monic oscillator into a quantum mechanical one: 

Observables The classical conjugated (generalized) coordinates for position and momen- 
tum, q and p, are replaced by self-adjoint operators Q and P (acting on some Hilbert 
space !tf). Any observable of the classical harmonic oscillator a = f(q,p) is re- 
placed by its quantum counter part A = f(Q,P) (which is hermitized if necessary). 
In particular, the Hamiltonian of the quantum mechanical oscillator is 

rj _ , mXS?Q 2 
n ~ 2m 2 ' 

Correspondence principle The relation between two quantum observables A = f(Q,P) 
and B = g(Q,P) of the oscillator system is expressed by their commutator [A,B] : = 
AB — BA. If a = f(q,p) and b = g(q,p) are the corresponding classical observables 
having Poisson bracket 

{a,b} = ^f(q,p)l p g(q,p) - l q g(q,p)l p f(q,p) 

then the correspondence principle holds: 

[A,B] = ih{a,b}t. (1.1) 

As a special case, Heisenberg's commutation relation emerges 

[Q,P]=ihl. (1.2) 

Time evolution If A = a(Q,P) is an observable that does not depend explicitly on time 
(i.e. does not depend on time other than via the time evolution of Q en P), then the 
classical description of the time evolution, ^ = {a, h} where h is the classical Hamil- 
tonian, is translated via the correspondence principle into the Heisenberg equation of 
motion 

having A (f) = e >tH/H Ae -itH/h as f orma i so i ut i n. 

In particular, ^ = P/m and ^ = -Q. 

States The states of the system are described by (normalized) vectors |oc) in the Hilbert 
space M . 
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Measurements Any single measurement of an oscillator observable A yields a value a 
in the spectrum of A (as self-adjoint operator) i.e. there is a normalized eigenstate 
\a) such that A\a) =a\a). Suppose for simplicity that the state of the system is 
completely determined by A, then these eigenstates of A form a complete set in H 
i.e. every state \Q can be decomposed as a "superposition" of eigenstates of A: 

|q I La \ a ) ( a \Q if the spectrum of A is discrete, ^ ^ 

1 Jda \a) (a\Q if the spectrum of A is continuous. 

Quantum theory states that when the system is in state |Q, the probability of mea- 
suring the value a for A, i.e. exactly a in the discrete case or within the resolution A 
of the measuring apparatus in the continuous case, is 

| (a\Q | 2 if the spectrum of A is discrete, 

■I""— A/2 I (^10 \ 2 daK A| (a\Q \ 2 if the spectrum of A is continuous. 

Directly after having measured the observable A, the system is in the eigenstate \a) 
belonging to the value a (collapse of the wave function). The probability of measuring 
another value a 1 directly after, is then zero, since a (a\a') = {Aa\a') = (a|A T ) a' = 
(a\Aa') = a' (a\a') which implies (as either a ^ or a' ^ 0) that (a\a') = 0. 

Expected value On average, measurements of A on identical systems, all in state |Q, will 
yield an average value of the observable A, which is called the expected value of A 
and is denoted by (A)^. 

Quantum theory states that (A)^ = (C|A|Q / (QQ, which can be justified from the 
probabilistic interpretation ( 1 .5) above. For instance, if A has a complete set of eigen- 
vectors and L, is normalized, then 

(QA\Q = La.a(Qa) ( fl |A|5) (51Q = E«£ <C|fl>flS fliS <5|Q = S>| (Qa) \ 2 

in the discrete case, while 

(QA\Q = Jda dd (Qa) (a\A\S) (5|Q = fda dd (Qa) a8(a -a)(a\Q 
= Ja\(Qa)\ 2 da 

in the continuous case. In both cases, we see that the possible outcomes of measure- 
ments a are weighted by their corresponding probabilities in state 

Maximal system of commuting observables: quantum numbers For the one dimensional 
harmonic oscillator, every eigenvalue q of the position operator Q is non-degenerate 
i.e. there is essentially one eigenstate \q) such that Q\q) = q \q): all vectors V|/ £ 9{ 
that satisfy Q\\t = q\\t are scalar multiples of \q). As a convenient consequence, there 
is an observable A such that measuring the observable means that we know that the 
system can only be in one state, viz. the (one and only) eigenstate \a) that corresponds 
to the measured eigenvalue a. In general, it is rare for a quantum system to have one 
observable that completely determines the system in this sense. More often than not 
there are more essentially different (i.e. linearly independent) eigenvectors for a given 
eigenvalue a i.e. the eigenspace of a consisting of all vectors in 9l on which A acts 
as multiplication by a is not one-dimensional. For example, the x-position of a two- 
dimensional oscillator does not uniquely determine its state. In such cases, one looks 
for observables {A,},- that commute pairwise: [A,',A ; ] = 0. For such observables, an 
eigenvector of one A; is at the same time an eigenvector of anotherA ; - [1, p. 207], [17, 
p. 30], and therefore a simultaneous eigenvector. In particular, such observables can 
be measured simultaneously. The simultaneous eigenvectors can be labeled by their 
eigenvalues a; ,A2, . . . with respect to the commuting set of observables A\ ,A2, . . . i.e. 
A,- |ai,fl2,fl3, . . .) — a, |ai,fl2,fl3, . . .). A set of commuting observables {A,},- is called 
maximal if there is no observable B such that [B,A,] = for all / (and, of course, B 
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should not be a function of the {A/},). If we have a maximal system of commuting 
observables {A,}, then the joint eigenvalues (a,),, enable a unique labeling of the 
joint eigenstates |oi,02, ■■■)■ These labeling eigenvalues 01,02, • • • are referred to as 
quantum numbers. In that case, measuring all the {A,},- simultaneously determines 
the state of the system completely, namely (01,02, . . .). E.g. in case of the two dimen- 
sional harmonic oscillator {Q x ,Q y }, {P x ,Py}, {Qx,Py}, {Qy,Px} all form a maximal 
system of commuting observables. 

1.2 Representations of the Canonical Commutation Rela- 
tions and the problem of uniqueness 

The above procedure does not provide a clue as to which Hilbert space one could choose 
and how the self-adjoint observables Q, P, and H of the oscillator system would look like. 
Starting from a set of generalized coordinates of the system, some general clues can be 
given, but as those clues have a high "who was first: the chicken or the egg" content, first 
quantization remains all in all a mystery. 

First quantization In general, first quantization of a system follows the following steps. 

1 . Take a set X representing the degrees of freedom of the system i.e. if the system has n 
degrees of freedom and generalized coordinates q\,...,q n with o, <G X t , then for X we 
can take X = n ( X ( = { (q\ , . . . , q n ) : o,- e X t for all i}. 

Examples: 

- if we are considering a free particle with 3 degrees of freedom and take for (01,02,03) 
the x, y and z position of the particle (or if we take for (#1 , #2 , #3) the x, y and z com- 
position of the momentum): X = R 3 . 

- if we are considering a spin ^-particle (so 1 degree of freedom q e { — 5, j}): X = 

1 2'2J- 

2. Form the Hilbert space 9i ~L 2 (X) consisting of the complex square integrable func- 
tions on X. 

3. The quantization of the system is given by: 

- states of the system at time t are (wave) functions V(/(X,f) e L 2 (X) with 

J x \\\i(x,t)\ 2 dx=l, 

- observables of the system are self-adjoint operators on L 2 (X). In particular, the 
observable g, attached to measurement of the i th generalized coordinate is diagonal: 
Qi is the multiplication by o,- in L 2 (X). 

Using the above hints, we now recall two familiar representations of the quantum mechan- 
ical harmonic oscillator: the Schrodinger and Heisenberg representation. 

1.2.1 Schrodinger representation 

As generalized coordinate, we take the position q of the oscillator. Then X = M, M — 
L 2 (M), and Q = M x i.e. the multiplication operator by x in L 2 (X). As Q is to be a self- 
adjoint operator, we must specify a domain and we could, for example, take a "maximal" 
domain: 

(Q\\i)(x) =x\\i(x) for ye Dq := e L 2 (R) : f R \xty(x)\ 2 dx < <*>} . 

For the conjugated coordinate P, one could take for example the derivation operator with a 
factor hj i thrown in (no excusing explanation offered here) and again a "maximal" domain: 
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WW = fsVW = -ifiVW for ¥ e ©p := { V e L 2 (R) : / K | ¥ '(x)| 2 < 
If one feels up to it, then one can check that, with the above definitions, Q and P are indeed 
self-adjoint and that the commutation relation [Q,P] = iht holds on the intersection of their 
domains. 

The Hamiltonian H=£-+ boils down to 



n 2m 7? 



+ 



2m 
t(a 2 x 2 



At this moment, we just want to note that there is a high degree of arbitrariness in the choice 
of Q, P, H and their domains. 

Also note that the symmetry between P and Q in the Hamiltonian H = |^ + mm 2 e is 
maximal when (units are chosen such that) m = CO = 1. Then H = ^ + 



1.2.2 Heisenberg representation 

As generalized coordinate we now take the energy e of the oscillator. Then X C [0, °°) and 
H is a multiplication operator on L 2 (X), but matters turn out to be not as simple as putting 
X = [0,°o) wAH=M x . 

To deduce Heisenberg's representation of the harmonic oscillator, most textbooks on quan- 
tum mechanics (e.g. [17, p.90-92]) follow the factorization method due to Dirac who rewrote 
the Hamiltonian as 



with 



A:= 




m(i>J 



and A 




(1.6) 



By formal calculation, the anti self-adjoint operators A^ and A are then shown to lower and 
raise the eigenvalue of H respectively, and are therefore called annihilation and creation 
operators. In particular, the set of eigenvalues of H (i.e. the quantum energy levels) turns 
out to be discrete, consisting of E„ = fti£>(n + ^),bgN:= {0,1,...} and the action of A^ 
and A on the corresponding eigenstates \n), n € N, is given by 

A f |n) = V^TT|«+1) (n = 0,l,...), 

A\n) = y/n\n-l) (n = l,2,...), (1.7) 
A|0)=0, 

i.e. the action of A on the state (of an existing oscillator system) with no energy quanta 
is to destroy the system (indicated by 0). Although Dirac's algebraic method does not 
use an explicit representation it is easy to conceive one by taking X = N, !H = £ 2 (N) = 
{(Co,Ci,...):ir=olC«| 2 <-}and 

|0> = (1,0,0,0,...) 

|1) = (0,1,0,0,...) 

|2) = (0,0,1,0,...) 



Using that Q = ^h/lm&iA* +A), P = v /fi/2mco(A t — A), H = h(£>(A f A + f 1) and the 
known action of A, (1.7) one obtains 



Q = 



2m(i3 
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1.2.3 Uniqueness of representations 

Having now (at least) two representations of the quantum mechanical harmonic oscillator, 
the question arises whether the are "the same" i.e. whether the "physical content" of both 
representations is the same. By the "same", we mean of course that the observable physical 
quantities calculated in either representation are equal. In quantum theory, the physical 
observable quantities of a system are the expected values (A) = (£|A|Q of the system's 
observables A in a given state of the system It can be shown then [1, p. 210-212] that 
representations have the same physical content if they are unitarily equivalent i.e. if states 
£ and observables A in the one representation are related to states and observables A' in 
the other by one and the same unitary map U: 

A i-> A' = UAU~ l , C^C' = ^C- 

Although perhaps not evident at this point, the Schrodinger and Heisenberg representation 
of the harmonic oscillator are indeed unitarily equivalent. We will show this later, after 
having developed theory which makes both representations simple examples of general 
"tools" to make representations (that are automatically equivalent). 

For the moment, however, we confine ourselves to noting that it has been rather coincidental 
that we have found two pairs (Q,P) that are unitarily equivalent indeed. To illustrate this 
statement, we will provide a reasonable definition of canonical pairs and show that such 
canonical pairs are not necessarily unitarily equivalent. 



A reasonable but failing definition of canonical pairs 

Let us start with the the following reasonable definition. 
Definition 1.1 (Heisenberg canonical pair) 

A pair (Q,P) of self-adjoint operators on a Hilbert space H is called a (Heisenberg) canon- 
ical pair if there exists a dense domain DcWso that 

(i) Q:D^D,P:D^D (i.e. D is stable), 

(ii) QPiSf - PQiSf = /(j) for all (j) e D, 

(Hi) Q and P are essentially 1 self-adjoint on D. 

Two canonical pairs (Q,P) on H and (Q' ,P') on tt' are called unitarily equivalent if there 
exists a unitary map U : 9l — » H' such that Q' = UQU^ 1 andP' = UPU^ 1 . (In particular, 
this imposes domain requirements: U~ lr Dp — "Dpi and U~ 1 Dq — D@). 

'i.e. although Q and P are not necessarily self-adjoint on the domain D, they do have unique self-adjoint 
extensions. 
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The following example, due to Nelson [15, p. 275], [20, Remark 3, p. 90 ; Problem 4, 
p. 94], shows that the above definition does not uniquely determine one pair of self-adjoint 
operators. 

Example 1.2 (Nelson) 

Let M be the Riemann surface of sjz and let M = L 2 (M) equipped with the Lebesgue 
measure (locally). Let 

13 ^ 13.13 

Q :=x,P :=-j-, and Q' :=x + -— , P' := -— , 
i ox i ay i Ox 

be defined on the domain D that consists of all infinitely differentiable functions with com- 
pact domain not containing the origin. 

Then (Q,P) and (Q 1 ,P') are canonical pairs that are not unitarily equivalent. In fact, the 
unitary groups generated by Q' and P' respectively do not satisfy the Weyl relations 1 .9 
further on, while the unitary groups generated by Q and P do. 



The Weyl- Von Neumann formulation of the canonical commutation relations 

Apparently, the canonical commutation relation as formulated in (1 .2), p. 9, does not deter- 
mine the position and momentum observable of a quantum system uniquely. To circumvent 
the domain complications associated with the self-adjoint character of Q and P, Weyl's idea 
was to formulate the canonical commutation relation (1.2) of Q and P in terms of bounded 
operators (which are everywhere defined on For this he used the one-to-one corre- 
spondence between self-adjoint operators and one-parameter, strongly continuous unitary 
groups. 

Theorem 1.3 (Stone's theorem) 

Let 9i be a Hilbert space. If A is a self-adjoint operator on 9i , then U t := e" A is a unitary 
operator, U t+S = U t U s for all s,(Gl (i.e. (U t ) t is a one-parameter group), and (j) e ^ and 
t — > fo implies U t § — > U tQ § (i.e. 1 1— > U t is strong continuous). 

Conversely, ifU t (t £ R) is a strongly continuous 1-parameter unitary group on H , then 
there is a self-adjoint operator A on H so that 

U t =e itA . 

See [15, p. 265-267] for a proof. 

Inspired by this correspondence, Weyl proposed to formulate the canonical commutation 
relation (1.2) of Q and P in terms of the associated unitary groups. Letting 

T, := e- i,p , and S s := e isQ , 
and formally applying the Campbell-Baker-Hausdorff formula (A.6): 

e\p(A + Bj) =exp(A)exp(B)e-2M] if [A, B] is a scalar, (1.8) 
one formally derives the Weyl relation: 

T t S s = e- ihst S s T t (s,t€R). (1.9) 

Taking this one step further, Von Neumann combined the two unitary groups into a two- 
parameter family of Weyl operators: 

W{t,s):=e ! r st T t S s = e- ! T st S s T t (j,/St), (1.10) 
that satisfy the Weyl relation as expressed by 

W(t,s)W(v,u) =e-^ {m - sv) W(s + u,t +v) (s,t,u,v G M). (1.11) 
Note that by the Baker-Campbell-Hausdorff formula 
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e ix(A+B) = e ixA e ixB £ - \ [iA.iB] ^ (A 6)) 

we have 

W(t,s)=e isQ - itP . (1.12) 

To compact notation and prepare for generalization further on, we reformulate the Weyl 
operators by writing z = t + is £ C 

W(z) = W(t + is) = W{t,s) = e isQ ~ i,p = e ^m-^(z)P = e ^' TtSs = e -%*S s T t . (1.13) 
The Weyl relation then takes the form 

W(z)W(w) =e-^ lm(Iw) W(z + w) {z = t + is,w = v + iueC). (1.14) 
Definition 1.4 (Weyl canonical pair) 

When we refer to a family of unitary operators W(z), z <E C, as Weyl operators, then we 
implicitly assume that the family is strongly continuous and that the Weyl relation (1.14) 
holds. 

The pair of generators Q and P of the strongly continuous unitary groups s i— > W(is) and 
1 1— > W(t) respectively, are called a Weyl canonical pair. 

In contrast to the Heisenberg form of the commutation relation, the Weyl form of the com- 
mutation relation does uniquely determine the two unitary one-parameter groups T t , f el, 
and S s , seR, or, equivalently, the two-parameter family W(z), z <E C ([15, p. 279], [11, p. 
32-35]). 

Furthermore, a Weyl canonical pair satisfies the conditions of (1.1) [14], [15, p. 275]. In 
conclusion, Weyl canonical pairs satisfy Heisenberg's form of the canonical commutation 
relation, but are in addition unique (up to unitary equivalence). For the cases we consider 
we will infer this uniqueness property from a general purpose tool that we will need and 
develop later. 

For the time being, let us note the following: 

- In the Schrodinger representation of the harmonic oscillator (jK = L 2 (M)), we have that 
T t := e~ ,,p is a shift to the right (over an amount fit), S s :=e ls ® is the multiplication by e' sx 
and the Weyl operators are W(t + is) := e ls Q- ,,p — e'~r st T t S s . The commutation relation 
between Q and P is expressed by the Weyl relation W(u)W(z) — e~'~r lm ( uz )W (u + z) for 
the Weyl operators. 

- Similarly, in the Heisenberg representation, a family of Weyl operators z i— » W(z) exists 
such that Q is the matrix of the generator of s W(is) in the energy basis and such that 
P the matrix of the generator of 1 1— ► W(t). 

- The operators Q' and P' in Nelson's example above do not satisfy the Weyl relation. 

As a consequence of the general theory we develop later, the Schrodinger and Heisenberg 
representation of the harmonic oscillator are unitarily equivalent, and the "representation" 
derived from Nelson's example is unitarily equivalent to neither of them. 

1.3 Coherent states 

The non-vanishing commutator of Q and P prohibits to measure both simultaneously with 
arbitrary precision. Striving to maximize precision in simultaneous measurement of (Q, P), 
one arrives at coherent states, which first appeared in Schrodinger's treatment of the har- 
monic oscillator [18]. Coherent states turn out to be closely related to states with minimal 
and equal uncertainty in Q and P. As we will see, coherent states furnish an context in 
which it is easy to establish the equivalence of the Schrodinger and Heisenberg represen- 
tation of the harmonic oscillator. Apart from that, coherent states form the link to classical 
states and, more importantly in our context, a natural entry point to second quantization 
which is to be discussed in the next section. In fact, we will be build all our theory on 
coherent states. 
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The uncertainty principle for Q and P 

Let us first recall the proof of the uncertainty relation AQ -AP>j. 

For an observable A of a system in state the uncertainty AA of A in state |\|/) is defined 
via 

(AA) 2 = ((A - (A) ) 2 ) = (\|/| (A- ( ¥ |A| ¥ ) ) 2 | ¥ ), 

i.e. the average deviation 2 of A from its expected value (A) = (\|/|A|\|/) in state |\|/). 

If A, B are two observables (self-adjoint operators), then their uncertainties AA, AB are 

related to their commutator. Indeed, as the quadratic polynomial 

f(a) := a 2 (A 2 ) +a(i[A,B]) + (B 2 ) = ((aA + ;B) ¥ |(aA + /B) ¥ ) = ||aA + ;fi) ¥ || 2 (a e C), 
is non-negative, its discriminant ((/[A,B])) 2 — 4(A 2 ) (B 2 ) must be non-positive: 

(A 2 )( J B 2 )>I|([A, J B])| 2 . 
Replacing A and B with (A — (A) ) and (B — (B) ) respectively, and observing that 

[A - (A) , B - (B)} = [A,B] , (as (A) and (B) are scalars), 
we find the uncertainty product relation: 

(AA) 2 (AB) 2 > i|([A,B])| 2 or equivalently, AA AB > i|([A,B])|. (1.15) 

Note from the above argument that minimum product uncertainty arises, when the discrim- 
inant of / is equal to zero, i.e. when / has one zero. In that case there is an a G C such that 
(aA + iB) = 0, where A = A - (A) and B = B - (B) . 

We now apply the above knowledge to the pair (Q,P). First, the familiar product uncer- 
tainty relation AQ ■ AP > -| emerges. Secondly, states of minimal (product) uncertainty in 
Q and P arise when (aQ + iP) |\|/) =0. 

Observe that if | ¥ ) is a state with (g)y = t and = s, then we can translate it via 
T_ t/h = e itp l h and S_ s/h = e~ isQ l h to obtain a state | ¥o ) with (g) ¥o = (P) Vo = 0. We then 
have (see A. 8) that \\\t) is a solution of 

(ag + iP) \Mf) = (a(Q -t) + i(P - s)) | ¥ ) = (*) 
if and only if |\|/rj) is a solution of 

(oce + iP)|Vo> =0. (**) 
Therefore, we can find the solutions of (*) by translating solutions of (**) in position and 
momentum. 

Let us find the minimal uncertainty states with (Q) = (P) = in the position representation 
(so P = P = ti/ijj and Q = Q is diagonal with eigen states Q\q) = q\q), q € R). The 
minimum product uncertainty condition then reads: 

(aQ + iP))y (q) = cu 1 yo(q) + fi^r 1 (q) = ( ? el), (1.16) 

aq 

where ¥ o(?) = (?| ¥ o) is the position-representation of state | ¥ o)- The solution of this 
differential equation is exp(— jj-q 2 ) which is normalizable if Re (a) > 0. Including an 
appropriate normalization constant (which we take to be real-valued), the solution reads 

VoteH^W-^ 2 ) (?GR,Re(a)>0). 

Then 

2 more precisely, it is the root-mean square deviation 
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(aq) 2 = (e 2 ) = ( V ie 2 | ¥ ) = /\ 2 l¥o(?)i 2 ^ = vW/? 2ex p(-^ 2 )^ 

(42) T 

2Re(a) ' 

and since \|/ has the minimum product uncertainty {AQ) 2 {AP) 2 = j,we have 
(AP) 2 = hRe 2 ( - a '> 

Therefore, the uncertainty in Q and P is equal, whenever Re (a) = 1 . The wave function of 
a state with (Q) = (P) = and AQ = AP equal to their (simultaneous) minimum \Jh/2 is 
therefore 

Vo(«) = ^exp(-^)expH^) (,6R). 

for some a = 1 + ;Im (a) € C. Observe that the choice of a does not change the probability 
density |\|/o(g)| 2 as it only adds a phase. The usual convention is to take a real i.e. a = 1 
and thus 

„2 



Yo(9) = -^=exp(-£-) («6R). (1.17) 



1 ^ 

If we take units such that mco = 1, we see from (1.16) and (1.6) that |\|/o) is a solution of 
A |\|/o) = i.e. |\|/o) = |0), the ground state of the oscillator which has a Gaussian form in 
position representation. 

We now derive the form of minimal equal uncertainty states with non-zero (Q) = t and 
(P) = s. As mentioned before, these are obtained by translating the position (via T t / h = 
e -aP/hy an( j me momen tum (via S s / h = e ls ®l h ) of the ground state |0) and they are precisely 
the solutions of 

([Q-t) + i{P-s))\y)=0, (1.18) 
which, after taking units such that mco = 1 in (1.6), boils down to 

, , , Z . , , , . \ . f Re (z) is the expected position, and 

A = —== \\\f) where z = (f + is) i.e. )(. F F ' (1.19) 

\/2h \ Irn(z) is the expected momentum. 

In position representation, equation (1.18) takes the form 

+ = zy(q) (q e R), 

which can be directly integrated to (including a normalization constant) 

Xf{q)=Ne- ilt & with |AT| = ^-^ ( 9 g R) 

or can be obtained by shifting the position and momentum of the ground state to t and s 
respectively. 

To motivate the choice of normalization constant N, we first observe that, with Vfo(<?) as 
in (1.17), the order in which the translations in position (e' s ® = S s ) and momentum space 
(e" p = T t ) are applied matters: 

(i./.or^t?)^^^-^, while 

(r, /ft os l/ftVo j(9) = 3==-^-^. 



In order to retain as much symmetry between Q and P as possible, the usual convention is 
to "translate in position and momentum space simultaneously" via our old friend 

e {isQ-i,P)lh = w{{t + is)/h) ^ e^T t/h S slh = e-n*S, /h T t/h . 
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Convention 1.5 

The minimal (equal) uncertainty state with position t and momentum s is defined as 

MinUnc(z) = e^ Q - itP ^ h \\t = W((t + is)/H)y . (1.20) 

In particular, in position space representation it takes the following form: 

(g|MinUnc(z))= e te / 2 * * L^ '^ (1.21) 

which is indeed the "middle way" between first applying T t j h and then 5 S /^ or vice versa: 
MinUnc(z) = e- ist / 2h S s/h T t/h \]f = e +ist / 2H T t/h S s/h \]f . 

Equation (1.19) tells us that, if in units are chosen such that mco = 1, MinUnc(z) is an 
eigenvector of A for eigenvalue z/ V2h. Things simplify in terms of the eigenstates of A. 

Convention 1.6 

Coherent states are defined as the eigenstates of the annihilation operator A, and labeled as 

ACoh{z)=zCoh{z). 

If (units are chosen such that) mco = 1, therefore in particular for the standard Hamiltonian 
H = (Q 2 +P 2 )/2, Coh(z) is a minimum uncertainty state of the harmonic oscillator, more 
precisely Coh(z) is characterized by 

i. expected position (Coh(z)\Q\Coh(z)) = V2ftRe(z), 

ii. expected momentum (Coh(z)\P\Coh(z)) = \/2Mm(z), 

Hi. equal uncertainty in Q and P such that the product uncertainty is minimal: AQ = AP = 



(1.22) 



since MinUnc(z) = Coh(z/v2fi.). 

Coherent states in Schrodinger representation 

In position representation, we have from (1.21) 

(<7|Coh(z)> = - j= e -(W5M 

„-{{\z\ 2 -z 2 ) ^_ , 

Also, from (1.19) and (1.20) 

Coh(z) = MinUnc(\/2fc) = W(V2hz/h) |0) = 

= W(-£=)\0). (L23) 

WW 

Furthermore, a calculation shows (see A. 9) 

(Coh(w)|Coh(z)) = Jdq(Coh(w)\q) (?|Coh(z)) = e'^ e'^ e Wz . (1.24) 
In particular, the ground state expectation of a Weyl operator is 

/o|W(-^==)|o\ = (Coh(0)|Coh(z)) =e-^\ (1.25) 



Since 

(9|Coh(fa)) = £*I^etol>/W ( q g R) 

we see, by the uniqueness of the Fourier transform on L 2 (R), that the set {Coh(w) :sgK} 
and thus the set of coherent vectors is total in L 2 (K). 

From a mathematical point of view, the coherent states Coh(z) are normalized versions of 
so-called exponential states: 



1 .3. Coherent states 
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e(z):=eil z l Coh(z), 
which read in the Schrodinger representation (use equation (1.22)) 



(9|Coh(z)) = ^e-(«/^ft" z ) 2 (q g R). 

Mathematically, exponential states are easier to use in calculations, which for example is 
expressed in their inner product: 

(e(w)|e(z)) 

(use (1.24)). The latter property and their completeness in fact determines exponential 
vectors up to a unitary transformation as we will see later (see section 1.14). 

Coherent states in the Heisenberg representation 

In the Heisenberg representation, the eigenvalue equation A |z) = z |z) can be solved by 
expanding \z) in the energy basis, 

|z>=I" =0 oc„|n) 
and rewriting the eigenvalue equation as 

l^ =Q za„\n) = z \z) =A\z) =Lk =1 a k Vk\k-l) 
so that a recursion relation for its coefficients emerges 

which has 

oc„ = ^a 

as solution. The coefficient ao can be determined from normalization: 
l = £|a n | 2 =£^|a | 2 = X|a | 2 , 

n=0 n=0 

so that the coherent states read 

| Z ) = e -w 2 /2 y J? \ n) <L 7) e -N 2 /2 gZ At |Q) _ (1 26) 

«=o v«! 

The coherent states in Heisenberg representation are total. Indeed, the linear span of { \z) : 
z € C} contains the vectors jc(z) := £n=o |n) which are total as follows from the lemma 
below. 

Lemma 1.7 (Exponential vector in Heisenberg representation) 

Define forz <E C, the exponential vector 



rc(z):=ir=0 7^l">- 

Then 

i. (%(z)\%(w)) =e lw (z,weC), 



u. 



^ r n(0] (=o = V^|n>. 



Idt" 
(see (A.12.1). 

Again the coherent vectors can be generated from |0) by the Weyl operators if units are 
chosen such that mco = 1. To see this, assume that mco = 1 and observe from (1.6) that 

isQ-itP = is(^/Kfl{A+A^ -it(^/h/2{A-Ai)/i} = (zA f -zA) , 

where z = t + is as usual. By the Baker-Campbell-Hausdorff formula, 

W( z ) = e isQ - itP = eV^i^-^) = e Vm^\-y/h/2u e -^ly/mzA\-y^iA\ 
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hence W(z/y/K/2) = e^V^e"^ 2 . Since e~ 1A |0) = (1 + I~=i =jjp) |0) = |0> +0, we 
see that 

W{z/^hj2) |0) = e ^ t e -^ e -jN 2 10) = e -iWV f |0) 



The distribution of the number of quanta A' = A ' A in a coherent state is given by the prob- 
abilities p n of measuring n quanta 

Pn = | («|Coh(z)> | 2 = | £„,«„, («|n') | 2 = |a„| 2 = J^V^I", 

which describes a Poisson distribution with parameter |z 2 |. In particular, the expected 
number of quanta and the variance thereof are \z 2 \ in a coherent state |z). 



1.3.1 Passage from physical to quantum probability 

In the formulas that appeared on stage to now, Planck's constant h figured as h/2, h, or 2H, 
and all of these are multiples of h/2. Therefore, formulas simplify if we set h/2 = 1 i.e. 
h = 2, and the latter is a convention often adopted in quantum probability theory. 
In particular, we then have the following reductions. 

1. The Weyl relation (1.14) reads: 

W(z)W(w) =e- ilm ^W(z + w) (z,weC). (1.27) 

2. The coherent vectors (1.24) in the Schrodinger representation (and as we will see 
later: also for cyclic vacuum representation of the harmonic oscillator) are obtained 
by applying Weyl operators to the ground state: 

Coh(z)=W(z)|0) (zeC). (1.28) 

3. The expectation of the Weyl operator in the ground state (1.25) in Schrodinger's rep- 
resentation is 

(0\W(z)\0) =e-2M 2 . (1.29) 

Again, we will see later that this holds for so-called cyclic vacuum representations of the 
harmonic oscillator. 



1.3.2 Different kind of representations 

Summarizing we have seen three kinds of representations for the quantum mechanical os- 
cillator: 

1 . The Schrodinger representation diagonalizes the position operator and may therefore 
be called the configuration space representation. This representation emphasizes the 
wave like character of the quantum mechanical oscillator, which is why it is also 
called a "wave representation". 

2. The Heisenberg (Fock) representation diagonalizes the energy operator. The corpus- 
cular nature of the quantum mechanical oscillator is evident is this representation, 
which is the reason why it is also called the particle representation. 



1.4. An apology for functorial quantization 



21 



3. The Bargmann-Segal representation is less well-known than the former two. This 
representation "diagonalizes" the creation operators: when we identify the phase 
space with the complex plane via q + ip <-> z, then states can be realized as analytic 
functions on the complex plane and creation operator becomes the multiplication 
with z. One could call this "complex wave" representation [3] [2, p. 23] the phase 
space representation. We will not work explicitly with this representation, but this 
representation is always lurking in the background, because quantization of fields is 
most elegant in this representation. The latter will become clearer when we look at 
the second quantization functor and its interpretation, see chapter 2. 

1.4 An apology for functorial quantization 

Although Dirac correspondence principle is a familiar and respected procedure to quantize 
the harmonic oscillator, some of its inherent difficulties could already be tasted. First of 
all, there is no direct recipe how to translate classical observables into self-adjoint opera- 
tors. Although physical reasoning and intuition often can provide formal expressions for 
the operators, one still faces the (sometimes formidable) problems of choosing a suitable 
Hilbert space and finding a domain on which the operator can be proven to be self-adjoint. 
In addition, sometimes more than one self-adjoint operator is a reasonable candidate for 
the observable and one has to choose the "right one" [15, p. 303]. A second and for our 
context worse problem is that the Dirac procedure breaks down for damped (so dissipative) 
oscillators, because it requires a Hamiltonian to generate the quantum time evolution. In 
this sense, the Dirac correspondence principle is only applicable to conservative (so at least 
energy conserving) classical systems. 

This section introduces another way of quantizing the harmonic oscillator to deal with the 
above shortcomings. This approach can be traced back to work of Simon [19, Chapter 1, 
paragraph 4]. The basic idea is to quantize the time evolution in the classical phase space 
via a functor (the second quantization functor F). When applied to the harmonic oscillator 
evolution, it reproduces well-known results, which gives us faith that the procedure makes 
sense in other situations as well. 

Besides mathematical elegance, the motivation for functorial quantization is its ability to 
deal with dissipative systems in a non ad-hoc manner. 

We start with a section to develop the bits and pieces that get the functorial machine going. 

1.5 A mathematical toolkit for Hilbert spaces and unitary 
maps 

Hilbert spaces and unitary maps are corner stones in quantum theory. The former, because 
the (pure) states of a system can be described by vectors in a suitable Hilbert space and the 
latter, because physically equivalent representations are linked by unitary maps between the 
representing Hilbert spaces. Moreover, the observables of a quantum system are generators 
of unitary one-parameter groups living on the representing Hilbert space. 
This section provides general tools to generate Hilbert spaces and unitary maps. As a par- 
ticular application of the tools developed, we will see that the Schrodinger and Heisenberg 
representation of the quantum harmonic oscillator are equivalent. 

1.5.1 Extension of linear bounded maps on total sets 

More often than not, the unitary maps that we will consider are difficult or tedious to 
describe in terms of their action on arbitrary vectors. A more suitable way then is to define 
the action of such maps on a particular set of vectors which is "total", i.e. such that it 
uniquely determines the action on an arbitrary vector. 
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Definition 1.8 

A linear subspace Mq of a Hilbert space 9i is dense in H if for every there exists a 

sequence^ e # such that - t, n \\ 2 = (% - ^ - ^ n ) -> 0. 

If we have defined a map Ao on a dense subspace Hq that is 

i. linear on H§, and 

ii. bounded (i.e. there is a constants such that || Ao^o 1 1 <^||£o|| for all £o £ #b), 

then the action of Aq can uniquely be extended to the whole of H by setting At,:— lim„ A^ n 
for any sequence (^„)„ in H§ with ||£„ — £|| — > (since the limit does not dependent on the 
sequence chosen by virtue of the boundedness). The extension A is automatically linear 
and bounded with the same constant: ||A£|| < K ||^| (£, e 

Definition 1.9 

A subset X of 9i is total if its linear span is dense in H . 

In a Hilbert space, the closure of the linear span of a subset X is X ±J -, where X L is the 
orthogonal complement of X which is defined as {rj e H : (r\\t,) = for all t, e x}. Con- 
sequently, 

Lemma 1.10 

X C H is total if only if for any r| e H 

(il|§)=0 (forallt,ex) => r| = 0. 

We will make extensively use of the following: 
Lemma 1.11 

LetX be a total subset of M and letAo be a map X — > M . 
Suppose that 

(i) Aq can be extended to a linear map on the linear span of X i.e. 

or, equivalently, £" =1 a, Ao(^i) = whenever £" =1 a,-^ = 0; 

(ii) Ao is bounded on the linear span of X i.e. there is a constant K > such that 

||I? =1 a,Aofe)ll < *||I*=i a&ll (» e N, a,- eC,^4 
Then Ao can be extended to a bounded linear map on all of M . 
Example 1.12 (Weyl operators and exponential vectors) 

Let X be a Hilbert space and suppose that there exists a Hilbert space T (3C) (a so called 
Fock space over %) that contains a total set £ {%_ ):={■£(/) : / e SC } sucn that 

(E(/)|Efe)> y(ac) =e^ (*) 

A set of vectors with property (*) will be called a set of exponential vectors. We will prove 
the existence of such a Hilbert space 7 later. 

We now show how a collection of Weyl operators (indexed by % ) can be defined on such a 

set of exponential vectors. 

For/e %_ set W(f) : E(3C) — > J (?c) via 

W(/) £ (g) :=e-</ls>-Hl/ll 2 £(/ + s) (g € E(JC)). 

Then 

IIW^lfa^^Olll^E^a^W^EfeOllW^E^)) 

= ^a i a ; -e-<^->-5ll/ll 2 e-</l^)-2ll/ll 2 e (/+Sil/+^) 

= ^ .a ( . a; . e -te.1/>-^ll/ll%-(/l^)-2ll/ll%ll/ll 2 +(/l^)+<sd/)+(«d^ 

= E <J a i ayeW = ||E?a*Efei)ll 2 - 
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Thus, W(f) is bounded (in fact isometric) on the linear span of the exponential vectors. 
In particular, || I? a,-£ fe)ll - implies \\W{f)\% a,-E (g t )]\\ = 0. By lemma (1.11), we 
therefore conclude that W(f) can be extended to a unitary operator on the whole of J (%.)■ 
A straightforward calculation shows that W(f)W(g) = e~ llm ^^W (f + g) on exponential 
vectors and hence on the whole of f (sc ) by continuity. 



1.5.2 Kolmogorov's dilation theorem for positive kernels 

A positive definite kernel on a set S is a map K : S x 5 — > C such that 

£ £ YjliKisj^i) > (n e N, . . . A«) e C" , ji, . . . ,5, e 5). 

'=1.7=1 

For example, if H is a Hilbert space with inner product (.|.), then K(t,,r[) :— (£|t|) is a 
positive definite kernel. Slightly more general: any map V : S — > # defines a positive 
definitekemelA'(s / ,j):=(V(s / )|V(j)>. Indeed, 

L"=ikj\K(sj,Si) = (e?=i W;)IE2=i WO) - IIILi WOII 2 >o. 

The following theorem shows that every positive definite kernel can be written this way. 
Theorem 1.13 

Let K be a positive definite kernel on a set S. Then, up to unitary equivalence, there exists 
a unique Hilbert space and a unique embedding V : S — > Hk such that 

(V(s')\V(s))=K(s',s) (s,s'GS) (1.30) 
{V(s) : s e 5} is total in . (1.31) 

Given a kernel ^ ona set S, a map V from S into some Hilbert space M satisfying (1.30) 
is called a Kolmogorov dilation of K. Such a dilation is called minimal if, in addition, (1.31) 
holds. For the proof we refer to [11, section 1.5, p. 30] 

The following examples show that many (constructions with) Hilbert spaces in quantum 
mechanics are Kolmogorov dilations of suitable kernels. 
Examples 1.14 

- Let S be a set and K(s',s) := 8/ s . Then Hk = £ 2 (S) and V maps S to the standard 
orthonormal basis of H defined by e s (t) := 8 4 , r , s, t e S. In particular, taking S = N leads 

to£ 2 (N). 

- Let S be the Cartesian product of two Hilbert spaces 3Q and 

•Ki x %2 = {(£1,^2) : £1 e Ki , \i e 0(2}, equipped with 
Ar(fti,fc),(TU,Ti 2 )) :=^i|fc>^-(Tli|Tl2> rf2 . 

Then = % <g> is the tensor product of % and SC2 and V (£1 , ^2) = £1 ® £2- 

- Let 5= 3C be a Hilbert space and/T (^,r() := (^|r)) 2 . Then^ = % ® s % is the symmetric 
tensor product of K and V (£) = £ (8) 

- Let 5 = ?C be a Hilbert space and /f(^,r|) := e™^ . Then is a positive definite kernel 
and 9(k = 7 {%.) is the (symmetric) Fock space over defined as 

and V (£) is the so-called exponential vector or coherent vector given by 
Exp© = l©X e M e « e M|p... 

- The last example will be of use later, when we want to view the damped harmonic oscil- 
lator as part of an energy conserving system. 

Let S = R and K %a (s,t) := e ~n\s-t\+Hs-t) with r| > 0. A minimal Kolmogorov dilation 
of this kernel is presented by 9l = L 2 (R) and V : t ^ v, e L 2 (R) defined as 
( e {T\+ia){x-t) ifx < f 

Vt(x) := < 
W 10 ifx>t. 
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The first kernel (8 4 y) above is clearly positive definite, and the last kernel (K^, a (s,t)) can 
be checked to be positive definite by verifying its dilation. It is perhaps less obvious that 
the other kernels are positive definite, but for this the following lemma helps: 

Lemma 1.15 (see A.14) 

Let Ki , K2 , Kt, , . . . : S x S — > C be positive definite kernels. 
Then 

i. If r\,r2 G [0, °°), then r\K\ + r^Ki is a positive definite kernel. 

ii. The point wise product Ki *K2 : (s',s) 1— » Ki(s',s)K2 (s' , s) is a positive definite kernel. 

Hi. if L5T= 1 K n {s' , s) exists for all (s 1 ,s) G S x S, then K : (s 1 , s) 1— > 1 K n (s',s) is a positive 
definite kernel. 

In particular, e K ^ s '^ is a positive definite kernel is if K is one. 

1.5.3 Cyclic vacuum representations of the Canonical Commutation 
Relations (CCR) 

As promised we would show that the Schrodinger and Heisenberg representation of the 
harmonic oscillator are equivalent. The deeper reason for this equivalence is that both are 
so-called cyclic vacuum representations of the canonical commutation relations over C, and 
are therefore automatically unitarily equivalent. As it is no extra cost, we prove the latter 
in a much more general setting: that of Weyl representations of the canonical commutation 
relations. 

As introduction, recall from (1.27) that Weyl's form of the canonical commutation relation 
of one harmonic oscillator can be represented with Weyl operators W(z)" =" exp(;'Im {z)Q — 
/Re (z)P) that satisfy the Weyl relation 

W(z)W(w) =e- ilm ^W(z + w) (z,weC). 

In terms of Heisenberg's commutation relation, two independent harmonic quantum oscil- 
lator systems are described by two canonical pairs (Qi,Pi) and (22,^2) satisfying [Qi,Pj] = 
ih5 ij ,[Q i ,Qj} = [P i ,Pj}=0. 
The same is accomplished by 

1 . describing the canonical pairs (Qi,Pi) and (Q2 ,Pi) by the corresponding two families 
of Weyl operators, W\ (zi ),zi G C , and ^2(22), Z2 G C and, 

2. capturing the independence ([21,62] = [61,^2] = [62,^1] = 0) by taking the tensor 
product of the two Weyl operators W\ (z\ ) <X> W^fe). Z\,Z2 G C, which are unitary op- 
erators on H\®tt2- 

These tensored Weyl operators satisfy 

(Wi (zi) <8> 1*2(22)) (Wi (wi) ® W 2 (w 2 )) = Wi (zi)Wi (wi) ® W 2 (z2)W 2 {w 2 ) 
= (e" iIm( " wl) Wi(zi+wi)) ® (e- ilm ^W l {z2 + w 2 )) 

-ilm I (zi,Z2)(wi,W2) ) 
= e V / W\{Z\ +W\)®W2{Z2 + W2). 

Reformulating in terms of W^(zi,Z2)) := ^i(zi) ® 1^2(22). (21,22) G C 2 and taking the 
usual (coordinate wise) addition and inner product on C 2 , we have 

W((Z,,Z2))W((VV1,W 2 ))=^^ 

for (zi,Z2), (wi,w 2 ) G C 2 . 

This motivates to say that a representation of the CCR of two independent harmonic oscil- 
lators "is" a single representation CCR over Ki := C 2 in the following sense: 
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Definition 1.16 

Let %_ be a complex Hilbert space. A representation of the Canonical Commutation Rela- 
tions (CCR) over %. is a map W from % to the unitary operators on some Hilbert space H 
such that 

W(f)W(g) =e- ilm ^W(f + g) (f,g e at), 

and 

t^W(tf% is continuous {j E %_ ^E tt). 

The map W is called a vacuum representation of the CCR if, in addition, there is a unit 
vector Q. E H such that 

(i2|W(/)|i2)= e -3ll/ll 2 (/€*). 

A vacuum representation of the CCR is called cyclic if {W (f)£l : f E 0^} is total in 9i . 

Existence and uniqueness of cyclic vacuum representations of a CCR is provided by: 

Lemma 1.17 (Existence and uniqueness of cyclic vacuum representations of the CCR) 

Let %_ be complex Hilbert space. 

Then there exists a Hilbert space H and a map W from 5£ into the unitary operators 11 {H ) 
on H such that 

(i) (Strong continuity) 

1 1-» W(t f% is continuous for each f E%_ and ^ e M , 

(ii) (Weyl rela tion ) 

W(f)W(g)=e- ilm ^W{f + g)forallf,gE %L, 

(Hi) (Cyclic vacuum vector) 

There exists a vector Q. E M that represents a vacuum state i.e. (£l\W(f)\£l) = 
e~ i nJ II and that is cyclic i.e. the set {W (/)£2 : / E %_ } is total in H . 

Moreover, if (W' ,£l') is another cyclic vacuum representation of the CCR overKi, then 
there exists a unitary map U : X -» !H' such that W'(f) = UW(f )U~ l (f E %_). 

(For the proof see [11, p. 49]). 
Remarks 1.18 

- As an application of the above, we see that the two-fold tensor product of the CCR over 
C ((Wi,9(i,£li), i = 1,2) is unitarily equivalent to a representation of the CCR over C 2 

(tf,W(zi,Z2),Q)): 

(Mi X 2 , Wi (z i ) ® W 2 (ii ) , n i ® n 2 ) ~ {X , W (z i , z 2 ) , ii) ) 

- A slight generalization of the above yields that a representation of CCR of n independent 
harmonic oscillators (which is an n-fold tensor product of the representation of the CCR 
over C), is unitarily equivalent to a representation of the CCR over C" = £ 2 ({ 1 , . . . , n}) 

- In this sense, the dimension of the Hilbert space 9(_ describes the number of independent 
oscillators. In fact, if S is a set (measure space) of certain cardinality, then %_ = L 2 (S) 
describes a system of as many independent oscillators as the cardinality of S is. 

For S = N (so 3C = ^ 2 (N)), we obtain a countable chain of oscillators, while for S = R 
(so 3C = L 2 (M)) we get a continuum of oscillators. 

- Representations that are not cyclic vacuum representation need not be unitarily equiva- 
lent, for example those associated with positive temperatures [4]. 

We will now discuss two generic methods to obtain cyclic vacuum representations: via a 
Fock space functor or via a Gaussian functor. 
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1.5.4 Cyclic vacuum representations of CCR via the Fock space func- 
tor 

Let Ki be a complex Hilbert space. By Kolmogorov's dilation theorem, we obtain the Fock 
space J (Ki) over Ki. The Fock space contains a total set, viz. that of the exponential vec- 
tors £ (/), / G Ki and these have the property that (£ (/)|£ (g)) = e^ g \ see the examples 
in section (1.14). Weyl operators on J (Ki) can be defined via 

W(f)<E(g) :=e-(f^W 2 t(f + g) (f,g €£(*))• 

(see 1.12). Since £1 :— £ (0) satisfies 

(Ci\W(f)\Ci) = (£(0)| e -Hl/ll 2 £(0)} = e -lll/H 2 e <0|0) =e -lll/ll 2 , 

we can view Q. as a vacuum vector. Furthermore, as W(f)£l — e~^^ <E (f) and the expo- 
nential vectors are dense in f (Ki), Q. is a cyclic vector. 

A contraction C : Ki — > 3C induces a contraction J (C) : J (Ki) — > J (X) that maps £ (/) 
to £ (Cf) for every / g Ki (Lemma 1.11). Viewing the Fock space as 

this map can be written as 

j(c) = ®: =0 c® n . 

Given an orthogonal projection P on Ki , we can define (A. 20) a self-adjoint operator df (P) 

on T (Ki) via 

e^ 7 ^ ■=T(e {kF ) (k € R). (1.32) 

This operator J J (P) is interpreted as the quantum random variable that counts how many 
particles the "question" P is answered "yes". For example, the total number of particles N 
equals df (I). 

The above procedure derives its strength from its generality, but remains rather abstract. 
Some concrete representations of Fock spaces provide some flesh and blood to it. 

The Hilbert space ^ 2 (N) as a representation of the Fock space over C 

As C is one-dimensional with orthonormal base { 1 }, an orthonormal base of C ® C is given 
by {101} and as this basis (= one vector) is already symmetrical, it is a basis of C ® s C 
too. In other words, C <g> C and C ® s C are one-dimensional and hence unitary equivalent 
to C. 

Therefore, 

f(C) = c©c©c 8 >ec®'©...~c©c®ce... 

= {(Zn)nm ■ Z n € C, |z„| 2 < °°} 

= £ 2 (E). 

Otherwise stated, the map Jt : C — > £ 2 (N), 
z^n(z) := (i'Z'TfPTJp---)' 

defines a minimal Kolmogorov dilation of the kernel K(z,w) :—e zw , since (k(z)\k(w)) — 
e 7w and the set {n(z) : z G C} is dense in £ 2 (N) (lemma 1.7). 

In particular, the Heisenberg representation of the harmonic oscillator is (a Fock represen- 
tation of) the canonical commutation relations over C. 
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The L 2 over Guichardet's space (A(R) 7j u) as a representation of the Fock space over 

L 2 (R) 

The Fock space over L (R) is formally 

<J (L 2 (R)) = C®L 2 (R) ® (l 2 (R)®. s L 2 (R)) ffi (l 2 (R) ®. v L 2 (R) ®. s L 2 (R)) 0... 

Note that L 2 (R) ®. S L 2 (R) ~ L 2 ymm (M. 2 ) i.e. there is a natural correspondence between L- 
functions on R 2 that are symmetrical: 

f(h,t 2 )= f{t2,h) (for almost all fi,f 2 e R) 
and symmetrical tensor products of L 2 -functions on R: 

>8 + 8®f\,. ^ f{h)g{h)+g{h)f{h) 



-){h,t 2 ) = 



since the "elementary" symmetrical functions i(/ ® g + g ® /), /,g G L 2 (R), are dense in 
the Hilbert space of symmetrical L 2 -functions on R 2 (A. 13). In the same vein, 

L 2 (R)®" ~L 2 ymm (W) = {/ G L 2 (R") : f(h,. . . ,t n ) =f(t x{1) ,. . .,t z{n) ) for all T e <>„}. 

Due to its symmetry, a symmetrical function on R" is determined by it values on { (t\ , . . . , t„) G 
R" : t\ < t2 < ■ ■ ■ < t n }. Hence, we can identify a symmetrical function on R" with a func- 
tion on A„(R) := {o C R : #(o) = «} i.e. A„(R) is the collection of all subsets of R that 
contain n (distinct) points. 
Using the correspondence 

(h,...,t n ) with t\ < ti < ■ ■ ■ < t„ < — ► the set of distinct points {t\,...,t n } 

we can transfer the Lebesgue measure X n on W to A„. 
Let us introduce the Guichardet space over R 

A(R) = {o G R : #(o) < °c} 

= U A »( R ) 

where 

A (R) := {0} 

A„(R) :={ocR:#(o) =«} (n > 1), 

i.e. A(R) is the collection of all finite subsets of R, and equip A(R) with a measure /j given 
by 

/!({«}) = 1 

V U„(R) = Ia„(R) • 

By the above, L 2 ymm (R") ~ L 2 (A„(R)), and therefore: 

^T(L 2 (R))=®: =0 L 2 (Rf" 

— ©n=0 ^symm(^") 

~C =0 L 2 (A„ 
cL 2 (u: =0 A„ 

= L 2 (A(R)). 

Under the identification 0~ =o L 2 (R)®" ~ L 2 (A(R)), the exponential vector ®" =0 f®"/Vn\ 
corresponds to the "product" vector %(f) : A(R) — > C given by 

n(f)(a):=U sea f(s) with 7t(/)(0):=l. 

Then 
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W)\n(g))=f n(f)(oMg)(c)v(dc) 

JAR) 



°° r 

„— (\J A„(R) 



„ = J.si<.v 2 <...<i„ 



= e (/l«). 



In particular, the vacuum state H is represented as the function that is 1 on {0} and 
elsewhere. 

The interpretation of an element of A„(R)) is a configuration of positions {s\, . . . ,s n } of n 
indistinguishable particles. 

1.5.5 Cyclic vacuum representations of the CCR via Gaussian func- 



We will provide the Gaussian representation of the harmonic oscillator, but will seize the 
opportunity to present it as a simple (sometimes degenerate) case of a more general con- 
struction which will be discussed in the next section. This general construction will be of 
use in the dilation of the harmonic oscillator which will be needed further on. We illustrate 
the finite dimensional case with E = E and the infinite dimensional case with E = Lj^(E). 

Finite dimensional case illustrated with E = E 

To emphasize the similarity in the construction for the finite dimensional case with the 
infinite dimensional case, we will treat the case E — E along the same lines as the infinite 
case, which may look at first sight a bit pedantic. 

Step 1: the basic ingredients We start with E = E as a (one dimensional) real vector 
space and its dual E', which is identified with E again: each linear continuous functional co : 
E — > E is of the form co„ (r) = r ■ u for some u £ E. In the general case (discussed later on), 
there is a Hilbert space H in between £ and E', E C H C E 1 . In this case, H — L|({1}) ~ E 
and now coincides with E — E as well as with E' ~ E. In particular, || r\\ 2 H = r 2 (r E E C H). 
The bilinear form linking E with E' is denoted by (-|-) : E x E' — > E, (r\u) = r-u. Second, 
we need a measure P on the dual E' ~ E such that J E ie'^ u ' ) d¥{u) — e^ilMltf, (r e E). 
From (A.2), it follows that the Gaussian measure P on E' = E with density 



tors 



y(x) =e"2^/V2i 



(1.33) 



with respect to the Borel-Lebesgue measure does the job. 
Summarizing we have 




(1.34) 



where it should be born in mind that 

1 . we are dealing with a triple of spaces E C H C E' which are here all identifiable with 
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2. the integration is over u e E' ~ R; the o-algebra on E' is that which makes all point 
evaluations 

P(r) : u (r\u) = r ■ u (u £ E') , 

r e E, measurable functions from E' — > R i.e. it is the Borel(-Lebesgue) o-algebra 
on E' ~ R; the measure P on £" is P(B) = J B e-i" 2 -^ (B a Borel-Lebesgue set in 

£" ~ R); 

3. the term e ,ru in the integrand is actually e'^"\ where (-|-) is the bilinear form linking 
E with E'. 

4. the exponential e~2 r2 is actually e^llMliR where ||r||g := r 2 is the Hilbert norm on 
H = R. 



Step 2: an isometry (3 : # — > Z^(£',P) By taking r = in (1.34), it is clear that P is a 
Gaussian probability measure on E', and (1.34) can be read as: 

f e^d¥ = e~^ (ret), 
Je> 

i.e. p(r) is a Gaussian random variable on E' with zero expectation and variance r 2 . In other 
words, ||P(r)|U (in = Var(p(/)) = r 2 = ||r|| 2 , and $ : E = H ^ Li(E') is an isometry. 



Step 3: complexification of P to pc The complexification pc of p is simply defined as 
Pc( r + is) = P(r) + 'P(^). The isometric character of p passes on to Pc: 

(Mr + is)\Mr + is)) L 2 cm = (PW + ipWIPcW + ipW)^^) 

= (PMIPW)4( £0 + (PWIPW)4 (£0 = '- 2 +^ 

= (r+iy|r+iy) c 
From now on, we will refer to P c with p. 



Step 4: exponential vectors For z = r + is e C, we define the exponential vector e(z) := 
e P(^)-2^ 2 Xo calculate the inner product of two such vectors, we consider 



,P(z)uP(w) 



4(£') 



,P(r+is)uP(p+i?)\ 



/4(£') 

gPM-B'PW| e P(p)+«P(<7)\ =E ( e P(r+p)+iP(g- S )\ 



Using that p is an isometry R — > L^(E') and using (A. 3), the latter term is seen to be 



exp 



\{{r + p) + i{q-s)} 7 



exp 



i{(r + ^) + (/ ? + /< 7 )} : 



exp 



j{z + w} 2 = exp[±z 2 + jw 2 +zw] , 



i.e. (e(z)|e(w)) L 2 (£ , ) = e ? 



Step 5: the exponential vectors form a total subset of h\ (E r ) Actually, by the unique- 
ness of the Fourier-transform, the subset of all e'PM : u i— > e'"', r e £ = R, is already total 
mL 2 c (E'). 
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The Gaussian representation of the CCR over C The representation obtained above is 
called the Gaussian representation of the CCR over C. 

In this representation, the action of the Weyl operator takes a rather simple form as 



W(t)e(z) (x) 



-m,-m 



e~2i'i ef+z 



e(x-2t) 



y(x-2t) 



y{x) 



1/2 



(x e 



where y is the density of the Gaussian measure P on R, see equation (1.33). Further, 



W(is)e(z) (x) = 



e-i e(is + z) (x) = e w ~ i z e ,sx (x e 



Concluding: apart from a mixing in of a factor y/y(x — 2t)/y(x),W(t) is a shift to the right, 
and W(is) is a multiplication with e lsx . Things take on a familiar form when dividing away 
a factor \fyjx): 



Equivalence of the Schrodinger and Heisenberg representation of the CCR over C 

The linear map U GS : Ll(R,P) -> L|(R) defined by 



(U G sf)(x) 



(x e - 



is isometric and surjective, hence unitary. Using this map to transfer the Weyl operators to 
Lf.(R), we see them in the following form: 



[W(t)f](x)=f(x-2t) 
[W(is)f](x)=e is *f(x) 



(x e - 



i.e. we obtain the Schrodinger representation. 

In conclusion, we have seen that the Heisenberg representation coincides with the Fock rep- 
resentation of the CCR over C which is unitarily equivalent with the Gaussian representa- 
tion of the CCR over C, as both are cyclic vacuum representations. Further the Schrodinger 
representation is unitarily equivalent with the Gaussian, and so with the Heisenberg's. 



Infinite dimensional case illustrated with E = Lr(R) 

Step 1: basic ingredients We illustrate the infinite dimensional case of the Gaussian 
functor with E = L^(R). To start with, we shortly stating the ingredients, and, after that, 
we will provide their definition, comment and illustrate. The basic ingredients are: 

1. a triple of spaces E C L\{T) =: H C E', and 

2. a measure P on E' such that 

f e , Vl to >JP(co)=e-3ll/lll (feE). (1.35) 
Je' 



We now give the details: 

- the middle space is a real Hilbert space L^(T) where T is a countable set (T C N) 
equipped with the counting measure (so that L^(T) = £ 2 (T)) or a continuum (rcl 
finite or infinite interval) equipped with the Borel-Lebesgue measure; 

- T can be interpreted as time; 

- E is a so-called countably Hilbert nuclear space. The details of the defining properties 
of this can be found in [8, Appendix], but for our concern, the following examples (to be 
used further on) will suffice: 

* E = H = E' = R = L^({1}) which we saw before; 

* E = S, the real- valued test functions of Schwartz on R, H = Ljjj(R), and E' = s' is the 
space of tempered real-valued distribution (also known as generalized functions); 
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- The existence of a measure as mentioned in (1.35) (under the appropriate conditions) is 
asserted by the Bochner-Minlos theorem (Hida); 

- The space E is dense in L\(T) with respect to || ll^z (r)5 

- E' is the dual of E 

- the bilinear form E x E' — > R that sends (/,©) to co(/) is denoted by (-|-) i.e. (/|co) := 
©(/) (The reason for using this order of E' and E in the bilinear form will become clear 
later); 

- the injections E — > H and H — > E' are continuous; ([8, Appendix]). 

- Note that, since / is real-valued, ||/||| = J T f(t) 2 dt. 

The elements of E induce point-evaluations on E'\ P(/), / G E, where P(/) is the map 
E' 9 CO ^ ©(/) := (/|co) e K. 

On we have the weak-* topology, a(E',E), i.e. the smallest topology for which all P(/), 
/ G £, are continuous. Let J be the o-algebra generated by the weak-* topology o(E',E), 
otherwise stated f is the o-algebra generated by the the cylinder sets 

A fu ...j n .B ■= {co e E> : ((/! | oo) , . . . , (/») G B} , 

where n>l,/i,...,/»e£ and B is a Borel set in R". 

Step 2: an isometry p : Lj|(r) -> i|(£',P) By virtue of (1.35), P(£") = 1 (take / = 0), 
so (Q := £',5 ) is a probability space. For / G E, p(/) : £2 = £' — > K can therefore be 
viewed as a random variable on this probability space. Moreover, (1.35) implies that the 
characteristic function of p(/) is 

/ £ ,e'W)j P ( a) ) =e -^ 2 ll/lll ((gR), 

i.e. P(/) is a centered (E(p(/)) = 0) Gaussian random variable with variance = 
f T f(t) 2 dt. Reformulating: 

HPCOH^r) = j £l ViffdF = E(P(/) 2 ) = E(P(/) 2 ) - [E(P(/))] 2 
= Var(P( / )) = ||/|| 2 |(r) , 

which implies that p is a linear isometry E — > E\{E'). Since £ is a dense subset of H = 
Lj^r), P can be extended to a linear isometry between the Hilbert spaces H = L\{T) 
and L\{E'), also denoted by p. As such, p automatically preserves the inner product. In 
particular 

Cov(P(/),P(g)) =E(P(/)P(g))-E(P(/))E(P(g)) =E(P(/)P(g)) 
HP(/)IP(^) H/I^ ( r) 

Step 3: complexiflcation of the map P to Pc : L^,(T) — * ^cC^') We d enne tne complex- 
ification of p as: 

Pc(/+fe) = P(/) + ipfe) (/,*e4(r)), 

That pc : L^(T) — > L^(E') is an isometry, can be seen from 

(Pc(/ + ^)IP C (/ + ^)) 4(£0 - (P(/) + 'Pfe)IP(/) + 'Pte)>4( £0 

= (P(/)|P(/))+7/(Pfe)|Pfe))-/(Pfe)|P(/))+/(P(/)|Pfe)) 
= Jr / 2 + JrS 2 - * /r */ + ' Jr /* = /r / 2 + It 8 2 

= (f +W '+ig)q,(T) ■ 

From now on we drop also the subscript C from Pc. 
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Step 4: Exponential vectors For f c E L^(T), we define the exponential vector e(/ c ) as 

e(f c ) = eW°')-*h#, (1.36) 
and we calculate the inner product of two such vectors. 

Suppose/,. :=f\+if 2 E andg c := g i +ig 2 EL 2 z (T) with f 1 ,f 2 , gi,g 2 E £r( 7 ')- Then 

^ e P(/c)|gP(gc)^ = ^ e P(/l+V2)| e P(gl+ig2)^ ^ 

= / e P(/l)+l'P(/2)|gP(gl)+iP(g2)\ = EeP(/l+Sl)+''P(S2-/2) 

\ 

which is of the form Ee' x+Y withX := $(g 2 — f 2 ) andF = P(/i +gi) centered real Gaussian 
random variables. Since p is an isometry L^{T) — ► Ljg(Zs'), we have that 

- Var(X) + Var(F ) + 2;Cov(X , Y) 

= -U T (g2-f2) 2 dt] + [f T (fi + gi) 2 dt}+2i[f T (g 2 -f 2 ) ■ {fx +gx)dt] 
= It [(/i + + '"(ft - f2)] 2 dt = J T [(fT+if2) + (8i+ igi)] 2 dt 
= JtUc + gc} 2 dt = J T J 2 dt + J T g 2 c dt + 2J T Jcgcdt 

so that by (A.3) 

E^Pfe-^+PW+sO = exp [i !t J 2 + i + /r ^ c 

Combining the above yields, 

m)K8c)) L 2 c{EI) =e fcl8c)L ^ (f c ,g c EL 2 c (T)), (1.37) 
whence the name exponential vectors. 

Step 5: Exponential vectors form a total subset of L}.{E') For f c = fi + if 2 E L^(T), 
the random variable |eP^ c )| 2 = e 2 P(/i) is of the form e Y where Y = P(2/i) is a centered 
Gaussian with variance 4 J T f\ (t) 2 dt. From (A.3), it therefore follows that 

e(/ c ) = eMM-ihMt) 2 * e L 2 (£')• 

Next, the set of exponential vectors is also total in Lq(E'), i.e. the linear span of the ex- 
ponential vectors is a dense linear subspace of Lg(E'). Actually, it suffices to consider 
real-valued test functions: 
Lemma 1.19 

The set {e'^ : g E E} is a total subset in L\(E'). 

<3(See [8, p. 116]). Note that £ is a dense linear subspace of L^(T). Now, let A be the complex 
linear subspace spanned by the e'^ g \ g e E. Since e ! P(«i) • e'Pfe) = g'Ptei+ft^ ^ j s actually an 
algebra. 

To prove that A is dense in Lg{E') is suffices to prove that if E L^E 1 ) is orthogonal to all 
random variables of the form 

exp[iP(E; *,■«)] , (f, GR, gi EE), 

then <\> = P-almost everywhere (P-almost surely). 

Step 1 

Fix n E N, gy , . . . ,g„ EE, and let T n E T be the o-algebra generated by X\ := P(gi), ...,X n := 
(3(g„). Suppose we have for every (t\,...,t n ) E R": 

f E ,cxp[ir i itkX k ]^dP = 0. 

Since E' E 7 n , the second equality in the following chain is valid: 

= J E > ^P [I" it k X k ] $ df = / £ , E (exp [£? ifcX*] | y„) dP 

= / F exp[l5f^]E($|^)rfP. 



1.6. Functorial quantization of the harmonic oscillator 



33 



where we have used in the third equality that exp [£" %XjJ is J„-measurable. Since a function is 
J n -measurable if and only if it is a function of (Xi, . . . ,X„) we can write E((j) | J n ) = f(Xi ,...,X„) 
for some measurable / : W — » K and conclude that 

i.e. the n-dimensional Fourier-transform of f(x\ , . . . ,x n ) with respect to the joint-probability dis- 
tribution of (X\ ,.. . ,X n ) vanishes on R". This means that f(x\ , . . . ,x„) = for P( Xl x„)~ amlost 
all (x u ... ,x n ) el" i.e. P[/(Xj , . . . ,X n ) / 0] = P (Xl) ...,x„) [/ 7^ 0] = 0. 
Conclusion: E(0|j„) = OP-a.e. 
Step 2 

Choose a system {gi,g2,---} in E that is complete in L^(T). Then the corresponding (J„ : = 
a{gi,...,g„} increase to f (i.e. Q n C (?„+i C 7 and J is the smallest o-algebra containing 

<l Recall that J is the smallest a-algebra such that all P(/), / e £, are measurable. 
Let f <E E. We prove that / is measurable with respect to the o-algebra generated by l) n g n , so 
that J is also generated by U„ Q n . 

Indeed, since E C L^(T) and the linear span of {g n : n <E N} is dense in E^(T), there exists a 
sequence /„ = Y?" a nigni of finite linear combinations of {g„ : n € N} such that ||/ — /„|| 2 — ► 0. 
In particular, f n (t) — > f(t) for almost all t € T, and / is measurable with respect to the o-algebra 
generated by {/,:»£ N}. Since each /„ is g^ n measurable, / is measurable with respect to the 
a-algebra generated by U n g„. E> 
Let Q := {A <E f : j A §d¥ = 0}. Then g is a o-algebra that contains each g n (by step 1 above), 
whence i.e. J A $dP = for all A e f ■ This implies that (|) = P-a.e. > 

1.6 Functorial quantization of the harmonic oscillator 

Instead of Dirac's quantization procedure, we can now outline a functorial quantization 
procedure that acts on the classical evolution in phase space. We take the simple classical 
harmonic oscillator as illustrating example. 

Step 1: Identify the classical phase space with a Hilbert space The simple harmonic 
oscillator has a position q and momentum p. Its phase space R 2 can be identified with C 
via (p, q) <-> z = p + iq- If we have n oscillators, the phase space corresponds with C" = 
£ 2 ({1,. . . ,«}). For a chain of oscillators, we have as phase space £ 2 ({1,2, ...})= ^ 2 (N), 
and a continuum of oscillators labeled by x € JR may correspond to L 2 (JR). 

Step 2: Identify the evolution in phase space with an evolution in the corresponding 
Hilbert space As the explicit time evolution of the harmonic oscillator is given by 
p(t) = -q(0) sin((0f ) + p(0)cos(cof) , 

q(t)= #(0)cos(cof) +p(0) sin(aw) , 
the evolution in terms of z reads 

z{t) = e iw z{0), 

i.e. the time evolution operator T t is the multiplication by e' m . 

In the same vein, n harmonic oscillators, each with frequency co,, can be described in time 
by the evolution (zi,...,z„) ^ (e iai 'zi, . . . ,e im "'z n ) in C". 

Step 3: Quantization of the phase space Next we need to specify how to obtain the 
canonical observables (Q,P), and in particular the Hilbert space they act on. 
This is done by taking a cyclic vacuum representation of the CCR over C. For example, one 
could take as Hilbert space the Fock space over C and end up with the Heisenberg represen- 
tation of the harmonic oscillator. Another, unitarily equivalent, choice would be to take a 
Gaussian representation which is unitarily equivalent with the Schrodinger representation. 
For n oscillators, one would take the CCR over C" = . . . ,«}) and for a chain of 

oscillators the CCR over £ 2 (N) is appropriate. 
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e mt z, can 



(re 



Step 4: Quantization of the time evolution The classical time evolution z 
now be quantized by describing a time evolution of the corresponding Weyl operators 
W(z) i— > W(e m z). In terms of the canonical observables Q and P, defined as the gener- 
ators of the strongly continuous unitary maps s W(is) and t \— > W(t) respectively, the 
evolution comes down to 

P(t) = -e(0)sin(cof)+P(0)cos(cof)," 
Q(t) = 2(0)cos(cof)+P(0)sin(cof), 

which corresponds to what we would expect based on the Heisenberg equation of motion 
for Q,P [17, 2.2.19, p. 83; 2.3.4, p. 95]. 

For n oscillators the evolution in C" described by W(zi, . . . ,Z n ) i-> W(e mt z\, ■ ■ .,e l(0nt z n ) 
reads in terms of corresponding infinitesimal generators 

Pi(t) = -Q^(0)sm(ai^t)+Pi(0)cos(<alt), , 

Qi{t) = 2,(0) cos (co,f) +P i (0)sin(co 1 f) , t 

In the remainder of this section we will discuss observables, states and wave-particle duality 
in the functorial quantization. 



(t 6i, i= l,...,n). 



1.6.1 Wave-particle duality in the harmonic oscillator 

Applying the functorial procedure outlined above yields the same quantum analogue of the 
classical harmonic oscillator as the usual Dirac procedure. Indeed, by choosing appropriate 
cyclic vacuum representations of the CCR, the familiar Schrodinger and Heisenberg repre- 
sentations are obtained respectively. States in the Schrodinger representation are described 
by (complex-valued) functions on M (= configuration space of position) that express the 
wave-like properties when measuring the position the harmonic oscillator, whence their 
name wave functions. Also, momentum states are described by wave functions. 
In the Heisenberg representation, states are (complex-valued) functions on N (= configura- 
tion space of energy) that relate to the particle-like properties found when measuring the 
the energy of the oscillator. 

Identifying the exponential states in Schrodinger representations with those in Heisenberg's 
representation, both representations are easily seen to be unitarily equivalent. 



1.6.2 Algebra of observables 

For a classical harmonic oscillator, the position and momentum observable q and p com- 
pletely determine the state of the system and therefore each observable o of the system can 
be viewed as a function o(q,p) of q and p. 

After quantization, q and p are represented as unbounded self-adjoint operators Q and P 
on a Hilbert space. A reasonable guess at the quantum observable O corresponding to 
o = o(p,q) can be made if o(p, q) admits a Taylor series expansion in q and p. In that case, 
a formal substituting q^>Q and p P in the Taylor series of o can be effectuated and after 
that, a domain for O on which O is well-defined and in addition self-adjoint can be looked 
for. 

From a mathematical viewpoint, the above way to obtain quantum analogues of the clas- 
sical observables is full of obstacles and pitfalls, since it is often hard to find a domain on 
which the Taylor series expansion can be proven to be self-adjoint. 

Therefore, Weyl took another approach using bounded operators, as the latter are well- 
defined on the whole Hilbert space and therefore mathematically easier to deal with. The 
main ideas are the following: 

i. Restrict to bounded observables of the system instead of unbounded ones. 

ii. If o is a bounded observable and therefore a function o(q,p) of q and p, then use 
its Fourier transform o(s,t) that satisfies o(q,p) = j fo(s,t)e lsg ~" p dsdt to define the 
quantum observable corresponding to the classical observable o as 
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= 0(Q,P) = [°° f d{ Sl t)e isQ - itP dsdt. 

J —oo J —oo 

We comment on the above ideas. 

The first idea is no restriction as the properties of unbounded self-adjoint operators can be 
encoded in bounded operators. For example, the event that a measurement of an observable 
O yields a value in (a,b) (i.e. the observed value o of O has a < o < b) is described by 
a projection operator Pi a ,b)- One can therefore imagine that the family of bounded self- 
adjoint operators {P( a ,y) : a < b in R} describes the observable O completely. 
A mathematical description of that is given by Von Neumann's spectral theorem. Before 
formulating that, let us describe its ingredients: 

Definition 1.20 

A projection-valued measure (better known as spectral measure in the physics literature) 
P is a map from the Borel o-algebra Borel(M.) to the algebra of bounded operators on a 
Hilbert space M , <B (H ), such that 

i. P(B) is a projection for each Borel set Bel, 

ii. P(B l nB 2 ) = P(A)-P(B) for Borel sets B h B 2 CR 
Hi. if Bi ,B 2 , ■ ■ ■ are disjoint Borel sets in R, then 

]im N ^Z l »P(B n )Z, = P(U n B n )Z ) ($e#). 

For each ^ e M , the map P^\ B ^ (^\P(B)\E,) is a measure on the Borel sets on R. If 
J* R \x\ 2 P^(dx) < oo then the integral J R xP^(dx) ~: f R x (^\P(dx)\^) exists (in M), and the set 
{£, : J R \x\ 2 (^\P(dx)\E,) < oo} is a linear subspace of . 

Theorem 1.21 (Von Neumann) 

There is a one-to-one correspondence between self-adjoint operators O on a Hilbert space 
9i and projection-valued measures P : Borel(R) — ► ^(M). 

Specifically, if O is a self-adjoint operator on 9i , then there is a projection valued measure 
P on •}{ such that 

i. the linear subspace Hp = e H : / R \x\ 2 (QP(dx) \%) < oo } j' s dense in H ; 

ii. O is well-defined on Hp and (£,\0\£,) = J R xP^(dx) for all \ e H P . 

Inparticular, ||0^|| 2 = (0^0$ = J R \x\ 2 ($|P(<fr)|$>. 
Symbolically, we denote this as 

= f R xP(dx). 

Conversely, given a projection-valued measure P on H and T) G H , we have that the se- 
Quencex\fi '.— P(-n,«)r| converges to ^, and is contained in the linear subspace Hp := e 
H : / R |x| 2 < oo } (A. 15). Hence, the latter subspace is dense in H . A self-adjoint 

operator O can be defined on Hp via the formula := f R x(^\P(dx)\£,) (and (%\0\r\) 

is from there defined by polarization) (see A. 15). 

Example 1.22 

In the Schrodinger cyclic vacuum representation of the CCR over C, W (is) = e' s @ is 
the multiplication operator \W(is)f\ — e' sq f(q) in L 2 (R) (see the end remarks of section 
1 .2.3). The form of the projection-valued measure of Q in the Schrodinger representation is 
P()(a,b) = l/ a m. This corresponds to the fact that in a normalized state \\t, the probability 

of measuring a value of Q that lies between a and b is / R l( a b -j (q)\ \y(q)\ 2 dq = \\\f(q) \ 2 dq. 
[17, p. 101]. 

Weyl's second idea, viz. defining quantum analogues via the Fourier transform, works if 
the bounded observable o — o(q 7 p) has a Fourier transform which is guaranteed if o is in- 
tegrable over the phase space (/ / \o(q,p) \dqdp < oo). Identifying as usual the phase space 
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with C (via (q,p) p + iq — z) and e ' s Q+ ltP w (t + is), we want the classical observable 
o(q,p) = j )o(s,t)e is i +it Pdsdt to correspond to 0(Q, P) = J fo(s,t)W(t + is)dtds. 
The latter integral symbolizes that some limit of (superpositions of) Weyl operators is taken, 
and we have to elucidate what kind of convergence is involved. It turns out that so-called 
strong convergence enables us to obtain interesting operators from (superpositions of) Weyl 
operators. 

Definition 1.23 (Strong convergence) 

A sequence of operators (T n )™ =1 in <B (X) converges strongly to an operator T e 2 (#") if 
T„^ -^T%(n^ oo) for each \ e X (i.e. lim„ \\T n ^ - Tt\\ = for all % e X). 

Strong convergence implies weak convergence: 
Definition 1.24 (Weak convergence) 

A sequence of operators (T n )™ =1 in <B (X) converges weakly to an operator T G 'B (X) if 
^\T n \r\)^^\T\r\)(n^^)foreach^T\^9{. 

To obtain interesting (bounded) operators from the Weyl operators and to do all the opera- 
tions that are commonly used in quantum mechanics, we need to take linear combinations 
(superpositions), compositions, adjoints and strong limits. 

Therefore, the collection of all bounded observables of the quantum harmonic oscillator 
system is mathematically best described by a Von Neumann algebra: 

Definition 1.25 

A von Neumann algebra on a Hilbert space X is a collection X of bounded operators on 
X that is closed under taking 

i. superpositions (linear combinations): A,B e 5\£ , OC, peC=>aA + pBG9Y;, 

ii. composition: A,B e X => AB e X , 
Hi. adjoints: A e X. => ^ e 9i , 

iv. strong limits: 9i 3 A„—>A <E <B(H) strongly ^A<E9i. 

For simplicity we will always assume that a von Neumann algebra in <B (#) contains the 
identity of S {H ): I e X. ■ 

Let «o be a subset of$(X). Then the von Neumann algebra generated by the elements of 
So, i.e. the intersection of all von Neumann algebra that contain So, is denoted by vX ®o}- 

The interesting operators that we can obtain via strong convergence from the Weyl opera- 
tors of the harmonic oscillator include all bounded measurable functions of Q, P and of N, 
and therefore in particular the spectral measures (projection valued measures) associated 
with Q, P, and N. 

The mathematical formulation of this is as follows. 

Definition 1.26 ( bounded measurable functions of a self-adjoint operator) 

Let O be a self-adjoint operator with projection-valued measure Po, and let f : M — > R 
be a bounded Borel function. Then for \|/ € X , B i— » (\|/|Po(B)|\|/) is a finite measure and 
I-oof( x ) (ty\Po(dx)\\\f) < oo. Therefore, an bounded operator f{0) can be defined via 

(W(0)\V) = (MPo(dx)\y) ( V e X), 

since matrix elements (Vi|/(^)|V2) can be defined via polarization (A.4). If f is real- 
valued then f{0) is Hermitean. In particular, if f = t B for some Borel set B, then f{0) = 
Po(B) i-e. the projection-valued measure ofOis obtained. 

With the notion of functional calculus above, we can now state 
Theorem 1.27 

Let (X ,Q) be a cyclic vacuum representation of the CCR over C (so a cyclic vacuum 
representation of the quantum harmonic oscillator), 
i. Let Q be the position operator and Pq its spectral measure. Then 

vX {e isQ :seR} = vX{PQ(a,b):a<b in R} = vX {f(Q) ■ f step function} 
= vX{f(Q):feL°°}. 
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ii. Let P be the momentum operator and Pp its spectral measure. Then 

v9{, {e~ i,p : t el} = vM{P P (a,b) : a < b in R} = vOt {f(P) : f step function} 
= v9i{f{P):feL-} 

Hi. Let H = h(N + \ ) be the Hamiltonian and let Pn be the spectral measure of the number 
operator N. Then 

v*l {s* : s e (-1, 1)} = v*[{P N (n) : n e N} = {f(N) :/gf} 
Cv9i{W{z):zeC} 

In particular, the von Neumann algebra generated by the Weyl operators, <W = v5\£ {W(z) : 
z € C}, contains all bounded functions of Q, P, andN. 

A key lemma to prove the above, is Weyl's previously mentioned idea of writing an operator 
0(Q,P) =// o(s,t)W(t + is)dtds. 

Theorem 1.28 (Proof in A.16) 

Suppose: 

(i) p : C — > C is continuous, 

(ii) W : C — > S (y{ ) is a strongly continuous family in i b{M), 

(Hi) themapz^> |p(z)| ||W(z)|| fromCinto [0,<») is integrable with respect to the Lebesgue 
measure on C. 

Then, T e v^{W(z) : z € C} C fl(tf), and ||r|| B(fi) < / c |p(z)| ||W(z)|| X 2 (Jz). 

(A version in which continuity of p is dropped can be found in [20, chapter 3]). 

Using the above lemma, we can prove the follow lemma, which is a step up to the first two 
statements of theorem 1.27. 
Lemma 1.29 (Proof in A.17) 

Let (a,b) c K. Then there exists a sequence of functions fa such that each fk is 

i. infinitely differentiable; 

ii. rapidly decreasing in all its derivatives i.e. x i— > x n f^ is bounded for alln,m e N; 
Hi- 0<fk<l; 

iv. fk coincides with t( a ,b) on JR\ ([a — j,a] U [b,b + j}) 

(Note that the first two properties tell that the fk a test functions in the sense of Schwartz). 
Then the Fourier transform of each fk is (continuous and) integrable, and 

S-Jk{s)e- isx ds = f k (x) 
By theorem 1.28, the operators 

P k Q {a 1 b) = j fk{s)W(is)ds, and 

/oo 
fk{t)W(—t)ds (the minus sign occurs due to the minus sign in e lsQ ~" p ) 
-co 

are in v3\£ {W(z) ■ z £ C} and so are their strong limits 

PQ(a,b) =s- ]iraPp(a,b), 

Pp(a,b)=s-1imP£(a,b), 
where Pq and Pp are the spectral measures of Q and P respectively (see 1.27). 
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Using the above result we can complete the first two statements of theorem 1.27. 
Proof of theorem 1.27 iji 

We will prove the following sequence of inclusions for the operator Q (the proof for P is 
analogously): 

v*L{f{Q) : f e L°°(M)} D v9i{e isQ : s E R} D v*L{P Q {a,b) : a < b in R} 

3 v*C {/(g) : / step function } D {/(g) : / G L°°(M)}. 

The first inclusion is trivial, since {f(Q) : / e L°°(M)} D {e"'e : s e R}. The second is 
the purport of lemma 1.29 above. The third inclusion follows from a standard measure 
theoretic argument. 
OLet 

Q ={Ae Borel(R) : 1 A (Q) e v9l{P Q (a,b) :a<bmR}}. 

We will prove that Q = Borel(R) meaning that 1b(Q), for each Borel set B, and hence also linear 
combinations of them, i.e. step functions are in \>9l {Pg(a,fe) : a < b in R}. 
First observe that the intervals (a, b) C R are contained in g as li a ,b){Q) = ^g( a !^) m me func- 
tional calculus for self-adjoint operators: 

(v|l(„.*)(G)lv) = = (MP(2(dq)\v) = [ {v/\P Q (dq)\\f) 

= <Y|P fl (fl,fc)|V> (V^)- 
Next, observe that Q is o-algebra: 

i. R e g , since every von Neumann algebra contains I 6 « (jf). 

ii. A e ^ R\A e <y , since 1r\a(2) = 1r(2) - 1a (2) = I - 1a (2) is a linear combination 
of elements in vJ\£ {-Pgfa, b) : a < £> in R}. 

iii. A\,A2,... e <J disjoint => U„A„ e g, since in L°°(R) the functions fy := £^ =1 1a„ con- 
verge point wise (almost everywhere) to / = lu„A„- Therefore the multiplication operators 
Mf n converge strongly to Mf in «(L 2 (R)): by Lebesgue's dominated convergence theorem, 
||/ng — fg\\2 ~ ¥ for every g € L 2 (R). Thus in the Schrodinger (cyclic vacuum) representation 
of the CCR we see that lu„A„ — / — s ~ lim/^ fy G vfA£ {Pq(cl : b) : a < b in R}. 

Concluding, g is a o-algebra that contains the open intervals of R and must therefore contain 
Borel(R).D> 

The last inclusion follows from the fact that every / e L°°(R) can be uniformly approxi- 
mated by step functions: J\ : Y^=-N n / N \\f\\" ' ~k\fe[n/N\\f\\~,(n+i)/N\\f\\«,)]> so that /n —> 
f strongly in the Schrodinger (cyclic vacuum) representation of the CCR (Lebesgue's dom- 
inated convergence theorem). □ 

The following lemma underpins the last statement of theorem 1.27. 
Lemma 1.30 

We can write s N as a strong integral with respect to Weyl operators: 
1 r e~2 — w 

^ = -/-j w(z)Mdz) (*e(-i,i)). 

Furthermore, 



P n = (n!) 



—s N 

ds" 



s=0 

where the derivative can be taken in the strong operator topology. 
Proof of lemma 1.30 and theorem 1.27 iii. 

By theorem 1.28, the integral over the Weyl operators exists and a verification shows that 
it coincides with s N (see A. 18). We can write = YX=o^P{k) (s£ (—1, 1), where P(k) is 
the projection onto the linear subspace spanned by the eigenvector e\ of N with eigenvalue 
k (i.e. the fe-particles subspace). Since the series s k is absolutely summable on (say) [—3,3] 
and the collection operators {P^ : k e N} is bounded (in the uniform operator topology), we 
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can prove by induction that is infinitely differentiable (in the uniform operator topology 
A. 19, and hence) in the strong operator topology and has as n-th derivative 

=IZ =n k-(k-i)-(k-n + i)- s k -»-P(k) (*e(-±,±)). 

This shows that v5\£ {s N :se (0, 1)} contains all {Ppj(k), k e N, and, a fortiori all functions 
in £°°(N) with finite support: / = £f =1 a^(k). By taking strong limits of the latter, we 
obtain all functions in £°°(N). □ 

Resuming: the quantum harmonic oscillator is completely described by its canonical ob- 
servables (Q,P), or equivalently, given the family Weyl operators W(z) — e ls ®~ " , z — 
t + is e C. All bounded observables of the quantum harmonic oscillator are bounded func- 
tions of Q,P and we can obtain these by taking the von Neumann algebra generated by 

W(z)=e is <2- itP ,z = t + is€C. 



1.6.3 States of the harmonic oscillator 

In quantum mechanics states of a quantum system are often given by a normalized vector 
\\t in the Hilbert space H i.e. quantum theory postulates that all physical information of the 
system is encoded in \|/ as follows: the expected value of a harmonic oscillator observable 
A is given by 

E V (A) = ( ¥ |A| ¥ ), 

and, more informatively, the probability of measuring a value of A in the interval (a\,a2) 
is given by 

P v ((fli,fl 2 ) - (y\PA(aua 2 )\y) =E v (P A (ai,a 2 )) (1.38) 

where Pa is the projection-valued measure associated with A (which lies in the von Neu- 
mann algebra generated by A, whence in 'W = v9\[ {Q,P}. 

For Py to be interpretable as a probability measure (as postulated by quantum theory), we 
need that 

1. P V (M) = 1, 

2. P v (5) > for all Borel sets S in M, 

3. P ¥ (K \ S) = 1 - P v (5) for all Borel set S, 

4. If 51,52, Sj, . . . are disjoint Borel sets with union U n S„ = S, then P v (5) = Lr^y^")- 

This is taken care of, since (Py = Ey o P A , P A is a projection- valued measure, see definition 
(1.20), and) Ey has the following properties: 

1. Ey(A) is unity-preserving: Ey(l) = 1, 

2. Ey(A) is positive: Ey(A^A) > for all bounded observables A, 

3. Ey is additive: Ey(A +B) = Ey(A) +Ey(B), 

4. E ¥ (A) is normal i.e. if Ai < A 2 < ... < A and A„ — » A strongly, then Ey(A„) | E ¥ (A). 
Furthermore, Ey is not only additive, but in fact linear, since we have Ey(AA) = A,Ey(A) 

(keC). 

It is clear that Ey can be defined on the whole of but for describing the state of 

the oscillator observables, we might as well restrict it to the von Neumann algebra of all 
observables of the harmonic oscillator, i.e. W = v9\[ {Q,P}. 

The above example of a state, induced by a normalized vector in !H , is called a pure state. 
Pures states often arise after a "filtration" (selective measurement) i.e. an experiment in 
which the system is exhaustively measured and is described by a (simultaneous) eigenstate 
of a complete (system of) observable(s). Normally, the state of the system is not completely 
known and must be regarded as a convex combination of pure states: \|/ = r„\|/„ where r n 
are positive weights reflecting the relative amount of state \|/„: £„r„ = 1 ([1, 14.1]). Such 
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a state is called a mixture and described by a density matrix [17, p. 177] [1, 14.1] i.e. there 
is a positive operator p such that (\|/|A|\|/) = 7>(pA), where Tr denotes taking the trace. In 
any case, the states of the form Ey(A) = 7>(pA) also satisfy the conditions 1.-4. above. 
Thus, if we choose to describe the observables of a quantum system with a von Neumann 
algebra, then we are led to the following definition of a state: 

Definition 1.31 

A state on a von Neumann algebra 5\£ a positive, linear, identity-preserving functional 
(j) : A — > C. Itis called normalifO < A„ | A andA„ A strongly implies that§(A n ) | (j)(A). 

The normal states are those that associate probability measures with projection valued mea- 
sures of observables. 



Chapter 2 



The quantization functor T for an 
arbitrary number of oscillators 
and dissipative time evolutions 

The previous chapter has shown that the functorial quantization procedure above repro- 
duces well-known results when used to quantize the harmonic oscillator. This gives us 
faith to apply it in a wider context, in particular to the damped harmonic oscillator. To this 
end, this section will discuss the general framework of the quantization functor T. 

2.1 Classical input: a phase space with a time evolution 
described as a Hilbert space with a contractive map 

In principle the general quantization procedure will not be applicable to all conceivable 
classical systems, but only to those that can be described by a system of oscillators (possi- 
bly after a suitably choice of generalized coordinates). This restriction is not that stringent 
as it may seem at first sight, since many systems near stable equilibrium can be well ap- 
proximated by such an oscillator system. A further demand is that the evolution in phase 
space has to correspond to a linear contraction, in the sense that will become more clear 
below. 

Phase space 3C To prepare an oscillator system for quantization, its phase space must 
be framed as a (complex) Hilbert space %, with as dimension the number of independent 
oscillators i.e. the degrees of freedom of the system. Thus if the system has n degrees of 
freedom (possibly n = °°) then the standard choice for the Hilbert space is ( 2 {{l, ■ ■ ■ ,«} 
(i.e. C" if n < °° and ^ 2 (N) if n is countably infinite). The independent oscillators then 
correspond with orthogonal one-dimensional subspaces in this Hilbert space i.e. the phase 
space (qk,Pk) of the k' h oscillator is C • e\ where the vectors {ek}k form an orthogonal base 
of 9C . Observables are linear (continuous) functionals ac — > C and states are vectors in %_ . 
E.g. an harmonic oscillator evolution in oscillator k can be described by the time-dependent 
state L,k{t) = e~ m et- On the other hand, the position observable of the \i h oscillator is 
^Re« e *|Q). 

Energy 1 1 1 1 2 The square of the norm of the Hilbert space 3C is interpreted as the energy 
of the oscillator system (possibly after a suitable transformation of the Hamiltonian). This 
is most easily seen in the standard choice of Hilbert space: every oscillator {qk,Pk) con- 
tributes \{q\ + p\) to the total energy of the system, whence the square of the Hilbert norm 
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\ Hkill + Pk) corresponds to the total energy of the system. 

Time-evolution Q The most ideal time-evolution map in phase space is that of a energy 
conservative system, and this corresponds in the above framework (»C, || ||) with the most 
ideal map in a Hilbert space: a unitary (i.e. norm-preserving) map. As we will see further 
on, an advantage of the quantization functor F is that it can deal with contractive maps in 
%i as well which opens the door to quantize dissipative oscillator systems. 
Summarizing, a time-evolution in the above sense is a family of contractions C t : %_ —* Ki . 
Here t runs through an index set T which denotes the times; T can either be discrete (N, Z) 
or continuous ([0,°°), K). 

The above framework is summarized in the following definition 
Definition 2.1 

A classical dynamical system consists of a Hilbert space %_ equipped with a family of 
contractions Q, t e T. 

2.2 Quantum output: a quantum system described by a 
von Neumann algebra of observables and a time evo- 
lution of completely positive maps 

The quantization functor F is specified by describing what it does to a Hilbert space 3C and 
to a contraction C. 

r(ac ) is a von Neumann algebra W^- with vacuum state § x First we take a represen- 
tation of the Canonical Commutation Relation over SC • This gives us a representation of 
the Weyl operators W(f), f E X, over some Hilbert space H . The Weyl operators encode 
for the position and momentum observables of the system of oscillators, but to obtain the 
algebra of all observables of the system (at least those that are determined by the position 
and momentum operators of the system such as the spectral measures of the positions and 
momenta of the oscillators), we take the strong closure and obtain a von Neumann algebra 
'Wx . States of the system are given by normalized positive normal functionals : <W X — ► C 
that are unity-preserving i.e. 0(1) = 1. 

For simplicity (and because it suffices for the situation we want to describe), we will assume 
that our representation of the CCR(^) is a cyclic vacuum representation, i.e. there is a 
normalized vector £l x in the representation Hilbert space H such that 

Being a pure state, 0^ is induced by a vector Q.^ € 9{ , therefore naturally defined on 
the whole of 2? (jH ) D Also note that all cyclic vacuum representations of r(ac) = 
(W x ,0^) are unitarily equivalent (A.21). 

A physical interpretation of choosing a vacuum representation is that all states of the sys- 
tem can be obtained by successively applying ladder operators (annihilation and creation 
operators) associated with / £ %_ to the vacuum state. In other words, a general state of the 
system is made of excitations of each of the constituting oscillators (which are labeled by 

/esc). 

T(C) is a linear strongly continuous, completely positive, state- and unity-preserving 
map between von Neumann algebras of type (^5^,030) Up-to-now, the contraction 
C has had the interpretation of a time-evolution in the Hilbert phase space 3C . This time- 
evolution could either correspond to that of a energy-conservative system or a energy- 
dissipative system. Later on, we will like to consider the damped harmonic oscillator, 
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which is dissipative quantum system, as a subsystem of a larger, enveloping, energy- 
conservative system and investigate relations between the two systems. This means that 
we must be able to shuttle "up" from a subsystem to its envelop (via an isometric inclu- 
sion) and shuttle "down" from the envelope to the subsystem (via an orthogonal projection). 
To accommodate this, we need the general formalism of quantizing a contraction between 
two different Hilbert phase spaces: C : — ► %2- 

We will first motivate the mathematical properties of F(C), for C : 3Q — > %2 a contraction 
between two Hilbert spaces. In the next section we will define its action on the Weyl 
operators, leaving the proof of its extension to the appendix. To relieve the burden of 
parentheses in notation, we will write Fc := F(C) below. 

As we want that linear relations between observables and taking limits of them in the von 
Neumann algebra of observables of system 1 ( W Xl ) are reflected in that of system 2 ( W X2 ), 
we require linearity and strong continuity: 

A 3 = aiA 1 +a 2 A 2 => r c (A 3 ) = ccir c (Ai) + a 2 r c (A 2 ), 

A„ — > A strongly in W Xl => F C (A„) — > T C (A) strongly in W X2 . 

Fc maps observables of system 1 ( r W- Kl ) to observables of system 2 { C W% 3 ) and induces 
therefore also a map from the linear functionals on system 2 to the linear functionals on 
system 1 : if (j)2 is a linear functional on system 2, then the induced linear functional on 
system 1 is r* c (fc)(A) = <h(r c (A)). 

The duality between F c and F^ can be interpreted as the duality of the Heisenberg pic- 
ture and Schrodinger picture in quantum mechanics. Fc reflects the Heisenberg picture of 
quantum mechanics: it describes how observables change. It is then natural to have the 
Schrodinger picture of quantum mechanics implemented by F c and this is possible pro- 
vided that this map describes how states change. At the very least, it should then map states 
to states, in other words, we require that Fc(§2) = <t>2 ° Tc is positive, identity preserving 
whenever (j) 2 is a state. Sufficient for that is that F c is positive and identity-preserving: 

A > T C (A) > 0, 

r c (i) = i. 

Moreover, we would like that normality of states is preserved, so that the resulting state 
can be interpreted as a probability measure. Sufficient for that turns out to be that Fc maps 
vacuum states to vacuum states: 

r c(<lW := *iv 2 or c = <iw 

The latter is not an unreasonable requirement, since the time-evolutions we consider (either 
that of energy-conservative or energy-dissipative systems) and the embeddings (in larger 
quantum systems) or projections (onto quantum subsystems) do not increase the energy of 
the system. From the interpretation of a vacuum state as the state of lowest energy level, we 
do not expect that F c yields a state of higher energy than the vacuum state when starting 
from a vacuum state. 

There is one complication with the above, however: positivity alone is not sufficient to 
make the functorial machine Fc work well. The oil needed is called complete positivity. 

Definition 2.2 

A linear map T between von Neumann algebras J4 and % is called n-positive (bSN) if the 
map T ® I maps positive elements of A ® M n to positive elements of % ®M n , where M n is 
the algebra ofnx n-matrices with complex coefficients. 
T is called completely positive if it is n-positive for all n e N. 

The rationale for this is the following. The map T* describes some time-evolution of (or, 
more generally, a physical operation on) a state of the quantum system described by $ 
yielding a quantum system A in some state. If a system is not affected by this operation 
then naturally the map T* is extended to a map sending states of ® ® HI to states of%.®%.. 
In particular, the corresponding (dual) map T ® 1^ : J? <g> 31 — > S ® must be positive. 
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It turns out that it is sufficient to consider finite-dimensional ^ only, i.e. von Neumann 
algebras over finite Hilbert spaces (C), which are the matrix algebras M n . 
This leads to the following definition: 

Definition 2.3 

Let Si and <B be von Neumann algebras with normal states 0^ and respectively. 
A physical quantum operation on a quantum system (also referred to as a quantum opera- 
tion) described by <B in state (j) 3 that yields a state (j)^ of some quantum system described 
by SL is a linear map T : (SL , 0^ ) — > (« , <|) s ) that is 

i. strongly continuous, 

ii. completely positive, 
Hi. identity preserving, 
iv. state preserving. 

We can now formulate the quantum version of a dynamical system. 
Definition 2.4 

A quantum dynamical system consists of a quantum system described by a von Neu- 
mann algebra and normal state equipped with a family of quantum operations 



Although the input and the output of the quantization functor T seem very different at first 
sight, they bear resemblance on an abstract level: both consist of "objects" (Hilbert spaces 
on the one hand and von Neumann algebras equipped with a normal state on the other) 
and there are "morphisms" between those objects (contractions from one Hilbert space to 
the other versus completely positive, strongly continuous, state- and identity-preserving 
maps between the von Neumann algebras). Such a collection of objects and morphisms 
(together with some properties, which we do not spell out explicitly [10] is called a cate- 
gory by mathematicians. Thus, F maps the category whose objects are Hilbert spaces and 
whose morphisms are contractions to the category whose objects are von Neumann alge- 
bras equipped with normal states and whose morphisms are completely positive, strongly 
continuous, state- and identity-preserving maps. 

As we will see below, T preserves compositions of morphisms and maps special types of 
morphisms in the one category to special types of morphisms in the other, which mathe- 
matically justifies naming it a functor [10] 

First, we will define the action of T in general terms i.e. on Von Neumann algebras gener- 
ated by arbitrary cyclic vacuum representations. 

Theorem 2.5 

Let %i, i= 1,2, be Hilbert spaces and let W; : Kj — > S (#,-), / ^ Wt(f) be cyclic vacuum 
representations of the CCR over Kj (with cyclic vacuum vector £2^ in the representing 
Hilbert spaces #i respectively). Let t Wo(%i) be the linear span of the Weyl operators over 



Xi- 

Assume that C : %_\ — > is a contraction and define To(C) : Tfo(^Ci) — » ^0(^2) by 



Then To(C) has a unique strongly continuous extension to an operationT{C) from r(9Ci) = 



T,\ (*,<!>*) ->(*,<!>*)>' el 



2.3 Mathematical formulation of r 




(2.1) 



More precisely stated, 
i. formula (2.1) determines by linearity unambiguously a map Wo(Xi) — * ^0(^2) and 
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ii. this map is strongly continuous, hence admits a strongly continuous extension F(C) 

from v9i{W (-Ki)} = c M(%i) to v^{W Q (n(2)} = < W{% 2 ), 
Hi. The required properties for F(C) to be an operation (see 2.3) are easily checked on 
WoilQi) and then pass on to <W (ffo) by strongly continuous extension. 

The proof of theorem 2.5 takes several steps and is deferred to the appendix B. 
The following theorem underlines the functorial character of F. One aspect of that needs 
some elucidation: the property of F to map special morphisms in the category of Hilbert 
spaces with contractions to special morphisms in the category of Von Neumann algebras 
with operations. The special morphisms in the category of Hilbert spaces with contractions 
that we consider are isometries, projections, and unitary maps. An isometry J is special 
in the sense that it allows to identify one Hilbert space (#~i) with a subspace of another 
(H) in such a way that all the Hilbert space structure of 9i\ (the linear structure and the 
inner products) is preserved. In other words: an isometry embeds a Hilbert space into an 
enveloping Hilbert space. The analog of this in the category of Von Neumann algebras 
with operations is an injective map that preserves all Von Neumann algebra structure (the 
algebraic structure and the adjoints), a so-called injective *-homomorphism. 
A projection P\ is special in the sense that it decomposes each vector in a Hilbert space 
in two components: one lying in a Hilbert subspace M\ and the other in its orthogonal 
complement #j . In other words, it allows us to pass from a Hilbert space onto a Hilbert 
subspace. If / = P\ : tt\ — > 9{ is the corresponding isometry embedding H\ into H , then P 
is characterized by Pi o / = 1^ . The analog of a projection in the category of Von Neumann 
algebras with operations is called a conditional expectation and is defined with respect to 
an injective *-homomorphism j as the unique morphism P such that Po j = I. 
Unitary maps are special because they allow to identify Hilbert spaces with each other. 
In particular, if unitary maps of a Hilbert space to itself are considered, they can describe 
reversible time-evolutions. The analog of unitary maps in the category of Von Neumann 
algebras and operations are called ^-isomorphism which are bijective *-homomorphisms. 
♦-Isomorphisms that map a Von Neumann algebra to itself are called ^-automorphisms. 
The latter can be used to describe time evolutions in a Von Neumann algebra. 
Summarizing: 
Definition 2.6 

Let {Sl,§n), (S , <j)» ) be Von Neumann algebras equipped with normal states. 

An operation j : (Jl,$ M ) — > that is linear such that j(A\A2) = j(A\)j(A2) and 

j(A)* = j(A*) is called a ^-homomorphism (here * refers to taking the adjoint). 

If j : — ► (®,<|>s) is an injective ^-homomorphism, then the conditional expectation 

with respect to j is the unique operation P : (® , ) — > (j? , (j)^ ) such thatPo j = l A (if it 

exists). 

An operation T : — ► (®,<|>s) is called an * -isomorphism if T is an bijective *- 
homomorphism and an * -automorphism if in addition (S = 

Thus, an injective *-homomorphism allows us to embed a quantum dynamical system into 
a larger one; a conditional expectations allows us to project from a larger quantum dy- 
namical system to one of its subsystems, and a *-automorphism can be used to describe a 
transformation (e.g. by time evolution) of a quantum dynamical system. 

Theorem 2.7 

Assume the setting of theorem 2.5 above. 

Then F preserves composition of morphisms: ifC\ : %\ — > %2 and C2 ■ %2 — > %2, are con- 
tractions between Hilbert spaces 9i\ , #2, #3, then r(CiC2) = T(Ci)T(C2). 
Further, F maps special morphisms in the category of Hilbert spaces with contractions to 
special morphisms in the category of Von Neumann algebras with operations: 

i. ifj : H\ —> 9li is an isometry, then F(J) : F(^H\ ) — > F(y{2) is an injective * -homomorphism. 

ii. ifj* :9i2—>9i\ is a projection, then F(J*) : r(#2) — » r(#i) is a conditional expecta- 
tion. 
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Hi. ifU : #1 — ► !H\ is a unitary map, then T(U) : r(#i) — > r(#i) is a * -automorphism. 



2.4 The usual approach to second quantization 

We conclude this part by relating the functorial approach to second quantization described 
above with the usual approach taken in physics literature. 

First of all, second quantization is used in physics when quantizing fields. Recall that a 
field is a map (observable) whose domain is a space-time continuum § x T and whose 
range V could be any set. Typically, § = M (one-dimensional Euclidean space) or § = R 3 
(three-dimensional Euclidean space) and T is a possibly infinite time interval (T = [a,°°), 
T = (— °°,a], T = [a,b], or T = K). Further, if V C R, then we speak of a scalar field and 
when V C R" with n > 1, then we speak of a vector field). An example of a scalar field is 
the displacement T) (z, t) from the equilibrium of a string, which is a function of the position 
along the string z and time t . As fields are systems with an infinite number of degrees of 
freedom, their quantization is more involved. 

The approach usually taken in physics literature can generally be outlined as follows [9, 
p. 107]: 

1. Let <j) = <j)(x,/) be the field whose dynamics we are interested in. 

2. Due to the continuous character of the space the classical Lagrangian L is a integral 
over a Lagrangian density L 

where is a vector with 4 components, ^(f,x), and V(j)(f,x) if we have a vector 

dt 

field on R 3 . 

3. Derive the conjugate field (field impulse) 
8L(f) dL 



n(t,x) 



«[^] 



where 1^ denotes a functional derivative. 



4. Invert the above to obtain -^(t ,x) in terms of n(t,x) and (j)(f ,x). 

5. Construct the Hamiltonian H as a function of the field (j)(f ,x) and field impulse n(t,x) 



•H ^f,x,(|)(f,x),7t(f,x)^ 



Note that in the above -^(f,x) is a function of §(t,x) and n(t,x) (and of t and x). 
Further, in the Lagrangian density is to be read as the vector with as first component 

di|> 



— (t,X, (j)(f,x), 7t(f,x)) 



and as second till fourth component V(j)(f ,x). Finally, H is a Hamiltonian density. 

6. The Hamiltonian operator is obtained by viewing (j)(f ,x) and (j)(f ,x) as a continuum of 
operators satisfying the equal time commutators 

[<j)(f,x), 7i(f,x')] -/8 3 (x-x') , 



while [0(f ,x),(j)(f',x / )] = [7i(f,x),7i(f',x / )] = 0. 
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Variants of the above sketched second quantization procedure involve going backwards- 
forwards between continuous and discrete versions of the system at steps 2 and 5 respec- 
tively. 

Often, obstacles arise when applying this quantization procedure, which physics tradition- 
ally bypasses "along the road": 

1. Naive substitution of operators for and % may lead to a non self-adjoint Hamilto- 
nian. Problems of these kind have to be solved ad-hoc by hermitization procedures. 

2. Nothing has been said about the Hilbert space on which the fields and n act, but in 
the particular case of the free scalar field one can fortunately do without (see below 
andB.9). 

3. Multiplication of operator fields at the same point may introduce new problems, such 
as infinite ground state energy levels. Problems of the latter kind have to be tackled 
ad-hoc by renormalization procedures (putting the system in a finite box with peri- 
odic boundary conditions, and see if the ground state energy level is proportional to 
the size of the box. If so, then one simply subtracts the infinity caused by the size of 
the system). 



2.5 Quantization of the free scalar field: the "usual" way 



i 

1 

y / 4rao(Fj' 



As a example to make the comparison between the "usual" approach to second quantization 
and the functorial approach easier, we will discuss the quantization of a (quadratic) real 
scalar field (j) which yields a self-adjoint (quantum) field. Our approach follows that of [9, 
p. 107-114] to a great extent, but we will only mention the main facts, leaving the details 
to the appendix (see B.9). 

The Lagrangian density of the quadratic free scalar field (j)(x,f ) reads 

i = I(30) 2 -ImV, 
leading to the Hamilton density 

H = j dx { I [tt(x) 2 + (g (*)) + m>(*) 2 } , 
which reads after Fourier transformation: 

H = h. fdk^(k)%(k) + [m 2 + k 2 ]^(k)avk ) 

After defining 

a{k) : = 

a\k) 
we obtain 

H=\P^dk CO(fe) a i {k)a{k)+a{k)a ! {k) . 

The vacuum state |0) has infinite energy and has to be renormalized by applying a Wick 
ordering of the Hamiltonian: 

H = f~ 00 dk <a{k)a\k)a{k). 

Excited states where only the oscillator corresponding with "mode" k contains an energy 
quantum can formally be created as a(k) |0), but have infinite energy. Therefore, it makes 
more sense to build "wave packets" 



(d(k)$(k) + in(k) 

C0(^)$ t (fe)-OT t (/t) 



+ 



|0):= f dkf{k)a\k)\Q) 

J oo 
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which are normalizable provided that 

- > ^Ola/^lo) = Jdk Jdk 1 fW)f(k) (Ojj(Jf )a*(k]0) 

at(*)a(f)+8(t-* / ) 

/■co 
d*l/(*)l 2 , 
-co 

i.e. the "density kernel" / has to be square integrable: / e L 2 (R). Coherent states are 
obtained as: 

Coh(/) :=K(f)ex P (fdkf(k)aHk))\0) , 

- 1 II fll 2 

where K(f ) = e II " "l 2 is a normalization factor. 
Thus, the Weyl operators are proportional to 

W(f)=exp(a}) (f € L 2 (M)). 

Thus, the usual quantization of the free scalar field corresponds to functorial quantization 
(via T) of the phase space L 2 (M) of an infinite number of oscillators, indexed by k e M. 



Part II 

Damped Harmonic Oscillator 
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Chapter 3 



The problem of irreversibility 



Some fundamental laws in physics, the Schrodinger equation, are invariant under time 
reversal. This has the following implication for a physical system that obeys these laws: we 
cannot distinguish (by no physical quantity whatsoever) whether time is running 'forward' 
or 'backward' . Another way of saying this is that in principle the course of events can be 
reversed: none of the changes that the system has gone through is permanent. Owing to the 
reversibility of its time evolution, such a system called a reversible dynamical system. 
An example of a reversible dynamical system is a classical harmonic oscillator. An har- 
monic oscillator is (fully) characterized by its frequency co, amplitude A (and mass m) and 
it has the following time evolution 

q(t) = A sin G)f , 
p(t) = d - 3 ^ L = mAcocosco/, 

If we replace in the above equations t by — t (which can be realized in practice by giving 
at time the oscillator the same speed initial speed but now in the converse direction), we 
again obtain a harmonic oscillator with the same frequency and mass. From a physical 
point of view, we obtain the same harmonic oscillator: the fact that the direction of the 
oscillation has flipped from 'positive' to 'negative' has no physical meaning because if we 
have two frames of references that can be obtained from each other by a rotation, then there 
is no physical way of distinguishing the one from the other. Note that the average energy of 
the harmonic oscillator over one period which is proportional to the amplitude, is conserved 
quantity in this case, and this conservation of energy is a feature of reversible dynamical 
systems. We may think of this as follows: the flow of energy stays within the system and 
cannot leave it. Therefore, the reversible systems are also called closed systems. 

There is, however, a major problem with time reversibility. In nature most courses of 
events are not reversible: who has ever seen a damped harmonic oscillator that started from 
rest position and began to oscillate with increasing amplitude (without applying external 
forces)? 

Since the reversible theories are too good to discard, physicists (from Boltzmann on wards) 
have taken up the challenge to explain irreversible behavior from the reversible theories. 
One point of view (the open system interpretation) is the following. First observe that an 
irreversible dynamical system looses energy to its surroundings and that we may think of 
this as follows: the flow of energy is not confined to the irreversible system, but leaves it 
and enters the system's environment. For this reason a irreversible system is called an open 
system. 

An example of an irreversible dynamical system is the classical damped harmonic oscilla- 
tor: 

q{t)=e-VAsmm, ) , , 
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The average energy of the oscillator over one period (roughly proportional to its amplitude 
e^A) decreases and the lost energy is found back as heat caused by friction and deposited 
to the oscillator's surroundings. 

If we take the environment large enough, then the flow of energy of the total system con- 
sisting of irreversible system plus environment, stays within the total system. Then the total 
system is a closed system again and it is expected to have reversible time evolution. 
In case of the damped harmonic oscillator, we could imagine that it is a point on a string 
on which a wave packet is traveling. As the (center of the) wave packet moves away from 
the point, the oscillation of the point is fading away with it. 

Summarizing, the open system interpretation states that a irreversible (open) system is a 
subsystem of a reversible (closed) system and that the time evolution and the observables 
of the irreversible system are 'projections' of the time evolution and the observables re- 
spectively of the reversible system. 

Let us now formulate the open system interpretation in a mathematical language for both 
classical and quantum dynamical systems. 



3.1 Dilation theory of dynamical systems 

Dilation derives from the Latin "dilatare", which means "to make broader". Dilation theory 
aims to give a mathematical description of embedding an irreversible dynamical system 

(s ATs,t) s ,te[o,°o)j into a larger reversible dynamical system (s ^(T^t),,^^^ . 
We can depict this symbolically as follows: 

T 'h<2 



s — s 

Here S and S describe the observables of the subsystem and the supersystem respectively. 
T SJ : S — > S is the irreversible time-evolution of S (hence s,t £ [0, °°)) and T s j : S — + S is the 
reversible time-evolution of S (hence fi,?2 £ R). Note that a description in terms of states 
can be obtained by taking the dual spaces S* and S* and translating the time-evolutions 
correspondingly. 

The embedding of irreversible system into the reversible one at time t\ is described by a 
one-to-one map j : S — > S that identifies S ~ j tl (s) as a subsystem of S . A corresponding 
map P t2 : S — > S has to exists that projects S to S such that P t o j, =l s for all tel. 
Furthermore, it must be possible to view the evolution T t in the irreversible system as a 
projection of an evolution in the reversible system: T tu , 2 = P t2 o T tht2 o j t[ . 
In the following we will restrict ourselves to the situation that the time-evolutions are 
translation-invariant: T s , s+t = To. t =: T, for all s,t £ [0,°°), T S:S+t = 7b jf =: T, for all s,teK, 
and moreover j t = j, P t =P i.e. 




3.1.1 Dilations of classical dynamical systems 

Let (9C,C t ) be a classical dynamical system (see 2.1), so a member of the category of 
Hilbert spaces with contractions. As isometries, unitary maps, and orthogonal projections 
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are the embeddings, isomorphisms, and projections respectively in this category (see the 
motivation before 2.6), a dilation of a classical dynamical system (HC,Q) is a quadruple 

(j , %. , Q , J* ^ ) where J is an isometry %, — > 2£ , /* : 3C — > %. an orthogonal projection, and 
Q, f e R a one-parameter unitary group. Summarizing: 



X. L - » 

7* 

3C - ' ?C 

Given an irreversible dynamical system, the question arises whether there exists dilations 
and if so, if there are more than one. For special classes of irreversible systems, general 
constructions for dilating them have been found. In general, if there exists one dilation, 
there are more which are not necessarily unitarily equivalent. To force uniqueness, so that 
the properties of the dilating system (S ,T t ) are determined by the irreversible subsystem 
(S,T t ), the concept of a minimal dilation is useful. 
A dilation is called minimal if J{lCj is dense in 

As example, we present a dilation of a "symmetrically" damped harmonic oscillator which 
we will use later. 

Theorem 3.1 (Sz. Nagy, Foias 1953; special case) 

Let ac = C be the phase space of one oscillator and consider the damped time evolution 

Q : C -> C, C,(z) = e (-?+ to >z, (r > 0). 

Unique up to unitary equivalence there exists a Hilbert space 3C with a unit vector v and 
a one-parameter group of unitary transformations T t =: U t on 9C such that the span of the 
vectors U t v,t € R, is dense in % and 

(v\U t v) = e (-y+ iw )>, (f>0). 

Therefore, we have the following minimal dilation of (C,C f ): 

r t=e (-l+i«>)x 
C — ► C 

~ u, ~ 

X — X 

The existence follows e.g. from taking Ki := L 2 (R 7 2ydx), U t := R t the shift to the right, 
and 



v(x) :-- 



ifx>0, 
e (i\-ua)x if x < o, 



while the uniqueness follows from Kolmogorovs dilation theorem (see the last example of 
1.14). The minimality follows from the bijectivity of the Fourier transform. For details see 
A.22. 

In practice several (unitarily equivalent) minimal unitary dilations of (Q) t can be useful. If 
3C = L 2 (R) and U t is the shift, then we speak of translation dilations of (Q) t . They differ 
only in the shape of v e L 2 (R) which must satisfy (in terms of their Fourier transforms) 

Particular solutions are v±(X) = l/(k — <0±z'r|). Here v+, which occurred above, leads to 
the incoming translation dilation and v_ to the outgoing translation dilation: 

( e -(x\+im)x ifx > 0j 

v_m := < 

10 ifx<0, 
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The former is more useful for the study of incoming fields and particles, the latter for 
outgoing ones. The unitary equivalence of these two unitary dilations, asserted by theorem 
3. 1, is implemented by the scattering operator S which in terms of the Fourier transform F 
can be written as 

7 -l „n \ X— co+iT) 



FM s F-\ *(*):= j^g. 



3.1.2 Dilations of quantum dynamical systems 

The embeddings, isomorphisms, and projections in the category of quantum dynamical sys- 
tems (2.4) are injective *-homomorphisms, ^-isomorphisms, and conditional expectations 
respectively, as has been elucidated in the the motivation before 2.6. Therefore a dila- 
tion of a quantum dynamical system ^(j? ,(()), 7^ is a quadruple ^ j, ,(()), T t ,P^j such that 

T t = Po T, o j, where j : — > (j? ,<|>) is an injective *-homomorphism, P the associated 
conditional expectation, and T t : — > t e R, is a family of *-automorphisms. 

Symbollically, we can depict this as follows 

(*,<>) — ^ (*,<>) 

j P 

(5,$) —5— 

Although not mentioned explicitly, it is understood implicitly that all morphisms are quan- 
tum operations, as we work in the category of quantum dynamical systems with quantum 
operations. In particular, j, P, and T t are strongly continuous. 

We will meet a quantum dynamical dilation in just a moment, when we quantize the 
damped harmonic oscillator. 



Chapter 4 

Quantum Mechanical Damped 
Harmonic Oscillator 



(re [OH), 



(4.1) 



A technical problem when quantizing the classical damped harmonic oscillator is that the 
roles of the position q and momentum p are not symmetric, because their time-evolutions 
differ more than a simple phase shift. We will quantize a symmetrical spiraling motion in 
phase space instead: 

q(t) =e _T|f (g(O)cosCOf +p(0)sincof), 

p(t) = e~ n ' (-q(0)sma>t +/?(0)coscof) 

i.e. we take m = 1 and discard the term proportional to r\e^ t in p(t) = (weak damp- 
ing). Although this may seem a crude approximation at first sight, it retains the essential 
features of a damped evolution and it can be viewed as a "Van Hove" limit of the damped 
harmonic oscillator 

In view of the above and the Heisenberg equations of motion, we formally define a quantum 
damped harmonic oscillator as a canonical pair (Q,P) obeying the time-evolution 



T,{Q)=e ^ (Q cos (Ot + P sin (at), 
T,(P) = e^' (-gsincof +Pcoscof) 



(te [0,oo)). 



(4.2) 



4.1 Quantization of the damped harmonic oscillator 

Identifying the phase space {(q,p) : q,p € R} with C, we can describe the damped spiraling 
evolution (4.1) as 

C,:C-C, z-eHl+^'z, (fe[0,oo)), 
and quantize it using the second quantization functor F to obtain 

r r :-r(c f ):r(C)^r(C) (fe[o,-)). 

Then 

T t (W(z)) - e^'-^Wie^+^'z), (4.3) 
and T, is defined on T(C), the whole von Neumann algebra generated by the Weyl operators 
associated to Q and P. 
If we define forx,y e M, 

xP-yQ 

as the Stone-generator of the strongly continuous one-parameter unitary group 

(W(X(x + iy))) Xm i.e. 

W(X(x + iy)) = e - iX ( xP -y® (X e R) , 
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then we have a natural way to let T, also act on the unbounded self adjoint operators 

xP-yQ. 
Indeed, for x,y <G R let 

T t (xP - yQ) = T t (ijL [w (X(x + iy))] ^ = J 



= i 



dX 
.d_ 
dX 



T t (W(X(x + iy))) 



using (4.3) 



(4.4) 



X 2 (e- 2 * - 1) -iX [Re (e im (x+iy))P-lm (e iat (x+iy) ) q] 



X=0 



= x[— sincof g + coscof P] +>'[coscor Q + sinatt P] 

Formally setting 

T t (xP -yQ)=x T, (P) -yT,(Q) (4.5) 

and equating the coefficients of x and y in (4.4) and (4.5) , we recover the equations (4.2). 
In this sense, T t is a formulation of the damped harmonic oscillator evolution in terms of 
Weyl operators that is equal to that given by (4.2) in terms of Q and P. 

Remarks 4.1 

Note that in a vacuum representation W: C — > U (tt ) of the CCR over C (thus the Gaussian, 
Schrodinger, and Heisenberg representations), the von Neumann algebra T(C) generated 
by all Weyl operators coincides with all of S (tt). 



4.2 Quantum dilation of the damped harmonic oscillator 

We now focus on dilating the damped oscillator evolution to an evolution determined by 
a Schrodinger equation. Taking the minimal dilation of (Q)t>o as in theorem 3.1 and 
applying the second quantization functor F to it, we obtain the following dilation of T t . 

r(C) — r(C) 

j p 

r(L 2 (K)) — r(L 2 (M)) 

Here T t = T(C,), j = T(J) with J : z ^ zv and P = Y(J*) : f ^ (v\f). Since f, := T(R t ) is a 
point wise continuous one-parameter group of *-automorphisms of the algebra F(L 2 (R) ) = 
S (L 2 (1R)), it is implemented by a unitary strongly continuous one-parameter group (U t — 

e itH )teR as 

f t (A) = U t AU*. 

So the evolution (T t ) t obeys a Schrodinger equation and the irreversible evolution 

(T t =Pof t oj) t 
is a 'shadow' of the reversible evolution (T t ) t . 

It is now a good moment to lend some physical flesh and blood to the mathematical skeleton 
depicted in the diagram above. 



4.3 Interpretation as quantum dynamical system: 
open system with (minimal) Bose dilation 

The semi-group (T t ) t expresses the time evolution of the damped harmonic oscillator, al- 
though not in terms of the unbounded observables Q and P, but in terms of their associated 
Weyl operators. T(C), which is the von Neumann algebra generated by the Weyl operators 
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associated to Q and P, is the algebra of all (bounded) observables of the damped harmonic 
oscillator. Implicitly, the state (j) c on T(C) is assumed, which describes that the damped 
harmonic oscillator is in its vacuum state. 

We may wonder how we could physically interpret this dilation of the quantum harmonic 
oscillator. From the discussion in section 2.4, it is clear that we have emerged (or coupled) 
the damped harmonic oscillator in a field of oscillators labeled by / <G L 2 (M.). Due to the 
field's commutation relations, this is a Bose field. Keeping in mind that the damped har- 
monic oscillator should model an atom loosing its energy to its surroundings by emission 
of light, the interpretation of r(L 2 (M)) as a Bose field of light quanta is readily made. The 
time evolution of r(L 2 (M)) is given by ^-automorphisms, and therefore reversible. These 
♦-automorphisms are a generalized way of saying that the time-evolution is given by a 
Schrodinger equation. 

Note, however, that we could have taken the Bose field too 'large' in the sense that some 
parts of the field do not interact with the oscillator. Removing the redundant parts of the 
Bose field, we are left with a quantum mechanical system called the minimal Bose dilation 
of the damped harmonic oscillator. 

The minimal Bose dilation is unique in the sense that any quantum mechanical system that 
evolves according to a Schrodinger equation and contains the damped harmonic oscillator 
in the sense that 

i. every oscillator observable can be recovered as a 'projection' of an observable of the 
total system onto the oscillator subsystem, 

ii. the time evolution of the damped oscillator is a 'projection' of the overall time evolution 
onto the oscillator subsystem and 

iii. the interaction of the oscillator with the surrounding system is via a Bose field, 
also contains an isomorphic copy of the minimal Bose dilation. 

Now r(L 2 (R)) and (T, = T(R t )) t are the observable algebra and time evolution respec- 
tively of the minimal Bose dilation of the damped harmonic oscillator. 
On r(L 2 (R)), the state L 2( R ) is assumed. 

As mentioned in ii., every observable of the damped harmonic oscillator is a projection 
of an observable of the total system (=oscillator and minimal Bose surroundings). Now 
j associates to an oscillator observable the corresponding observable of the minimal Bose 
dilation, while P projects observables of the minimal Bose dilation onto oscillator observ- 
ables. 

4.4 Interpretation as quantum stochastic process 

No one who has been exposed to quantum mechanics seriously has failed to notice that 
the physical content of quantum mechanics is of stochastic nature. We cannot predict the 
outcome of a physical quantity in a single experiment, but we can predict the probability 
distribution of an observable. The latter is the relative frequency with which each of the 
possible outcomes of measuring that observable occurs during an (infinitely) repeated se- 
ries of measurements on identical systems. Such a collection of identical systems is called 
an ensemble, and arises e.g. because of an identical preparation procedure. The merit of 
quantum theory is that it postulates rules how to obtain and calculate with those probabili- 
ties, and that these rules turn out to aptly predict how small scale physical systems behave. 
Formulating these rules more comprehensively on an abstract level leads to quantum prob- 
ability theory in which the concepts of quantum probability spaces and quantum stochastic 
processes enter stage. 

Quantum probability space Recall that a quantum dynamical system (a , <j), T t ) consists 
of two parts: a Von Neumann algebra A with normal state (j), and a family of quantum 
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operations (T t ) t . It is the first part that can be interpreted as a quantum probability space as 
soon as we have identified what the events and their probabilities are. 

Definition 4.2 

A stochastic experiment is said to be modeled by a quantum probability space if 
the following is the case. The experiment can be repeated arbitrarily often and sufficiently 
many of the orthogonal projections E in A have an interpretation as an event i.e. a state- 
ment about the occurrence of certain outcomes of the experiment that can be tested by 
observation. The relative frequency with which the event encoded by E occurs is given by 

m- 

The projection is the impossible event and I is the sure event. Many events relate to the 
outcomes of observables of the system. If A is an observable andP^O) its projection- valued 
measure, then Pa ( (— °°, s] ) is the question whether a measurement of A yielded an outcome 
that was less than or equal to s. 

In the picturesque language of Mackey [ 1 2] events may be considered as questions that one 
can ask the system. 

In contrast to classical stochastic experiments, as envisaged by Kolmogorov, we allow ad- 
justment of the observation equipment between the trials, that is, we are allowed to ask new 
questions. When asking different questions, a quantum peculiarity arises. If two questions 
E,F are compatible i.e. [E,F] — 0, then those questions can be asked together in the same 
trial, and EF denotes the event that both E and F occur, and EV F := E + F — EF is the 
event that either E or F occurs. 

Incompatible questions can never be asked together. They can be asked in different trials 
by readjusting the apparatus in between. Inside one single trial it is sometimes possible 
to ask incompatible questions one after the other, but then the fact that one question was 
posed, may influence the answer to subsequent questions. 

If the questions 2?i,2?2> • • • are asked in each trial and in this order, the asymptotic frac- 
tion of the trials in which they are all answered 'yes' is 

d?(EiE 2 ■ ■■E k -iE k E k -\ ■ ■■E 2 Ei). 

The generality of the concept of a quantum probability space is illuminated by the fact that 
a (classical) probability space is included as special case (which we shall not prove): 
Theorem 4.3 

Let (Q,E,P) be a probability space and let Mf be the multiplication by f e L°°(£2,E,P). 
Then the algebra 

J4 := {Mf : f e L°°(i2,E,P)} 

is a (commutative) von Neumann algebra of operators on L 2 (i2,E,P), and (j) : Mf i— » J fdP 
is a faithful normal state on A. Conversely, every commutative von Neumann algebra A 
with a faithful normal state is of the above form for some classical probability space. 

Taking / = Is for S € E, it is clear that (j)(/) = P(5) as should be. The question arises how 
random variables are encoded in terms of (commutative) Von Neumann algebras. 

Quantum random variable Classically, a random variable with values in (R,Borel(R)) 
is a measurable function X : (£1, E) — > (K, Borel (R)). We can alternatively characterize X 
by the o-algebra embedding 

X' 1 :R->I.:S^X- l (S):={®eQ.:X(a)£S}. 

The probability distribution Px of X is then the induced measure on R: 

P x - Pol- 1 . 

In fact, giving X is equivalent to giving the injective *-homomorphism 

jx :L~(M, Borel (R),P X ) ^L°°(i2,E,P) :/h/oX, 
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because = §oj x where and are be the (standard) choices for the states on L°°(K, Borel(R) 
and L°°(i2,E,P) i.e. the expectation with respect to P x and P. 

Definition 4.4 

A quantum random variable on a probability space (si , 0) is an injective *-homomorphism 
j from some other von Neumann algebra $ to Si which maps L s to L, . The state \|/ := o j 
on % is called the probability distribution of j. 

We say that j takes values in if S = L"(Q',L'). The random variable j is then 

completely determined by the projection- valued measure 

In particular, we see that quantum probability generalizes the events and random variables 
as known in classical probability theory. 

Quantum stochastic processes 
Definition 4.5 

A quantum stochastic process is a family (j t )teT °f random variables (J?,0) — ► (-#,0) 

indexed by time. Here, T is a linearly ordered set such as Z, K, N or R + . If T — M. or M + 

we require that for all A E A the curve t i— » j t (A) is strongly continuous. 

If T is a group, say 7L or K, then the process is called stationary provided j, = jo o T t for 

some representation 1 1— » T t ofT into the automorphisms of (A,0). 

Remarks 4.6 

If a classical process is stationary, X t and X s differ by an automorphism i t of the underly- 
ing probability space: X t = Xo o x t . These remarks generalize to quantum probability, as 
quantum stochastic processes only differ by an automorphism of Von Neumann algebras 
(=quantum probability spaces). 

The idea behind dilation (open system interpretation) is that we are observing a subsystem 
with observable algebra Si of a larger environment with algebra A that we cannot see. 
According to the Heisenberg picture, the smaller algebra is moving inside the larger one. A 
question E e A can be asked at time t only if E G jt{&)- If to < h < ■ ■ ■ < t n is a sequence 
of times, and E\,...,E n a sequence of events in Si, then 

?(i'r, {E\)ji 2 (E 2 ) ■ ■ ■ h n -x (En- 1 )jt„ {Entity ■ ■ ■ jt 2 (E 2 )j H (E i )) 
is the probability that E\ occurs at time t\, E 2 at time t 2 , ■ ■ ., and E n at time t n . 
Can we calculate this probability if we only observe the subsystem (i.e. if we cannot eval- 
uate and j t )7 

Two-time-probabilities Suppose that for all s e T there exists a conditional expectation 
P s with respect to j s , that is, not only do there exist embeddings j t of the subsystem into 
the environment, but also the corresponding "projections" that associate observables of the 
total system with observables of the subsystem. 

Then the probability for F to occur at time s and E at time t > s can be written as 

$(j s (F)j t (E)j s (F))=$(FP s (j t (E))F)=$(FT s , t (E)F), 

where T sf — P s o j, is an endomorphism of (-#,0), the transition operator from time s to 
time t . In other words, the two-time probability can be "calculated" from information that 
is available within the subsystem. 

Multi-time-probabilities The above reduction to the subsystem succeeds for more than two 
time points if there also exist conditional expectations onto the algebras 

*(—,»] := v# {./,(.») :*<*}. 



and moreover the Markov property holds: 

*<* => P{—*]Ut(A))Cj s (x). 



(4.6) 
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Theorem 4.7 

Let (j t : ( A , (j)) — > (j? ,(())) be a stationary Markov process with conditional expectations 
P t . Then the transition operators form a monoid: 

0<s<t<u =*> T s jT t]U = T^ u . 

In particular, if the process is stationary, then T t := Tqj = T s ^ +t satisfies 

T s T, = T s+t {s,t>0). 

In the latter case, (r f ) f >o is known as the dynamical semigroup induced by the stationary 
Markov process. Conversely, the process (j t )teT is called the quantum Markov dilation of 
the dynamical semigroup (T t ) te j. 

The situation is symbolized by the commutative diagram 

(*,<>) — 5-> (*,<>) 
j p 
(*,$) — (5,$) 



Embedded classical processes In case j?oC^, and J4.Q c A are commutative Von Neu- 
mann sub algebras so that the following diagram (with the appropriate restrictions of the 
states and operations) 

T, 



0*0, 4>) 
j 

(So,?) 



(■ao,4>) 
p 

(5 ,$) 



commutes, then there exist classical probability spaces (R, $ and (Q,E,P), such that 
L-(J?,«,/i) — L°°(R,-B, M ) 



L-(£1,E,P) 



L"(Q,L,I 



In this case, we can often identify a stationary classical Markov process or chain (X,),, as 
follows: 



T,:f^(x^f f(y)P,(x,dy)) 



L°°(Q,E,I 



r,:F^Fox, 



/>:y^E(r|X ) 

L°°(Q,E,P). 



where 

i. 7} : corresponds to an automorphism i t of the underlying measure space (Q,E,P): 
7](T) = y o t,; in particular: X, = Xo o x, is stationary random process on (£1, E) with 
values in (R, S); 

ii. the embedding ; maps / e L°°(#, 2 to /(Xq) e L°°(i2,E,P); 

iii. the conditional expectation P is closely related to the classical conditional expectation 

y^E(y|x ) 

i.e. the conditional expectation given the process at time 0; more precisely: 

P(Y)=geL-(R,V, t i) E(Y\Xo)=g(Xo); 
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iv. the maps (T t ) t form a semi-group of transition probabilities; if P t (x, dy) is the probabil- 
ity that the process (X,) t is in the area dy at time t given that it started fromx at time 0, 
then 

(T t f)(x) = J f(y)P,(x,dy) (x€R,fe L~(/?, *,,,)); 

v. the equality in the upper right corner of the diagram: ¥,(f(X t )\Xo) = (T t f)(Xo) = 
I f(y)Pt(Xo,dy) is also intuitively obvious, as the expected value of f(X t ) given that 
the process started at Xq is the average of all possible f(y) weighted by the probability 
that X, would arrive at y at time t , when it started from Xq. 

4.5 Stochastic behavior of the oscillator 

Let us return to the quantum dilation of the damped harmonic oscillator (section 4.2), where 
T, = r(Q), f t := T(U,), j = F(J) with J : z h-> zv and P = F{J*). Due to the functorial 
character of F (B.8), this is automatically a quantum Markov dilation. Furthermore, the 
functorial character of F also implies that we can split T t as 

T t = F(e^+ ia >) = F{e iw )F{e-^). 

The operator F(e~' m ) is the automorphism of the quantum harmonic oscillator studied 
previously. Let us now look at the 'dissipative' part F(e^'). 

Proposition 4.8 

For < c < 1 the operator F(c) leaves the abelian subalgebras generated by I and any of 
the operators xP — yQ invariant. In particular its action on the algebra 

d:={f(Q):feL°°(R)} 

is given by 

2 

r(c)(/(g)) = -7=3^ / f(cQ+x)dx. (4.7) 

Proof 

Obviously, F(c) leaves the linear span of {W(Xz) : A- e K} invariant, and then also its strong 
closure. Putting f(Q) = e' y ® we see that the right hand side of (4.7) equals 

( 1 r e~^ z ^)e ixy dx\e icyQ = e~^ 1 ~ cl ^ 2w ^ 

\y/2n(l-c 2 ) J — ° ) 

which is equal to the left hand side by the definition of F. The theorem follows from the 
strong continuity of F(c) . □ 

We recognize the semigroup T t of transition operators restricted to Q, as the transition 
operators of a diffusion on M with a drift towards the origin proportional to the distance to 
the origin. 

4.5.1 The driving field 

It is interesting to look at the dilation of the semigroup 7] for co = in the incoming dilation. 
Then we obtain an abelian subalgebra of T(L 2 (M, 2r|)) by putting 

Q:=v*C{T,oj(e is Q):t,seR}. 
In the Gaussian representation of the algebra A = r(L 2 (M)) this is just 

Q,:=f t oj(Q)=J^e^dB s , 
where B s is a Wiener process with variance 2r|. 

Indeed, let <£(/) for / e L 2 (R,2r() be defined as the generator of W(Xf) i.e. e m ^ = 
W(Xf). Observe that <t>(f) is a random variable that has a normal distribution with mean 
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and variance ||/|| 2 in the vacuum state (because its characteristic function is A, i— > (W(Xf)) ■ 

1 1 2 II f\\2 

e -i K 11/11 ). Furthermore, the Weyl relations imply that <£>(/) and <t>(g) are compatible if 
(f\g) is real, and independent if (f\g) = 0. In particular, if % = L 2 {R,2r\dx), then 



B, :-- 



*(l [0/] iff>0, 
*(-l[,, ] iff<0, 



is a stochastic process with compatible normally distributed independent increments having 
variance 

$ K ((B t -B s ) 2 )=2t\\t-s\, 
i.e. (B t ) t is a Wiener process. 
Theorem 4.9 

The process Q t := <t>(v t ) satisfies the integral equation 

Q, - Qs = -r| f Q u du +B t - B s (s<t), (4.8) 

J s 

which is the integral version of the stochastic differential equation 
dQ t = -r\Q t dt + dB t . 

So we find an embedded Ornstein-Uhlenbeck process in our Markov dilation. 
Proof 

The following equality between functions in % = L 2 (K,2r|) holds 

U,v- - U s v- = -T) // U u v-du + 1 M {s<t). 
Acting with 4> on both sides of the equation yields (4.8). □ 



4.5.2 Quanta 



We now concentrate on another abelian subalgebra of T(C) = %(?{), namely 

K :={Eo«^:a6r(N)}, 
where E n is the orthogonal projection onto e„ in the Heisenberg representation ttneis = 
£ 2 (N, ^f). So 5\£ is the algebra of all diagonal matrices in Heisenberg's matrix mechanics. 
This time we do not put co = 0. Let (j)„ denote the state A i-> tr(£„A). Write H = 2N + 1. 
Proposition 4.10 

The algebra 5\£ is invariant for T t and 

$„(T t (s N )) = (l-e-^(l-s)) n . 
This is the probability generating function of a pure death process [7] with the generator 
TO • 



L = 2r| 



1 








Proof 

We can write s N as a weak integral over the operators W (z) (lemma 1 .30, appendix A. 1 8) 



s N = - f e-^M 2 W (z)X 2 (dz 

K JC 

where A-2 denotes the two-dimensional Lebesgue measure on 



(4.9) 



This relation can be ex- 
plicitly checked by taking matrix elements with respect to coherent vectors. Application of 
T, to both sides of (4.9) yields for all u, v £ C, 

(n(u)\T t (s N )n(v))=e av ( 1 - e - rv ( 1 - s ». 
Since this expression is not sensitive to the relative phase of u and v, the operator T t (s^) 
lies in 9t . The statement is proved by taking coefficients of (uv) n on both sides. □ 
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4.5.3 Emitted quanta 

Finally, let us have a look at what happens outside the oscillator while it is cascading 
down its energy spectrum. Since we are interested in outgoing quanta at positive times, we 
consider the outgoing translation dilation of (T t ) t >o- The Weyl algebra of T(L 2 (R)) has a 
suggestive interpretation in its Fock representation L 2 (A(R) . Let N t be the number operator 
dj (P Vt ) that counts the excitations of the oscillator at time t . 

As C\[ ~ £°° is invariant under the action of T t = r(e^ T|+!t0 ^), we can obtain a dilation of 
T t : 5\£ — > 5\£ by restriction: 



r(ac) r(ac). 

However, since for different times s and / the functions v s and v r are neither orthogonal 
nor parallel (i.e. multiples of each other), the one dimensional projections P Vs and P V) do 
not commute. A fortiori the associated number operators N t := dj (P Vt ) do not commute 
either, since 

e iXN, e i^N s _ rp [g{kP vt rp ^ e ipPvs ) — j ( e OJ'vte iljPv ') ^ f (e i,jPv * e'^t ) = e^e ^' . 
Nevertheless, we can identify a classical process as follows. For each / £ R, we can con- 
sider three number operators: 

N t :=dj (P Vt ) the number of quanta in the oscillator, 

M t :=dj (M i[f ^ ) the number of quanta that have not yet left the oscillator, 

K t :— df {Mi, f| ) the number of outgoing quanta that have left the oscillator. 
Note that in this choice of representation for the dilation (the outgoing dilation), M t —N t i.e. 
the number of incoming quanta, is not given by a multiplication operator, but the number 
of outgoing quanta is (hence the name outgoing representation). Note furthermore that the 
operators M t and K s do commute for all t , s £ R, hence generate commutative von Neumann 
subalgebras of T(ac). 

For positive times the number of incoming quanta M t — N t has expectation in the states 
of the form 6oP(9ef' (N)) that we consider. So we may expect that replacing N t by M, 
would lead to an embedded classical Markov chain. 
Proposition 4.11 

Fort > 0, wehave the following commuting diagram involving commutative von Neumann 
algebras: 

M — - — ► M 



p=nr) 



where M is the von Neumann algebra generated by T t (s M ) (t > 0, s £ (-1, 1)). Further- 
more, <M = 5\£ . 

Proof 

For all z, u £ C, s £ [0, 1] we have that 
(%{u)\Por{U t )o (s m )\k{z)) = (%{u)\s M '\%{z)) 

= L:=olA nm ^\v^{o)\ 2 ^{do) 

= lT n= ^uzY ■ ± /a„(m) (n« =1 *W^ i V { rj ) dn... dr n 
= exp(2r)Mz/ oo e- 2l l'"/['H (r)dr) 

= exp(uz{[-e-^%-[se-^]T)) 
= exp{uz{l-e- 2 V(l-s))), 
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and the latter equals (k(u)\T,(s n )\k(z)). 

The equality of von Neumann algebras fM and 9\£ follows from the fact that M[ ()00 j is a 
strongly continuous limit of linear combinations of P Vs , s <t □ 

Finally, since K t + M t is equal to the total number of quanta, which is the same constant 
for all times (given a state), we conclude that a quantum is emitted precisely at the mo- 
ment that the oscillator makes a downward jump. Moreover these jumps are made at inde- 
pendent exponentially distributed random times. These phenomena turn out to be natural 
consequences of damped harmonic motion in a non-commutative description. 



Epilogue: other models for 
the damped harmonic oscillator 



It took some time after the conception of quantum theory before it became clear that dis- 
sipation was not only a macroscopic phenomenon, but that there were also microscopic 
systems with dissipative behavior. For example, the quantum mechanical analysis of ra- 
diation in a microwave cavity needs to incorporate dissipation, since cavity losses are not 
negligible (Senitzky). Other examples occur in nuclear physics: nuclear fission, giant res- 
onances, and heavy ion collisions (see p. 7 of Dekker for references). 
Since then the interest in developing quantum mechanical models capable of describing 
dissipation emerged. The damped harmonic oscillator (dho) functions as a standard exam- 
ple for quantum mechanics of open dissipative systems, much like the harmonic oscillator 
is the standard example for quantum mechanics of closed systems evolving via unitary 
time-evolution operator. A steady flow of articles on quantization of the dho began in the 
1940's and has not faded away today. The 1980 review of Dekker mentioned already 232 
articles on the subject (see Harris), but concluded there was still no generally accepted 
way of quantization. Some consensus seems to have grown, owing to the use of the word 
"canonical" in some more recent articles. Based on a non-random, non-systemic literature 
search, we give an impression of some of the major approaches to quantize the dho, with- 
out claiming completeness in any sense, however. (For a completer historical survey up to 
1980, see Dekker). 

One of the first approaches could be summarized as based on quantizing the classical equa- 
tion of motion. The basic idea in this approach is to start with the classical equation of 
motion for a damped harmonic oscillator, find the Hamiltonian or Lagrangian which leads 
to these equation of motion and then quantize by conventional formal methods. 
One exponent of this approach is the so-called (Bateman)-Caldirola-Kanai model (Bate- 
man, Caldirola, Kanai, Huang and Wu) which derives quantum mechanics from a "dissi- 
pative Hamiltonian". This Hamiltonian was actually proposed earlier by Bateman, but in 
a classical context (Bateman 1938). This approach has the attractiveness of providing an 
exact solution, in essence because the classical equation of motion has an exact solution 
and formal quantization merely has the effect of converting the classical variables into op- 
erators. The solution is therefore formally the same, but must now be interpreted as an 
expression in operators. A critique on these models is that the commutator of Q and P 
fades away i.e. Heisenberg's uncertainty principle is violated. 

A slightly different way to quantize directly the classical equation starts from a "dissipa- 
tive lagrangian" that give raise to the classical equation of motion of the dho (see note 
[25] of Huang & Wu and see Papadopoulos). The idea here is that the propagator of the 
Schrodinger equation contains all quantum mechanical information of the system (like the 
Hamiltonian does) and that this propagator can be constructed from the Lagrangian. For 
a quadratic Lagrangian of a particle, the propagator can be given explicitly. Papadopoulos 
extends this idea to Lagrangians with not only time-dependent total energy, but also a time- 
dependent mass and developed a method to evaluate the corresponding propagator K by 
solving a certain differential equation for a matrix related to the Lagrangian. The dho fits 
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in this framework as its Lagrangian contains a time dependent coefficient for the velocity 
which can be interpreted as a time-dependent mass. 

A second approach uses an interaction Hamiltonian and applies perturbation theory. This 
approach tries to understand damping as the result of an interaction of two types of systems. 
One is a rather simple system (the undamped ho) that we want to describe completely, and 
the other is an extremely complicated system (the loss mechanism) which we want to de- 
scribe very incompletely. We are interested in the loss mechanism only to the extent that it 
affects the behavior of the oscillator and even in this respect the interest is limited: we con- 
cern ourselves only with the lowest order interactions which produce loss, since description 
of higher order interactions requires a more detailed knowledge of the loss mechanism than 
that provided by a single dissipation constant. Underlying the usual concept of dissipation 
is the understanding that the loss mechanism is affected only slightly by its interaction with 
the oscillator, whereas the oscillator may undergo large changes due to the loss mechanism. 
Thus the loss mechanism may be treated by a perturbation approach, but the oscillator may 
not. 

Common to the above two approaches is that the description of the environment to which 
the oscillator looses its energy is restricted to a dissipation constant (and possibly temper- 
ature). 

Approaches that construct an environment of the dho also exist. These in fact close the 
system, i.e. make a realistic or artificial embedding in a larger system that preserves energy. 
This way, Hamiltonians that describe a total, conserved energy can be obtained. 
An example of this line of thought is the so-called doubling the degrees of freedom ap- 
proach. In fact also this idea can be traced back to a Hamiltonian that was coined by 
Bateman, the so-called dual Hamiltonian (see Bateman). The idea is that the damped oscil- 
lator is coupled to its time-reversed image oscillator which absorbs the energy lost, so that 
the energy of the whole system is conserved (i.e. the system is closed). In fact, since the 
phase space of the whole system describes the dho and its image, the degrees of freedom 
are effectively doubled. Another way of looking at this is that the adding a time-reversed 
oscillator restores the breaking of the time-reversal symmetry. 

In the earlier attemps to elaborate this idea, difficulties arose such as that the time evolu- 
tion leads out of the Hilbert space of states, but later a satisfactory quantization could be 
achieved within the framework of quantum field theory (Feshbach et ai, Celeghini et al, 
Blasone et aJ). A more recent quantization based on the doubling the degrees of freedom 
approach was given by Banerjee et aJ. He obtained separate Lagrangians for the oscillator 
and its time-reversed image, although these Lagrangians were bound to have complex- 
valued parameters. In fact, these Lagrangians are complex-conjugated. This implies that 
resulting Hamiltonians are so-called "pseudo-hermitian", which nevertheless is sufficient 
to diagonalize the individual Hamiltonians via a generalization of the Dirac-Heisenberg 
treatment of the harmonic oscillator i.e. by forming ladder operators, a number operator, 
so that all eigenstates of the Hamiltonian can be obtained by applying a suitable ladder 
operator to a ground state. 

The doubling of degrees of freedom approach has the conceptual disadvantage that the en- 
vironment to which the dho is coupled (i.e. the image of the oscillator in which time is 
running backwards) is rather artificial. More realistic environments are constructed when 
forming unitary dilations such as in this thesis. In fact, this thesis presented a Bose quanti- 
zation of finite number of oscillators via an explicit dilation of contraction evolution T t . 
A Fermi quantization of a infinite dimensional phase space of the dho (i.e. an infinite num- 
ber of oscillators) is discussed in Latimer. The latter can be regarded as a logical extension 
of the doubling procedure. In fact, the contraction evolution T t of a single dissipative oscil- 
lator space can be "embedded" in a orthogonal evolution on the doubled phase space (see 
Latimer) and by embedding this space in an even larger space (of L 2 (R) type) a unitary 
dilation can be constructed. This has an advantage over the doubling the degrees of free- 
dom approach in case of an infinite dimensional phase space. The doubling the degrees of 
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freedom procedures only closes the system at each point in time and, if the phase space is 
infinite dimensional, at each point in time a rep of the canonical anti-commutations rela- 
tions must be made, while Latimer's method closes the system at once for all future times. 
Therefore, only one representing Fock space is needed and the vacuum vector is stationary. 
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Appendix A 

General 



A.l partial integration 

Definition A.l (Indefinite integral) 

Suppose F and f are functions whose domain contains a (closed bounded) interval [a,b] C 



Lemma A.2 (Partial Integration) 

Let f and g be integrable over a closed bounded interval [a,b] C R. Suppose F and G are 
indefinite integrals on [a,b] of f and g respectively. Then: 



Proof 

The proof is a straightforward application of Fubini's theorem. 

fa F S ■= J{Fg)[ a .j,](x)dx == f F [atb] (x)g [afi ](x)dx 

= f( F ( a ) + J ' f[a.x]{t)dt)g [a ^(x)dx = F(a)(G(b) - G(a)) + /(/ f[ a ^\(t)g[ a ,b](x)dt)dx 

NOW, use that / M (t)g [aM (x) = f [a>b] {t)g [a , b \ W l[a,*] (0 = f[a,b] (0S[a,i>] Wl[(,i] M = f[a.b\ {t)g[t, 

and apply Fubini (allowed since / \f [a h] {t)g [aM (x)l [t ^ (x)\dtdx < J \f [a>b] {t)g [afi] (x) \dtdx = 



Then F is called an ^indefinite integral of / on [a,b] iff 
P x f{t)dt = F{y) - F(x) (x,y e [a, ft] with* < y) 



F{b)G{b)-F(a)G(a)= / Fg+ / fG 




Ia\f\Ia\8\<-) 



to obtain: 



laFg ■= f(Fg) [aM (x)dx ==F(a)(G(b)-G(a)) + ff [aM (t)(fg [tM (x)dx)dt 
= F(a)(G(b) - G(a)) + Jf WM (t)(G(b) - G{t))dt 
= F(a){G(b) - G{a)) + G(b)(F{b) - F(a)) - J f [a , h] {t)G(t)dt 
=:F(b)G(b)-F(a)G(a)-J a b fG 



□ 



A.2 Integration formulae 



From contour integration, one obtains 
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[°° e -^+zx dx = JK/Ha , a g ( 0) oo) , Z € C) . 

J-oo V a 

Using partial integration 

Jx 2 e- ax2 dx = $-{x/2a)de- a}? 

= \-{x/2a)e- ax2 \ -f-e- ax2 /(2a)dx 

or by differentiating the above with respect to a, one obtains 

(ae(0,oo)). 



A.3 Joint centered Gaussian distributions 



LetX and Y be joint centered (i.e. with E(X) = E(F) = 0) Gaussian random variables i.e. 

X ~N{0,a 2 ), 
Y~N(0,b 2 ), 
Cov(X,Y)=pab, 

where p is the correlation between X and Y, and a and b are the standard deviations of X 
and Y respectively. 

The covariance matrix V of (X,Y) reads: 



V = 



a 2 pab 
pab b 2 



, with det(V)=a 2 fe 2 (l-p z ) 



If det(V) ^ (i.e. p ^ ±1), so that 

b 2 -pab 
-pab a 2 



det(V) 

the joint probability density of (X,Y) is given by 



p(x,y) 



1 



2jt v / det(V) 
1 



2jty / det(y) 
1 

27i v /ciet(y) 



exp 
exp 



1 1 



2 det(y) 



(fe 2 x 2 + aV-2patoy) 



exp[-(px 2 + a>' 2 -2yxy)] 



(A.l) 



where 



a : = 



2det(V) 



2det(y) 



y:= 



pab 
2det(V) 



Using that 



J — c 



-r^: 2 +^: _ z 2 /4r 



(re(o,oo), ze c), 



(A.2) 



A.4. Polarization formula 
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we calculate 

E(e ix+Y ) = r e ix+ yp(x,y)dxdy 

J — CO J —CO 



iet(V) /-= 



27tV^et(y) 



27t v /det(y) V P 



7C -4r 



e -{a-fmy 2 e {i+W)y d y 



1 



27ty / det(y) V P 
1 1 



jc - i 



71 e (l+''T/P) 2 /4(a-7 2 /P) 



2v^eT(y) v/^P - f 

" p - a + 2/y 



• exp 



a-y 2 /P 
p - a + 2/y 



4(ap-y 2 ) 



exp 
exp 



4(ap - y 2 ] 
b 2 - a 2 + i(2pab) 



(P^±l), 



where we have used in the fourth equality that a — y 2 /p > (since this is equivalent with 
1 > p 2 ), and the one but last and last equality that ap -y 2 = 1 /(4det(V)). 
Observing that b 2 - a 2 + i(2pab) = Var(F) - Var(X) + ;(2Cov(X , Y)), we can summarize 
this as 



E(e^ +F )=exp 



-Var(X) + Var(F ) + 2/Cov(X , Y) 



(A3) 



A.4 Polarization formula 

We want to show that the inner product is actually determined by the norm and, more 
generally, that a linear operator on a domain is actually determined by its "matrix-elements" 
on that domain. 

Lemma A.3 (Polarization formula) 

Let <D be a linear space and let B : 2? x 2? — > C be a map such that 

i. for each r| e 2? , the map ^h^B(^,T|) (t, e D) is complex-linear, 

ii. for each ^ e D, the map r) B(£,T|) (r\ e D) is linear, 

ThenB{^) = \Yl =0 {-i) k B{Z, + i%Z, + i k r\) for all ^ e ©. 
Examples A.4 

- If (-|-) is an inner product, then (^|r|) = iY.l=o(—') k \\^, + ' kr l\\ ■ 

- If A is a linear map, then (Z,\A\r\) = \ EjU(-i)* (5 + '^^^ + 

Proof 

Since Z? is additive in both arguments, we have 

B(^ + t 1 ^ + t 1 )-B(^-t 1 ^-t 1 ) = 2B(^t 1 )+2B(t 1 ^). 
Likewise, but now using in addition that B(r\,it,) = iB(r\,t,), while B(ir\,t,) = -iB(r\,£,) 

B($ + tt\£ + ty)-B$-tt\£-tt\) = 2W($,i\)-2W(%Z,). 
Therefore, B&T\) = \ Z 3 k=0 (-i) k B& + i*r\£ + 



□ 
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A.5 The operator expansion formula 

Let A be a self-adjoint and B be a bounded operator. Then 

exp ( ixA )B exp ( — ixA ) 

= B + ix[A,B] + M [A, [A,B}} + ... + ^ [A,[A,[A,...[A ,B]...]]] + . . . 

n times n times 

Let /(x) := exp(kA)Bexp(— ixA). Then 

/'(*) = ^gjjr = iA-exp(aA)Bexp(— ixA) +exp(aA)Bexp(— ixA) ■ (—iA) = [iA,f(x)]. 
From this, we derive next 

f{x) = d _m = m = [Wi igaj = [Wi 

and, with induction, 

/W W = [;A, = J/A, [iA, [/A, . . . [iA ,f(x)] . ..}}}. 

n times n times 

This means that in a (formal) Taylor series around x = 0: 

/(*) = /(0) +f(0)x + f (0)| +/ (3) (0)| + • • • 

(ix) 2 (ix) 3 
= B + ix[A,B] + ^~[A, [A,B}} + [A, [A, [A,B]]] + ... 

In particular, if [A,B] e C, then e' A Be~' A =B + i[A,B] or [e ,A ,B] = [A,B]e iA . 

A.6 The Campbell-Baker-Hausdorff lemma 

Let A and B be two self-adjoint operators whose commutator [A,B] commutes with both A 
and B (this is e.g. the case if [A,B] is a scalar). Then 

exp(;x[A+B]) = exp(ixA)exp(ixB)e~^ [mB l (A.4) 

Indeed, it suffices to show that the right hand side satisfies the differential equation that 
defines the left hand side (Stone's theorem). In other words, we need to check that 

d_ ^gixA e ixB e -\j?[iA,iB]^ _ ^ _|_ -g^ixA^xB g -^x 2 [iAJB] _ 

The latter is verified by direct computation. By the product rule, 

d_ ^ e ixA e ixB e - \x 2 [iA,iB] s j _ ^yixA g ixB ^^[iA.iB] _|_ £ ixA (jfiyixB e ~\i?[iA t iB\ 

+ e ixA e ixB (-x[iA,iB] )e~ M*.*] . 
Since [A,B] commutes with both A, B and thus e ,e , the latter term is 

(-x[iA , iB] )e bA e bcB e-^ 3? ^ . 
For the middle term use that the operator expansion theorem of above yields that 
e ixA {iB)e- ixA = iB + ix[A, iB] + + + . . . 

i.e. 

e ixA (iB) = iBe ixA +ix[A,iB]e ixA . 

Summarizing all terms then yields (iA + iB)e ixA e' xB e~^ x2 ^ iA '' B \ which was what to be veri- 
fied. 

(Alternatively, let g(x) := exp(xA)exp(xB), differentiate both sides with respect to x and, 
using the operator expansion formula above, show that ^j^- = (A +B +x[A,B])g(x). Inte- 
grate with respect to x like an ordinary differential equation.) 



A.7. Relation of the representation of in position and momentum representation 
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A.7 Relation of the representation of |\|/) in position and 
momentum representation 

Pj h generates translations in position, —Q/h translations in momentum (see e.g. [17]). 

A.8 Coherent states are eigenstates of A 

Recall that he operator e~ ltP l h implements a position translation over t to the left (e~" p / h \\f(q) = 
\\f(q + t)) and e' s Ql h implements a momentum translation over s to the left: e ,s Q/ h \\f(p) = 
\)f(p + s). 

In order to show that 

(a(Q -t) + i(P-s))\\]f) = =► [aQ + iP] e - is Q +itP | ¥ ) = 

IVo> 

and, vice versa, that 

[aQ + iP] |Yo) =0 =► (a(Q -t) + i(P-s)) e isQ - itP | ¥o ) = 

W 

it suffices to show (for the first statement) that 

[aQ + iP]e-*Q/n+i' p /n = e -isQ/n+itP/h ^ Q _ t - ) + ;(p_ 

^ e -(-isQ+itP)lh{ a Q + iP ] e (-isQ+itP)/h = a (g _ t ) + i(p _ y ) ) 

which can be verified by applying the operator equality e A Be~ A = B + [A,B] that holds for 
A, B with [A,B] eC: 

[{isQ-itP)/h,aQ + iP] = isa[Q,P]/h+{-it)i[P,Q]/h=-as + -it 



A.9 Inner product of coherent vectors 

(Coh(w)|Coh(z)> = J dq(Coh(w)\q) (?|Coh(z)) 

J V V%h ) \/%h 

= C 1 € 1 . e -(q/V2h-w) 2 e -(<,/V2h-z) 2 dq 

Vnh J 

e -^(\w\ 2 -W 2 ) e -^(\z\ 2 -Z 2 ) r g 2 2q(z+W) ,2,-2+, 

= = / e T v/2» [Z+ >dq 

= vTn e V^ exp lvV) i^) 

= e -^\«\ 2 --v 2 ) e -{{\z\ 2 -z 2 ) . e -(z 2 +w 2 ) e2 -(z+w) 2 
= e -^ 2 -« 2 ) e -{(\z\ 2 -z 2 ). e -{z 2 e -^ 2 ) e z« 

= e -2-\ w \ 2 e -h\A 2 e w 



A. 10 Coherent states in the Schrodinger representation 

We determine the form of coherent states in Schrodinger representation by solving the 
eigenvalue equation A |\|/) = z |v(/) in position representation and then evaluate the action of 
the Weyl oper ators on them. This action takes a particular simple form if h/2 = raco = 1 
i.e. ^/nm/2h = 1 /2, see (A.7), (A.8), (A.9), . 
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A.10.1 Coherent vectors in position representation 

Since 



/ mco / iP \ I mco 

— \Q + — - ) = \hrQ- 



1 



-jp 



2h V & ' m&J V 2ft \/2ftmco 
the differential equation for V|/(g) = (q\y) in position representation reads: 
«<?¥(<?) + jEiV(q) = zMq) with a := y/^j[, 



i.e. 



\^(q) = -(2a 2 q-2az)y{q). 
Integrating this equation, we obtain 

y(q)=Nexp(-(a 2 q 2 -2azq)j (?et), 
where the normalization constant N satisfies 

\N\- 2 = J |exp(-(ocV -2azqf)\ 2 dq = J exp Re (-(2a 2 q 2 -4azqfjdq 



j exp (-(2. 



a 2 q 2 - 4aRe 



(z)q))dq 



n /[-4aRe(z)] 2 



:<x-' eXp ( 4 2a 2 i 



) 



For the convenience of N taking a a simple form, we admit a phase factor exp ^/Im (z)Re (z)j , 
^fexp(-I(|z| 2 +z 2 )). 



so that N 



All in all, we then have 



Coh(z)^^exp(-i(|z| 2 +z 2 ))exp(-(aV-2az< ? )). 
A.10.2 The action of W(z) on the vacuum state 

From (A. 5), 

Coh(0) = -{/^ exp(-aV), 
and from (1.13), 



(A.5) 



W(z)^e-^'S s T t 



so that 



W(z)Coh(0)=exp(-f S f)5. s r r ^ 2 f -exp(-aV) 

= ■ exp(- ^st)exp{isq) exp (-a 2 (q - f ft) 2 ) 
= ■ exp(- yif + isq + -a 2 q 2 + 2a 2 thq - a 2 t 2 ti 



On the other hand, if we define 



Z - := \l2( 
then from (A.5) 



/mco 



tah+— (z = t + iseC), 
2a 



(A.6) 



A.10. Coherent states in the Schrodinger representation 
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Coh(z.) = [Coh(f + w)] fH 



■tah,si->s / '(2a) 



2 -^- exp (f 2 + /if )) exp (-a 2 g 2 + latq + 2iasq)^j 



/«(oS,sMs/(2a) 



exp (f 2 a 2 ft 2 + Zi^)) exp (-a 2 g 2 + 2a 2 htq + isqjj 
Concluding: 

W(z)Coh(0)=Coh(z.). 

A. 10.3 The action of W(z) on coherent vectors 

By the Weyl relation W(u)W{z) = e'^^W (u + z), the action of W(u) on Coh(z) can 
be evaluated as 



W( M )Coh(z) =W{u)W(z')Coh(0) =e- i 'i lm{uz ' ) W{u+z')Coh{Q) 
= e -^^-) Co h( u .+z). 



where 



z* = ^— + Has = \\% ( * — + is^/mu>\ (z = t + is € C) 
ah V h\Jm& > 



(A.7) 



(A.8) 



i.e. z^z, and z^z' are each others inverses. 

Concluding, the Weyl operators reshuffle the coherent vectors and add a phase factor. A 
particular simple form arises if mco = h/2 = 1: W(t<)Coh(z) = e~'7 Im (" z >Coh(M + z). 

A. 10.4 Coherent vector algebra 

The vacuum state expectation of the Weyl operators are given by 
(Coh(0)|W(z)|Coh(0)) = (Coh(0)|Coh(z.)) 



= J (^J^ e -o?q^ ^^2^ e -i(\z.\ 2 +^) e -(aV-2az.g)^ dq 
= y / 2Z e -i(lz.l 2 +4) J e -(2a 2 q 2 -2az.q) dq 
= \J?fe-^ z -\ 2+ & • y^exp((2az.)7(4-2a 2 )) 
= ex P (-i|z.| 2 ) (zeC), 
which enables us to calculate the inner product between two coherent vectors as 

(Coh( M .)|Coh(z.)) = (W( M )Coh(0)|W(z)Coh(0))S 
= ^Coh(0)|W( M ) f W(z)|Coh(0)^ = (Coh(0)|W(-«)W(z)|Coh(0)) 

= e-'^ Im (Coh(0) | W(z - u) |Coh(0)} = e^ lm l z '-"- 1 2 

= £ ,4lm(i; z ) £ ,-^|z.| 2 -|«.| 2 +Re(77.z.) = e -\\z.\ 2 e -\u.\ 2 e u,z. 



(A.9) 



(A. 10) 



since if we write u — y + ix and z — t + is, then 

u.Z. = yj§ (yVm©- j=) \[\ (tVHim+ - | (bf (mco) + ^] + i[ys - fx]) , 



i.e. |lm(i<z) = Im(w.z.). 
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A. 11 Exponential vectors 

Exponential vectors are no more or less than "unnormalized" coherent vectors i.e. 

e(z) = e ^l 2 Coh(z) ( = 5) {/^exp(-iz 2 )exp(-(aV-2az< ? )) 

V 7t 2 V ) (AU) 

= y ex P ( 2a ^ - _ a2q2 ^) ■ 

From the above, the inner product between two exponential vectors is then (e(u)\e(z)) = 
In particular, if mco and h/2 = 1 are taken to be 1 (so that a = j), then 

£(z) = ^fe eXp ( Z9_ ^ 2 -?'? 2 ))- 

A.ll.l The coherent (exponential) vector are total in ^schro = L Z (R) 

If / e L 2 (R) satisfies (e(z)|/) = for all z e C, then 

= <E(-&)|/) = /" (e~^ f{q))e^dq (t e R). 

2 2 

In particular the Fourier transform of g e - a ~<? vanishes and, by the injectivity of 
the Fourier-Plancherel transform on L 2 (R), therefore 

e~ a2q2 f(q) = in L 2 (R) i.e. / = in L 2 (M) . 

A. 12 Coherent states in the Heisenberg representation 

To calculate the action of W(z) on coherent vectors, we write the coherent vector Coh(t<) 
in terms of exp(uA^) and we write W(z) in terms of A and A T . 

Step 1: Writing W{z) = W(t + is) = exp(wg - HP) in terms of A andA f 
Using the relations (1.6) between A and A 1 " with P and Q, 

Q = \[^ {A +Ai) = h {A and p = \fW {A ~ A+)// = ah{A At)// ' 

we have that 

isQ - itP = (A + A f ) - ;raft(A - A*)/i = z.A + - z.A 

with z. = f aft + ^ as above. 
Therefore, 

W(z)=exp(z.A + -z.A) (z. := + -fa) = tah+ &) (A.12) 

Step 2: Coherent vectors 

From the relation A 1 " |n) = Vn+ 1 |n + 1), it follows that 

CO f| 

e(z):=exp(«A + )|0H E^l")- 

The "exponential vectors" e(z) are closely related to the coherent vectors Coh(z), since 
Coh(z)=e-W/ 2 e(z). 

Step 3: The action ofW(z) on the vacuum state |0) 



A.12. Coherent states in the Heisenberg representation 
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If A, B are operators for which [A,B] commutes with A and with B, then the Baker-Campbell - 
Hausdorff formula e A e B = e A + B e i\ A M holds. Using that, we derive 

W(z)exp( M A t ) = exp(z.A t -z.a) exp( M A t ) = e^'")/ 2 exp((z. + M )A f -zS) 

= e -2.")/2 e -(z.+«)?. eX p((z. + M )A f ) exp(-z.A) 

= e -^"-2k-l 2 exp((z. + M )A + ) exp(-z.A) . 

Since exp^— z.A^ = 1 + z.) n A") and A (and its powers) annihilate |0), we have 
exp^— z.A^ |0) = |0). Using the Baker-Campbell-Hausdorff lemma, 

W(z)e(«) = exp^z.A 1 " -z.a) exp( M A f ) |0) = e~^*"~^' z *' 2 exp^(z, + t/)A^ exp(-z.A) |0) 
Therefore, 

W(z)Coh(0) = W(z)e(0) =e-2W 2 e ( z .) = Coh(z.). 

Step 4: The action ofW(z) on coherent vectors 
Using the Weyl relation, it follows as in (A.7) that 

W( M )Coh(z) = e- ! 7 Im ( s ')Coh( M . +z). 

The vacuum state expectation of the Weyl operators can be calculated as 



(Coh(0)|W(z)|Coh(0)> = (Coh(0)|Coh(z.)) = (e(0)\ e -^ e(z.) 

= e4W 2 (zeC), 
and as in (A. 10), the inner product between two coherent vectors is then verified to be 
(Coh(«.)|Coh(z.)> =e-H z 'lVl"'lV' z v 

A.12.1 Coherent and exponential vectors are total in #H e i S = £ 2 (N) 

As expected by formal differentiation, we have 

and this can be verified as follows: inductively define for / G R 

e(°)(0 := e(r) = Y ±= \n) = (l, -i=, 4=, 4=, - • ■)> 
W W „=oV«! Vl! V2! V3! 

e(" + ')(Q : =Um eW W- eW ^. 

S - t 

and prove by induction that 1 1 e w (f ) - VnT | n) \ \ 2 — > if f -> for all n e N. 
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A.13 L 2 (Rr"^L 2 symm (W) 

(See [14]) We treat the case n — 2; the case n > 2 is a straightforward generalization thereof. 
Step 1: L 2 (R) <g> L 2 (R) ~ L 2 (R 2 ) 

Let 4>:L 2 (R)(X)L 2 (R) — > L 2 (R 2 ) be the linear map that sends / <g> g to the function (f i , f 2 ) 
/(fi) ■ g{t 2 ). Let #b be the linear subspace of L 2 (R) consisting of the linear combinations 
of indicators of measurable sets, 

#b = {Li CXilfij : n G N, a,- G C, B t measurable }. 
Then H§ is dense in L 2 (R) , and as consequence itfo <8> #b is dense in L 2 (R) ® L 2 (R) while 
its range under <t> is 

#i := {Lf^'lA^B, : N e N, a,- G C, A,-, B, measurable } 
which is dense in L 2 (R 2 ). 

Moreover, <t> is isometric as map from yfrj <g> tt§ onto tt\ , since 

||4>(l Bl <g> 1b 2 )|| 2 = Jo So l 2 Bl (t i ) • 1 2 B2 {ti)dh dt 2 = KB i ) • HBi) = 1 1 1b, 1 1 2 1 1 U 2 1 1 2 ■ 
Therefore, <t> extends to an isometric map from L 2 (R) ® L 2 (R) onto L 2 (R 2 ) . 
Step 2: L 2 (R) <g>. s L 2 (R) ~ L 2 ymm (R 2 ) 
The map P : L 2 (R 2 ) -> L 2 ymm (R 2 ) defined by 

(PF)(t u t 2 ) := f ha)+ffe.'i) (F e L 2 (R 2 ) , ( fl ,f 2 ) G R 2 ), 
is surjective i.e. every symmetrical L 2 -function on R 2 can be obtained in this way. 
Concluding, L 2 (Rf* ~ L 2 (R 2 ). 



A. 14 Positive definite kernels 

The verification that the (infinite) sum of positive definite kernels is a positive definite 
kernel is straightforward from the definition. 

We now show that the point wise product of two positive definite kernels is positive definite. 

Let K\ and K 2 be positive definite kernels on S, take n G N, s\ , . . . ,s n G S, and Xi, . . . , X n G C. 
We need to show that Y,ij'XjKi(sj,Si)K 2 (sj 7 Si)Xi > 0. 

First, observe that the n x n matrix with elements A := K\ {sj,si) is positive definite, since 
K\ is positive definite (X'AX > for all X G C). As a consequence, there exists a matrix B 
such that A = B 2 and B' = B. In particular, A Jt = YX=\B jk B ki = ££ =1 %B fa . 
Therefore, 

Y,kjKi(sj,Si)K 2 (sj,Si)Xi = Y,kjAjiK 2 (sj,Si)Xi = £ {XjB kj ) ■K 2 (sj 7 s i ) ■ (B ki Xi) 

hj hj hj,k 

=j:{xfK 2 (sj, Si )xfn>o, 
k } 

since K 2 is a positive definite kernel, and (Xf^ :— B ki Xi) n i=l are n sequences of scalars. 



A.15 Projection valued measures 

is a dense subset Given a projection-valued measure P on 9{ and r| G H , we have 
that r)„ := P(-«,n)r| converges to r| since U„(— n,n) = R and P(R) = /. Furthermore, for 
fixed n and Borel set B: 

P n „(B) = (T\n\P(B)\r\ n ) = (P(-n,n)T)\P(B)\P(-n,n)T]) = (r\\P{-n,n)P{B)P{-n,n)\r\) 

= (Tl|P(Bn(-n,n)h)=P T1 !(_„,„) (B). 
This means that 

/ R |x| 2 P n „(Jx) = f |x| 2 P„(<&) < » 2 P^((-»,«)) = n 2 (Tl|P(-n,n)h) < « 2 ||il|| 2 < ~. 
Hence ^ := G ^ : J R \x\ 2 (^\P(dx)\^) < 00} is dense. 



A. 15. Projection valued measures 
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Hp is a linear subspace Furthermore, for fixed B, P := P(B) is a projection and the 
mapAf(^,r|) := (£|P|t|) is an inner product except for the fact that = implies that 

^ € P# (i.e. the null space of P). This however is still enough to show with the standard 
argument that N satisfies the Schwartz inequality: |JV(4,T|)| < ^ N (£,,£,) ^N(r[,ri). In other 
words, | (t,\P\v\) | < a/(^|P|^)-\/(ti|P|ti). The product of the latter is of course less than the 
product of the largest i.e. 

I &\p\t\) I < max^TO , (n|Ph)} < + (nlPh) • 

As a result, we have that 

(§ + ti|P(B)|^ + ti) <2^\P(B)\^) + (Tl|P(B)h)) (B G Borel(M)). 
so if ^,r) e then ^ + r| G and is a linear subspace. 

The formula = / x(£,\P(dx)\Q defines a linear operator on Hp For ^ e Hp, we 

have that 

hence the right hand side is well defined on Hp. 
For^,r) e Hp, (^|0|r|) is defined by polarization: 

<^|o|ti> := ^£(-0*^+^101^+^) = \ t(-0* / **W^)- 

To prove that (^|(9|r| + Q = (^|<9|r|) + it therefore suffices to show that the measure 

?ELo( _i ')*^+(*(ti+0 
coincides with the sum of two other measures, viz. 

iELo(-0*JW and ?£LoH)^+«*c- 

That however is clear since for each Borel set B 

llLo(-0^(ri + Q( B ) = 3i:Lo((-0 fc ^+^(n+QTOI^+^(Ti+Q> 

= (^|P(B)|(T 1 + Q) 

by polarization, and likewise 

iELo(-0*^( B ) ^ &\ p ( B M and ?ELo(-0^ + ^(B) = <W)IQ • 

O is self-adjoint The argumentation proceeds smoother if we realize that 

coincides with the notion of integration of x with respect to the complex-valued measure 
B^(l\P(m) i-e. (§|Oh) =/x<^(dx)h). 
This way, we see that for ^,r| e = ©o 



ftlOfc) = | *<Tl|P(dx)|$) - / x(T]\P(dxm = J x(P(dx)^) = J x{Z,\P{dx)\l\) 

i.e. (^|<9|r|) = (0£|t|) which shows that D * 3 2?o and that the restriction of O* to D 
coincides with O. We now show that actually, 'Do* = ®o- To this end, suppose that £ e 
© *\£> i.e. 

(SlOfo) = <0*$fo) forallTie^., 

but /^^(^P^jc)^) =oo. In particular, %± 0. Since T) n := P{—n,n)0*\ € #jp = £>o (see 
the argument why Hp is a dense subset), we have by taking r| = r|„ in the above equation 
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(£,\0\r\ n ) = (0^h„) (allneN). 
For n — > oo the right hand side converges to (0*^|6>*^) (since 0*t, e H and r|„ — > <9*c;) and 
is therefore bounded. The left hand side evaluates to 

/OO /"DO />ft 

x&\P(dx)P(-n,n)\0*Z > ) = *<$|P(<k)0*$) 
-oo J oo J 77 

x{OP{dx%\$ = x(%\0\P{dx)%) = / x/ x'(^|P(Jx')|P(rfx)^) 

= r x ^i,\p{dx)^)= r j(%\p(dx)\%), 

J —n J — n 

and tends to °° if n — > oo. Contradiction. Thus, there cannot be a £ e £>q* \ ©o. 



A.16 Proof of lemma 1.28 

For neN, let □„ be the square in C ~ R 2 with corners at the points (±2", ±2"). Thus, 

the length of □„ is 2 • 2" and DoCDiCHC.... We can think of □„ as filled up 

with small squares of length {\) n , which have corners at the points (j) n • (k,l) with k,l — 

0,±1,±2,...,±2 2 ". 1 

Let <£ n be the collection of these squares. 

What we will use of this construction is the following: 

i. □„ is an increasing sequence of squares which fill up M ~ C, 

ii. □„ is the union of all squares R in fP„, 

iii. for each R <G (B n , the distance of two points within R is less than \/2(|)", 

iv. if n > m then each square R e <£ m is the union of squares in 2>„ that lie within it. 
For n € N define the step functions 

T n := £ X 2 (/?)-p(z ff )-W(z R )= /"[£ p(z«)W(z ff )l ff (z)]x 2 (*), 

where zr G R are chosen such that 

|p(a)|||W(a)||W< / |p(z)|||W(z)||X 2 (*) («€!»,), (&) 

JR 

(for z/j we could, for example, choose the point in R where the continuous function z i— > 
|p(z)| || W(z)| attains its minimum over/?). 

Then T„ e ® (#), and we prove that 

i. for all £ e # , (r„^)„ G N is a Cauchy sequence in therefore its limit T% := lim„ Tj^ 
exists; 

ii. for all ^e^, (T\\T\)=f c {T\\W(z)%)v(z)Wdz). 

Step 1: for each the map z i-> p(z)W(z)£ from C into H is continuous 

Simply, because z i— > W(z)£ as well as z i— > p(z) are continuous. 

Step 2: for each the sequence (T„^)„ in H is Cauchy 

Let "% G H and m > n > N in N. We split up T m £, — T n ^ in three terms: 

T m \-T n \= E MR)P(zr)W(z S )1- £ X 2 (/?)p(z«)W(z«)^ 

ffes> m *ea> B 

= £ ^^pywyu £ x 2 (fl) P (z s )w(z^ 

scc\rjjv Rcc\n N 

+ ( I ^(«)p(z#fe)- E X 2 (fl)p(z*)W(z*)k 



'The centers of these squares are (j)^ 1 ■ (±(2k+ 1), ±(2/ + 1)), with k,l = 0, 1, . . . ,2 2 " - 1. This are in total 
2 2 " ■ 2 2 " = 2 4 " squares in each quadrant, so 2 4 " +2 squares of volume 2~ 2 " in total (which equals the volume 2 2n+2 

of a). 
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and we show that the norm of each of these terms vanish if m > n — ► °°. Actually the first 
two terms will vanish because of the integrability of z i— > |p(z)| ||W(z)||, and the last term 
will vanish because of the uniform continuity of z i— > p(z)W(z)£, on compacta. 
To make this precise, let e > 0. 
First we take N so large that 



[ \ ? {z)\\\W{z)\\Mdz)<t/m\\) 

Jc\a N 

(Lebesgue's dominated convergence theorem). 

With this choice of N we have for the norm of the first term: 



(#1) 



£ X 2 (R)p(z R )W(z R )Z, 

Rcc\n N 



R€(P m 

k~cc\n N 



< m [ v \ P (z)\\\w(z)\Mdz) < m f P (z)\\w(z)\Mdz) 

Ja m \a N Jc\o N 
<e/3, 

where the second inequality uses the choice of Zg, see (&). 
Likewise the norm of the second term is < e/3. 

Now, we consider the third term. Since z i— > p (z)W (z)i; is continuous, it is uniformly con- 
tinuous on the compact set and so there exists an M > N in N such that for all n > M 

z,z! e D N , 
z,z' Gflwifh 

Suppose now that both m and n are larger than M. 

Since m> n, each square R e <P n is a union of disjoint squares R e <P m . Further, since 
m > n^> M >N,Dn is the union of all R e "P n with C and likewise D N is the union 
of all R e rP m with C Thus, 

£ WpfeMz^- £ X 2 (fl)p(z«Mz^ 
Rca N 



p(z)W(zK-p(z')W(z%\< 



3X 2 (a N ) 



(!) 



fi6fP„ 



Rca NRcR 



Rca NRcR 



< £ £ ^W||P(^)W(^)-P^)W( Z ^| 

«CDjv ScR 



rco NrcR 
<e/3 



Rca N 



X 2 (R)=X 2 (n N ) 



Step 3: for all lim m (^7^) = / c p( z ) (l\\W{z)%) 

Let t,,r\ e # . First we mention thatz — > p(z) (r||W(z)^) is integrable, since it is a continu- 
ous function (by the strong continuity of z — > W(z)) dominated by the integrable function 

z^p(z)llw(z)Hliyihii. 

Next, we split / c p(z) (r\\W(z)£,) X 2 (dz) - (r\\T n %) in three terms: 
f c p(z)(r\\W(z)l)X 2 (dz)-(r\\U) 



L 



p(z)(r\\W(z)l)l 2 (dz)- £ fc(*)p(»)ft|W(»)$> 

C\Djv /; Gi > n 

RcC\Da? 



+ 



/ 

Ja N 



p(z)(r\\W(z)l) £ p(z K )(T!|W(z^)l«(z) 



Rep„ 
Rcn N 



X 2 (dz) 



84 



Appendix A. General 



and show that each of them vanishes if m — > °°. 

To this end, let e > 0. First take N so large that / c ^ Djv |p(z)| \\W(z)\\ X 2 (dz) < e/(3 ||£|| ||tiJ|. 
Then the first term is in absolute value < e/3. The same holds for the second term, since 
the absolute value of each summand is majorized by ||r||| |p (^) I II^WIIM^z) by 
definition of the zrs. To handle the third term, we use (as in step 3) the uniform continuity 
of z — > p(z) ||W(z)|| to find anM > N so large that for n>M: 

z,z' e U N , ) 

z,z'eflwith \\p(z)W(zK-p(z')W(z%\\< , , . (!2) 

Re¥„ J 3||r||| A 2 (U/v) 

Next, we rewrite the integrand of the third term as: 

p(z)(i\Wz)Z,) ln N M E p(z R )(r\\W( ZR )l)l R (z) 

RS1'„ 

Rca N 



E [p&mm -p(zr)(t\\w(z R k))mz 



R£l'„ 
RcD N 



£ (Tl| [p( Z )W( Z )^-p(«)W(^])- 1,(2) 



RGS>„ 



which implies that for n > M the absolute value of the integrand is bounded by e/(3X2(0^)) 
(use (!2)), and hence the third term is bounded by e/3. 
Concluding, ^\T%) = lim„ (i\\T n %) = J c p(z) (l\\W(z)%). 



A.17 Proof of lemma 1.29 

We only discuss the claim for Q. The claim for P is proved analogously. The existence 
of a sequence fk with the desired properties follows from [16, p. 93-94]. As the Fourier 
transform maps the space of Schwartz test functions S into itself bijectively [16], P^ a b ^ = 

J c fi c (s)W(is)ds exists strongly, lies in the von Neumann v3\£ , and for any state \|/: 

Taking the Schrodinger representation as cyclic vacuum representation of the CCR, we 
have that Q is the multiplication by q (so that W(is) = e is @ is the multiplication by e' sq ) and 
(q\\\t) = \)f(q), so that 

(v\Ptajb)\v) = f-~Ms) (Me isQ \v)ds = S-~S-Jk{s)e isq \V{q)?dsdq 

= J-~fk(q)W(q)\ 2 dq ^ I b a W{q)\ 2 dq 

which equals (\\t\P(a,b)\\\t) i.e. the probability that a measurement of Q gives a value in 
(a,b) when the system is in state \\t. Hence Pg(a,b) lies in v?\£ too. 

If / = EJU a-iHaubty then /(g) = «il( fl ,A)(2) = I"=i ^PqM) lies in v5\£ ■ 
Finally, if / S L°° (R) , then there exist a sequence of simple functions such that 1 1 / — fa \ | «, 
(see [16]). As a consequence of the dominated convergence theorem, the corresponding 
sequence (of multiplication operators) will converge to (the multiplication by) /, and hence 

{/(e):/er}c^. 



A. 18 Proof of lemma 1.30 

We show that 

- 1 1+1 Izl 2 

Takes e (-1,1) and let 



A.19. The uniform derivative of a absolutely summable operator-valued series 



85 



P«:=5^«P(-2**SI*I 2 ). 
and let W(z), z € C, be the Weyl operators on C. Then that the conditions of theorem 
1.28, are satisfied, so the strong integral on the right hand side makes sense. To show that 
it coincides with s N , we evaluate its "matrix elements" with respect to the total subset of 
exponential vectors. 
Take «,veC. Then by theorem 1.28 

(e(«)| [f c P(z)W(z)Mdz)} \e(v)) = f c P(z){e(u)\W(z)\e(v))X 2 (dz) 

= J p(z)(E(u)\[e-^-^ 2 ]e(z + v))\ 2 (dz) 

= f e^p(z)e-^ 2 e-^ + ^X 2 (dz). 
Jc 

Now 

P(Z> -s) n(l-s) 

Further, writing z = x + iy, we calculate 

- (z\v) + (u\z) = -zv + uz 

= [— v + 77]x + i[v+ u]y . 

Using this and the fact that for w € C and a e (0,°°) f™ aa exp(—ax 2 + wx)dx = a&xp(w 2 /Aa) 
we resume calculating 

(e(u)\[J c P(z)W(z)X 2 (dz))\e(v)) 

= - —— e-^ + ^ v+ ^ x dx- e-^+^+^dy 

(1—S)KJoo Joo 

= e <»|v) exp [ • -477v ] = e'™ v . 
Since e„ = ^-e(x), 

(e„| (^p(z)W |e m ) 

we see that f c p{z)W(z)X 2 (dz) is in vJ\£, and that (e„| J c p(z)W(z)A, 2 (t/z)e m ) = s"8 nm i.e. 
/ c p(z)W( z )^(dz)=^. 

A.19 The uniform derivative of a absolutely summable operator- 
valued series 

Let f(z) = Y.k a kZ k be a series that is absolutely summable on an interval containing 0, and 
let P k be a bounded series of operators. Since f(z) and its derivatives are all absolutely 
summable (f(z) is a convergent power series on a disk around in C and all its derivatives 
exist on that disk and are hence absolutely convergent too), the series 

f(z) := £ \a k \s k P k 
k=i 

f^(z):='£k---(k-n + l)\a k \s^P k , 

k=n 



86 



Appendix A. General 



and 



k=i 

AW ( s ) : = £ jfc • (jfc - 1) • • • (k - n + \)a k s k - n P k 

k=n 

exist. 

After defining 



F^{s):=^a k s k P k 



k=l 



s :-- 



■lim 



fW(f)-F w W 



<-»s f — s 
we can prove by induction A^"' = F^ for all n. In fact, 

F(")(t)-F^(s) , ,i 



f-5 

CO 

£ k---(k-n + l)a k 



fk—n ^k— n 



t-S 



k=n+l 

oo 

<supj|P„||- £ fc---(fc-n + l)-|Ot|- 



fc=«+i 



(Jfc -«)/""- 1 P k 

' -{k -n)s k - n - 1 



sup||P„| 

n 



A.20 The count operator d^F (P) 

Given an orthogonal projection P on ?C, we can define a self-adjoint operator df (P) on 
via 

e iM7 {P).- 7 ( e «'W>) (A, e R). (A. 14) 

Indeed, P is self-adjoint and therefore e'^ p is a unitary map %_ — > = ||/|| for all 

/€ As a consequence, J (e 1 ^) is unitary map J J (^),as ||/i <g> ft®...® /„|| = 

1 1 /ill II/2II ' 1 1 /nil for all /1, . . . ,/„ e 3C. As one-parameter semi-group, X<-^ T (e' Xp ) it has 
a generator <f J (P) on J (ac). 



A.21 Representations of T(%i ) are unitarily equivalent 

Suppose anc l i^'x^x) ®t Q two representations of T(%,) over # and re- 

spectively, then both of them contain a cyclic vacuum representation of the CCR over 3C, 
hence there is a unitary map U : H' -> H such that W'(f) = U~ l W(f)U for all f ^ %_. 
In other words, the linear map S {H) — > S (#"')> A i-> A' = U~ l AU sends {W(/) : / G 3C } 
bijectively to {W'(f) : / € 3C}. Moreover, A A' is strongly continuous and therefore 
extends uniquely to the strong closures of {W(f) : f € %_} and {W'(f) : / € ?C} i.e. it is a 
strongly continuous, bijective map from to . 



A.22 A minimal dilation of the classical damped harmonic 
oscillator 



Proof 

Existence: Then for all t > we have 



A.22. A minimal dilation of the classical damped harmonic oscillator 



87 



(v\R t v) = [ e^+^e^-^-'hjdx 

J— oo 

= e (-ri+to)< f° e 2v 2 ydx 

J — oo 



_ e (-Tl+/CO)t _ 

Furthermore, the linear span of the vectors R t v is dense in L 2 (R,dx) = L 2 (R,2ydx). To 
see this, it suffices to show that {R t v : t <G R} x = {0}. When we realize that the Fourier 
transformation extends to a unitary map J : L 2 (R, dx) — ► L 2 (R, dx) 

rN 

(J g) (k) := L 2 — lim / e~ ,to g(x) afx for almost all k e R 

N~>°°J-N 

that satisfies 5" T t T~ l — M e -i K , we see that we can just as well prove that 

{0} = {7 T t T - 1 J v : t e M.} 1 - = {e^ 7 v : t e K} x . 
To this end, let f e L 2 with / _L e" x J v for all t e R. Then 

/OO /"CO 
f(k)(e i,x fv){k)dk= / f{k)e- itk (fv){k)dk forallreE. 
-oo J — oo 

Apparently, the Fourier transform of f ■ j~ v vanishes, so / • J v = a.e.. Since ( J v)(fe) = 
11 _ i . ( 1 m+<:) ^ for almost all keR, this means that / = in L 2 (R). 

Uniqueness: We claim that any minimal unitary dilation (SC,V, Ut) of (G) r >o induces -in 
the notation of theorem 3.2- an Kolmogorov dilation {9{ = %i,V{t) =U t v) of the follow- 
ing, positive definite, kernel on R: 

( e (-n+i<o)(t-s) if t > St 

K{s,t) - | e (-i-to)(,-») iff <,. 

Before establishing this claim, let us first see how it implies the uniqueness: let (2C,v,t/ r ) 
and v,U t ) be two minimal unitary dilations of (C f ) f >o, then -assuming the claim- they 
are both Kolmogorov dilations of one and the same kernel. By the Kolmogorov Dilation 
theorem 3.2 there exists a unitary map U : such that for all t el: UV(t) = V(t) i.e. 

UU,v = U t vJov all t e R. Taking f = we have Uv = v so that for all t , t' e R : UU, (U t >v) = 
UU t+t iv — U t+ frv = U t U t iv = U t U (U t /v). Since the span of {U t v : t e R} is dense, we may 
conclude that U t — UU t U~ l . So: any two minimal unitary dilations of (Q)t>o are unitarily 
equivalent. 

We now establish the claim: let (sc,v,J/ f ) be an minimal dilation of (C t ) t >o- Obviously, 
(s,t) i ► (U s v\U t v) is a positive definite kernel on R with Kolmogorov dilation (5/ = Ki,V{t) - 
U t v). Using the fact that (U t ) t eR is a unitary one-parameter group, we see that for t , s € R : 

( e (-i\+ita)(t-s) if t>s 

(U s v\U t v) = (U s v\UsU,- s v) = ( V \U t . s v) = |___ = e( _ ii _ to)(j _ () ^ ^ 

So K(s,t) = (U s v\U t v) is a positive definite kernel and V : R — > 3C , V(f) = f/>v is an 
embedding such that (V(j)|V(0) = K(s,t) and yV(R) = XL. □ 
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Appendix B 

The 2 nd quantization functor T 



B.l Complete Positivity 

When working in representations, there is a simple criterion for complete positivity: 
Theorem B.l (criterion complete positivity in representations) 

Let j?i C <B(H\), At. C ® (#2) be Vb/i Neumann algebras, and letF be a linear map A\ — » 
j?2- Then T is completely positive if and only if for all N £ N, aJJ £1, . . . e #2, and aii 
Ai,...,A„ e ili: 

I^ =1 (^|r[A^]^) >o. 
B.2 Strong operator topology 

Lemma B.2 

Let X be a dense subset of a be a Hilbert space H . Then F a — » F in <B(9{) strongly if and 
only if || (F a - F)t\\ for all^e X and sup a ||F a || S (^) < °°. 

In particular, if H is separable, then the strong operator topology is metrizable (and conti- 
nuity coincides with sequential continuity) on operator norm bounded subsets of <B (#"). 

<3 Replacing F a by G a := F a — F, we have sup a ||G a || < 00 if and only if sup a ||F a || < °°. In other 
words, we may assume that F = 0. 

If: suppose that | F a i; 1 1 — > for all % <E X , and sup a 1 1 F a 1 1 B 1^ \ < °o. We now prove that 1 1 F a r\ \ \ — > 
for all n e # . To this end, let T| e and e > 0. Take a i; e X with ||r| — £|| < e, and take oc(e) 
such that ||F a £|| < e for a >~ a(e). Then for a >- a(e): 

||F«ti|| < ||F„(ti-5)|| + ||F„5|| < (supJlFxH) h-^||+e 

< (sup a ||F a || + l)e 

Since e can be chosen arbitrarily, it follows that ||F a T||| — ► 0. 

a 

Only if: Suppose ||F a T||| — > for all r) e 9l . Then, sup a ||F a || < °° by the uniform boundedness 
principle, and trivialiter ||F a ^|| — » for all i; e ;t . D> 

B.3 Fock space 

Lemma B.3 

Suppose w\ , W2 , ■ ■ ■ , w m £ %i,Wi^0 (all i). Then there exists au£ %_ such that (u | w,-) 7^ 
(all i). 

As a consequence: if f\ ,...,/„ are n different vectors in 3C, then there exists au € X such 
that (u\fi) (u\f n ) are n different (complex) numbers. 
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OLet Hi :={f e X ■ (/|vf,-) =0}, i = l,...,m, be the ortho-complements of C • w; in %. All of 
these are closed (linear) subspaces of X , and none of them contains an open ball (i.e. a subset 
B(e,h) := {g € X '■ \\g — h\\ < e } f° r some e > and h e X)- Therefore their union cannot be all 
of X- 

<lThis is in fact an application of Baire's theorem in finite dimensional complete metric spaces. 
Indeed, {/,■ := X \ Hi, i = 1 , . . . , m are open, dense subsets of X ■ Their intersection is open too, 
and, since only a finite number of open sets are involved, the intersection is also dense in X , 
whence non-empty. That means that the union of the Hi cannot be all of X ■ > 

Take u <E X such that u £ Hi (i = 1, . . . ,m). Then (w|w,) ^ for all ('. 

The corollary follows by taking as w, all the differences /j — /; with k^l.t> 

Lemma B.4 

Let X and H be Hilbert spaces, and suppose that for each f E X a vector e(f) € 9{ is given 
such that (e(g)\e(f)} H = exp((g|/)„). Then {e(f) : f G X.} is linearly independent. 

<Let f\,---,f n be n different vectors in X, and suppose there exists scalars d\, <x n such that 

E?=iM/i) = o. 

Choose a g € X such that (g|/i) , . . . , (g\f n ) are n different complex numbers (use lemma B.3. 
Then: 

UOieMtWi)) = (e(tg)\nOie(fi))=0 (t e R) . 
Since the system of exponential functions, {t i— > exp(fz) : z G C}, is linearly independent, we 
conclude that ai = . . . = a„ = 0. > 

Let £ (x) ■= : / € X } be the set of exponential vectors in X ■ 
Lemma B.5 

'E(X) is an independent total system off (x.)- 

<The independence follows from the lemma above, for the totalness observe that ^-£(f/) = /(g> 
...®f (elaborate). r> 

B.4 (Cyclic) (vacuum) representations of the CCR over a 
Hilbert space %. 

Lemma B.6 

All cyclic vacuum representations are unitarily equivalent 

The norm on « ) can be expressed in term of the vacuum state. 

Lemma B.7 

Expression of inner product, and hence the norm in terms of vacuum state for cyclic repre- 
sentations 

(F 2 £l\F 1 Q,)=(0{F*F 1 ) 

As a result of the above, a cyclic vacuum representation of the CCR over a Hilbert space 
X enables one to express operator norm, strong, and weak convergence in Tq(x ) in terms 
of the vacuum state. 
Lemma B.8 

Let (W,Q.) be a cyclic vacuum representation of the CCR over X- Then: 
i. Ai -> A in T (X) strongly if and only if 0)(f* (A ( - - A)* (A, - A)f) -» for all F e 
To(x). 
Lemma B.9 

Let (W,£l) be a vacuum representation of the CCR over X in u{j{) (not necessarily 
cyclic), i.e. W : X. — ► U (M) is a unitary representation of (X,+) and Q. e tf. is a unit 
vector such that (£l\W(f)Q.) =e"2ll/ll 2 for all f G X- 
Then {W(f) : f e X} is linearly independent in u{9<). 



B.5. T-functor 
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OLet e(f) := e^ f ^W(f)£l (f € X)- We will reduce the independence of {W(f) : f € x} to that 
of {e(f):f ex}. 

First, a calculation shows that the e(f), f € X, satisfy (e(g)|e(/)) = exp({g|/)) (f,g e X) and 
are therefore linearly independent. 

Now suppose that f\, ...,/„ are n different vectors in X such that £" (X,W(/ ( ) = for some scalars 
a!,...,a„. Then: 

0= (E?a i WC/i))o = E;(a i e-ilWII 2 ) • e^W{fi)Cl 

= Z\{a i e-W>f). e {f i ). 

By the independence of the e(f), (f <E X), a,e _ ill^H = 0, hence a; = for all ;'.E> 

Corollary B.10 

Let (W, Q.) be a vacuum representation of the CCR over X in 11 (tf). 

Let Tq{x) be the linear span of{W(f) : f G X.} in Then every A eT (x) has a 

unique decomposition A — a,W (/)) . 

B.5 T-functor 

Given a contraction C : 3Q — > X2 and a representation of the CCR over Xi, we can use the 
above corollary, to define a linear map T(C) : Tq{x.\) — > ro(3Q) by 

r(c) : w(f) ^ e -2(ll/ll 2 -llc/ll 2 ) W (c/) (/eaci) . 

We set out to prove: 
Theorem B.ll 

Let (Wi,£2i) a/id (1^2,^2) he cyclic vacuum representations of the CCR over Xi and Xi 

respectively . 

Then 

i. The map T(C) is a completely positive map Tq(Ki ) — > To (^2). 

ii. The map T(C) is continuous with respect to the operator norm topology, and it has 
therefore a unique extension to the C* -algebra generated byFo(Xi) and To{X2)- 

Hi. The map T(C) is continuous with respect to the strong operator topology, and it has 
therefore a unique extension to the von Neumann- algebras T(Xi ) and F(X2) generated 
byFo(Xi) andFo(X2) respectively. 

iv. By continuity, both extensions inherit the complete positivity. 

The proof of theorem B.ll rest on two pillars. 

First, in a cyclic vacuum representation one can express convergence in terms of the vac- 
uum state (see lemma B.8). Second, we can express the state of the system after the action 
of r(C) in terms of the state of the system before (see lemmas B.19). 

We now provide the details in a series of lemmata. To lighten the burden of lengthy expres- 
sions in the calculus of Weyl operators, we fix some notation, and derive some calculation 
rules which will thereafter be used without further reference: 
Lemma B.12 

Let (W, £1) be a cyclic vacuum representation of the CCR over %_ . Define 

a)(/):=ex P (-I||/|| 2 )e[0,oo) (f e x) , 
o(f\g) : = exp(-flm (f\g)) e {z e C : \z\ = 1} 

Then for f,gJ l J 2 ,fc,U E X: 

i. (Q|W(/)Q> =©(/); 

ii. W(f)W(g) = c(f\g)W(f + g); 

Hi. o(g\f) = c(-f\g) = o(/| -g) = c(f\g) = c(f\g)- 1 ; 
iv. a(/i +/2I/3 +/ 4 ) - a(/i I/3) a(/i |/ 4 ) a(/ 2 |/ 3 ) a(/ 2 |/ 4 ); 
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v. «>(/! +/ 2 )a(/i|/ 2 ) = a)(/ 1 )(0(/2)exp(- (fy\f 2 )); 

vi. co(-/) = (0(f); 

vii. W(f)* = W(— /), si/ice W is unitary. 

viii. C0(/i+/ 2 +/3+/4) 

= C0(/i)(0(/2)(0(/ 3 )(0(/4) 

xexp(-Re (/ 1 |/ 2 ))exp(-Re (/ 1 |/ 3 ))exp(-Re (/i|/ 4 )) 
xexp(-Re (/ 2 |/ 3 })exp(-Re (/ 2 |/ 4 ))exp(-Re (/ 3 |/ 4 )); 

ix. W(/i)W(/ 2 )W(/ 3 )W(/ 4 ) 

= |/ 2 )a(/i |/ 3 )o(/i |/ 4 )a(/ 2 |/ 3 )a(/ 2 |/ 4 )a(/ 3 |/ 4 )W(/i +/ 2 +/ 3 +/ 4 ); 

Conclusion: for /i,/ 2 ,/ 3 ,/ 4 e 3C we have 

(i2|w(/ 1 )w(/ 2 )w(/ 3 )w(/ 4 )|i2)-nt 1 co(/i-)n 1 < ! - < ,< 4 ex P (- (/;■!/,)) 
= co(/i)<»(/ 2 )<»(/ 3 )<»(/ 4 ) 
exp(- (/i|/ 2 )exp(- (/i|/ 3 ))exp(- (/\|/ 4 » 
exp(-(/ 2 |/ 3 ))exp(-(/ 2 |/ 4 )) 
exp(-(/ 3 |/ 4 )). 

< (viii) Observe ||£*/i|| 2 = (lf =1 /}| J^, /,) = E (J (//I/;) = E? ||/d| 2 +Il<,<;<4 [</<'l/j> + (fj\fi)] 

Lt\\fi\\+'Li<i<j<4^({fi\fj))- 

(ix) Working from the right to the left we find 

W(h)W(f 2 ) [w(f 3 )W(f 4 )] =W(f l )W(f 2 ) [o(f 3 \f 4 )W(f 3 +f 4 )] 

= o(/ 3 |/ 4 ) W{h) [o(/ 2 |/ 3 +f 4 )W(f 2 +/ 3 +/ 4 )] 

= 0(/ 3 |/4)0(/ 2 |/ 3 +/ 4 ) [o(h \f 2 +/ 3 +/ 4 )W(/j +/ 2 +/ 3 +/ 4 )] 

= o(/i|/ 2 )o(/ 1 |/ 3 )o(/ 1 |/ 4 )o(/ 2 |/ 3 )o(/ 2 |/ 4 )o(/ 3 |/ 4 )W(/i+/ 2 +/ 3 +/ 4 ). 

V „ " „ " v ' 

For the conclusion observe that o(/|g)exp(— Re (f\g)) = exp(— and then combine (viii) 

and (ix). E> 

We now return to our map T(C) : rrj(9Q ) — * Tq{%2)- 
To begin with, we calculate T(C) [Wi (f )Wi (/)] : 

Lemma B.13 

Take f J' e %±. Then: 

r ( c)[ Wl (/)w,(/)] = m ^cf^flf)^ {CfWCf) ■ ,B1) 

In particular, 

r(c) [Wi (fywx (/)] = • w 2 (c/') w 2 (c/) , 

where 

is a positive definite kernel. 

<]By the above calculation rules we have co(/i +/ 2 )o"(/i |/ 2 ) = co(/i)ff)(/2)exp(— (/] |/ 2 )), so that 

r(c) [w, <f)Wi (/)] = a(/'|/)r(c) (/'+/)] 



^TWf {mW2{C/+Cf) 



c» 2 (C(/' + /) o(C/'|C/) 2 

o) 1 (/)co(/)exp(-(/|/)) 
co 2 (C/')co 2 (C/)exp(-(C/'|C/}) 



W 2 (C/')W 2 (C/). 
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The second part of the lemma follows by observing that W(f')* = W(—f) and that both 

(f fl, ; ^ (/X/) 

(J ' J> co 2 (C/')a) 2 (C/) 

and 

(/',/) ^exp(+(/'|(7-C*C)/» 
are positive definite kernels, and so their product is. t> 



B.6 Complete positivity of r(C) 

As a result the map T(C) is completely positive: 
Lemma B.14 

LetFi,...,F JV er (3Ci), a/id^i,...,^ e # 2 . T/ien 

I^i(^|r(c)[^]^) >o, 

i.e. r(C) : ro(?Ci) — > Tq{%2) is completely positive. 
Proof 

To facilitate the reasoning, we employ the following notational trick: we write F p — Y,k=n p +i a kW(f k ) 
for suitable Q=:n\< ... < n^+i =: K, and let andr^, for k G {n p + 1, . . . ,n p+ i}, be iden- 
tical to 
Then 

I M (^|r(c)[F;F p ]^) 

= I^=i (y r(c) [( e vh +1 a/ w(/o) * (e£J +1 «)]?,) 
= E^iE^ + iE^ + iW^-ia/JO-<%|w(c/ I )*w(c/ t )Tii t ) 

= l£,=i<^oc* •!,(/},/*) • (W(C/ ; )ii/|W(C/,)iu> , 
Now, the matrixM= ((VK(C/;)r|/|VK(C/i : )ri (: ))^ ;=1 is positive definite. Letting N = N* be 
the (positive) square root of M, we can write: 

(W(Cf,)m\W(Cf k )r\ k ) = LLi(NN*) kl = i m Awv,; ; = z m W km , 

and hence: 

= Ll l=1 aiak-L(fi,f k ) ■ (Zi=xN^N km ) 



Lk,i (oc;A^ /m ) • (a k N km ) -L(f h f k ) 



>0, 



because L(-, •) is a positive definite kernel. □ 
Corollary B.15 

LetFi,...,F JV er (3Ci) andGi,...,G;vero(aC2). T/ie/i 

I^ =1 G*r(C)[F;Fj G p >0, 
symbolically: 

F*Fi...F*F N Gi 
Gi,...,Gjv-r(C) •-. • : >0 
f;Fl..f;f n G n 
<3Let £ G #2- Setting i;„ := G„£, yields 

(41 (l^=i g; r(c) [f;f p ] g p )^) = e m <^|r(c) [f;f p ]^> > o , 

by the complete positivity of T(C) . D> 

An elementary consequence of the above is: 

Lemma B.16 (Schwartz inequality for T(C) : r (Xi) — » r (?C2)) 
LetFer (^). Then [r(C)(F)] * • T(C)(F) < T(C) [F*F)] . 

<Take N = 2,Fi=I,F 2 = F, Gi = T(C){F), G 2 = -I in B.15, and do not forget that F{C) [F*] = 
[r(C)F]*.> 
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B.7 Continuity of T(C) 

Lemma B.17 

Let (Wi,Q.i) and (W 2 ,£l 2 ) be cyclic vacuum representations of the CCR over %\ and % 2 
respectively. Denote the corresponding vacuum states by C0i and 0) 2 respectively. 
Take f J' € K\, andg',g e Xi- Then: 

co 2 (w 2 (g') (r(c) [Wt (f)Wi (/)] ) w 2 (g^j 

= C o 1 (CV)co 1 (C^) eXp( - fel(/ - CC 

x 001 (Wi (CY) Wj (/') Wj {f)Wi (C*g)) . 

< Filling in the above expression and using the result of the above lemma, we find 

0)2 (w 2 ( 8 ") (r(c) [Wi {f')Wi (/)] ) w 2 (g) 

(0! (/')«! (/)exp(- (/'I/)) 



Q)2(C/')co 2 (C/)exp(-(C/'|C/)) ; 
x co 2 (g')a)2(C/')a)2(C/)co 2 (g) 

xe -{g'\Cf')) e -{g'\Cf)) e -{g'\g) e -{Cf\Cf) e -{Cf\g) e -{Cf\g) 
0>l(/')K)l(/)a)2(g')M2(g) 

xe -(/'l/> e -<s'lc/'> e -(g'|c/> e -te'|g> e -(c/'|g> e -(c/| g > 
co 2 (g')co 2 (g) exp(- (g'\g) 



coi (C*g')<»i (C*g) exp(- {C*g'\C*g) ) / 

X CO! (CV)Q)i(/')C0i (C*g) 

x e -<cVI/'> e -<c*g'|/> e -(c«g'|c*g> e -(/'|/> e -(/'|c*g> e -(/|c«g> 

CQ2(g')co 2 (g) exp(-(g'|g)) \ 
co 1 (C*g')co 1 (C*g)exp(-{Cg'|Cg)) J 

x to! (w, (cy ) iv, (/') Wi (/) m (c* g ; 

where in the one but last equality we have divided and multiplied by C0i (C*g')<Bi (C*g)e~( c f \ c f) , 
and in the last equality we have used lemma B.12 for the identification with co^ (•). D> 

Corollary B.18 

Take f J' e % x ,andg\g <= % 2 . Then: 



£2 2 



iW)(nO[W)W)])w 2 (s) 



n 2 



= K(g',g) ■ (Cli W* (C*g') W* (/') Wi [f)Wi (C*g) 

where 



K «>* = ^i^ Btp(+(g,|(/ - 0C>) ^ ) ( *' * e %) 

is a positive definite kernel. 

< Since IV* (A) = fc), we only need to replace g' and /' in lemma B.17 by — g' and — /' respec- 
tively. Since ©(—/?) = co(fc) the factors co(-) do not change, however the factor exp(— (g'| (/ — CC*)g}) 
becomes exp(+ (g'| (/ — CC*)g}). 
As the positive definiteness concerns: both 

C0 2 (g') (02(g) 



(g',g) ' 



co 1 (C*g')coi(C*g) 



B.7. Continuity of T(C) 



95 



and 

(g',g)^exp(+(g'|(/-CC*)g» 
are positive definite kernels, and so is there product. C> 

The following tool is crucial for establishing the operator norm and strong continuity of 

r(0:ro(s&)->r (sc2). 

Lemma B.19 

Let G e To(Kj) be fixed. 

Then there exist G\, . . . ,G m £ rn(3j) (only depending on G) such that 

(^(g^(t(c)[f*f])g^ = f>i(c£ [f*f] G k ) (fer (aa)) ■ (*) 

Alternatively formulated: 

Let be in the dense linear subspace Tq(%2)^2 of H2. Then there exist Jji, . . . ,\ m in the 
dense linear subspaceTo(%x)Q.\ of 9f\ (only depending on £,) such that 



k=l 



and 



(t(c)[f*f})^ 



#2 



^2 k=\ 



(Fer (3Q)) 



which implies by the Schwartz inequality that 

2 



[T(C)F]% 



^ k=\ 



(fer («a)). 



(%i) 



(%2) 



(%3) 



< Proof of (*) 

Looking at corollary B.18 while keeping in mind that 0>i,Cfl2, and T(C) are linear, we see that for 

eacfiF = if =1 a,w(/,)er (a:i) 



(02[Wi(g')(r(C)[F*F])w 2 ( g ) 



= K(g',g) • a)! (w*(C*g') F*F Wi (C*g)) . 
Consequently, if we fix G = £? =1 P 7 W 2 (g 7 ) £ r (a^) then 



(#) 



(02\G*{r(C)[F*F]jG 

= £P;%,ft) • «»l (Wi (C*g/)* [F*F]W, (C*g,-)) 
J,/ 

Since (g'.g) h-» K(g',g) is positive definite, the m x m matrix (K(gf,gj))™ j, =l is positive definite: 
let A' = N* = (N(gf ,gj))"j f-i ^ e ' ts (positive) square root. Then 

= LT= 1 N(8k,8f)N(g k ,gj). 
Plugging this in (#) yields 

co 2 (g* (r(C) [F*F] ) G) = E? =1 coi {pi [F*F] G k ) , 

where 

Gk = L7VjN(g k ,gj)-W l (C*g j ), 
only depends on G = YJJ $jW 2 (gj) (via C and K(-,-)). 
Proof of (%1- 3) 

(%2) is a mere restatement of (*), %1 follows by taking F = I (realizing that T(C)I = T(C)W(0) = 
I), and (%3) results from (%2) by using the Schwartz inequality (see B.16): 

|| [r(c)F]^|| 2 = [r(c)F]* [r(c)F]4) < ($1 (r(c)[F*F])i;) 
= ||(r(c)[F*F])^|| 2 . 
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Corollary B.20 

The map T(C) : TofXi) — ► ^(2(2) is a contraction with respect to the operator norm on 
2 (#1) and <B (#2) respectively, i.e. 

l|r(c)F|| sK2) < \\F\\ Hf(i) (Fgr («a)). 

In particular, T(C) can be uniquely extended to an operator norm continuous map between 
C* -algebras generated by To and To (2(2) respectively. The resulting extension is again 
completely positive. 

<LetF er (aa). We prove that ||[r(C)F]£|| < ||F^|| e # 2 ) in two steps. 
(Step!: ||[r(C)F]^|| < ||F|| ||£|| for % in the dense linear subspace r (9C2)n 2 of rt 2 -) 
Take i; £ Fo(3C2)f22. With i;j , . . . ,i; m as in lemma B.19 we have (by (%3) of the same lemma) 

II [r(c)F]^|| 2 < Ef=i ||^|| 2 < ll^ll 2 -ir= 1 llC'll 2 = I|f|| 2 ||y 2 . 

(Step 2: \\[T(C)F}^\\ < \\F\\ totalise X 2 ) 

Let % e H 2 and let e > 0. Take ^ € r Q (^)0. 2 with ||£- £o|| < e. Then 

1 1 [r(c)F] %\\<\ I [r(c)F] (5 - So) 1 1 + 1 1 [r(c)F] ^ 1 1 

< ||r(C)F|| WS-bW + \\F\\ \M < \\T(C)F\\E+ \\f\\ (ii^o-SII + \\U) 

<(||r(C)F|| + ||F||)e+||F||||^|| , 
and since e is arbitrarily, the proof is completed. E> 
Now for the strong continuity. 
Corollary B.21 

The map T(C) : TolXi) — > ro(5C2) can be uniquely extended to a map between the von- 
Neumann algebras r(?Ci) and r(9yj2) generated by To(?Ci) a/id To(?C2) respectively. The 
resulting extension is again strongly continuous and completely positive. 

< First recall lemma B.19: 

there exist m £ N, and %\,...,% m such that 



Sero(?(2)n 2 : 



[noF^H^jxjFiy 2 (Fer (% 



(&) 



Definition of the extension ofT(C) to r(sCi) andr(aC2) 

Let F e r(3Cl). We first define T(C)F on the dense linear subspace ro(3£2)£2 2 of 2f 2 and extend 
it to all of 9l 2 . 

Let i; 6 Lo(5C2)n 2 , and let ^ m as in (&). Next, choose a net (F«) a in T§(%x) that strongly 

converges to F. Then 

|| [r(C)(F a -F a ,)]S|| 2 < E™=i ll(F„-F a ,)4|| 2 -0 (a, a' -co) , 
so ( [r(C)F a ] i;) is a Cauchy net in 5/2, having a limit T|. 

If iG\)\ were another net in L ( )(3Ci ) converging strongly to F then again by (&), 
|| [r(C)(F a -G,)]^|| 2 < y%=\ \\F a ik-G£ k \\ 2 - (OA—) , 

since both (F a ^) a and (Gj^i;^ converge toF^ for k = 1, . . . ,m. Therefore, the net ([r^GjJi;)^ 
also converges to T|. 

Furthermore, since £f =1 ||^|| 2 = ||f;|| 2 , we see that 



lim„ [r(C)F„] 5|| 2 =lim«|| [r(C)F a ] \ \ | 2 < £JL 1 lima 1 1 F a ^ 



,2 



=Lr=i||^|| 2 <i|F|| 2 -ir=i||^|| 2 = ll F ll 2 -^n 2 ' 

Conclusion: we can define a map F(C)F : Tq(%2)^-2 ~ * ^2 by setting 
' [T{C)F\^ = Uma[r(C)Fa]% ($er (s&)). 

whenever (F a ) a is a(ny) net in ro(3Ci) that strongly converges to F , 
and F(C)F can be uniquely extended to all of tt 2 , with 

ll r (C)F|| 3(;tf2) <||F|| a(;tfi) (FeHad)), (&D 
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and 

there exist m e N, and %\ , . . . , i; m such that 



(&2) 

[r(c)F^|| 2 <vj =1 ||Fiy 2 (Fer(xi)). 



Strong continuity ofr(C) : F(Xl ) -*■ r(%) 
Suppose F a — > strongly in $ (Wi). Then 

||F„1|| 2 ^0 (5ero(9Ci)n 1 ),andsup a ||F a || sKi) <co. (&3) 

To establish that T(C)F a — > strongly in l B(H2), it suffices to check (see lemma B.2) that (i) 
sup a || [T(C)F a ] \\ m{Jb < co, and (ii) || [T(C)F] i;|| 2 - for all % e r (aC2)£> 2 . 
Now (i) follows from (&1) and (&3) above, and (ii) from (&2) and (&3). 

Complete positivity ofT(C) : F(X\ ) -» r(ffo) 
By lemma B.14, we have 

Ep, e =i(^l r (c)[^^]^> >o (Fi,...,F ll ero(«a) > §i,...,§Ye^2), 

and since T(C) (and the operator multiplication) is strongly continuous, this inequality extends to 
Fi,...,F„er(9Ci).> 



B.8 Properties of the T-functor 
B.8.1 The functorial character of T 

The map F from the category of Hilbert spaces with contractions to the category of quantum 
probability spaces (von Neumann algebras with a normal state) with quantum operations 
(completely positive, state and 1-preserving maps) is functorial as expressed by the follow- 
ing theorem. 

Theorem B.22 

LetCi : %\ — > %2 and C 2 : X2 — > be a contractions. 

ThenT{CiC 2 ) = T(Ci)T(C 2 ). 

<3 A one-line writing exercise shows that 

r(d) (V(c 2 )f) =r(CiC 2 )F 



for F = W(f) with / e Xl- By linearity of r(d)r(C 2 ) and T(CiC 2 ) this equality extends to 
F e r ( )(aCi ) and by strong continuity also to F e T(9Ci). > 



B.8.2 r(C) if C = 7 is an isometry 

Recall that a linear, adjoint-, product-, (and unity-)preserving map between Von Neumann 
algebras is called a * isomorphism. 

Theorem B.23 

Let J : %_\ — > %2 be an isometry. Then j := F(J) : r(9Q) — » T(aC2) is a/i injective *- 
homomorphism, so a quantum random variable. 

< We first prove that r(7) is an *-homomorphism. 

Let /',/ € 3d. Since 7 is an isometry, we have (Jf'\Jf) = (f'\f), and ||y/|| 2 = ||/|| 2 , so to(/) = 
co(7/) (here we use the notation co(/) = exp(— j ||/|| 2 ) of lemma B.12 again). 
Using this in formula (B.l) (lemma B.13, p. 92), we see that 

T(J)[W(f)-W(f)] =W(Jf')-W(Jf) = [T(J)W(f')] ■ [T(J)W(f)] , 

and 

r(j)[w(fr]=r(j)[w(-f)] = $zf} J w(-jf)=w(jfy = [t(j)w(/)Y , 

so that by linearity is an *-homomorphism on Tq(Xi ), and by strongly continuous extension 
also on r(3Ci). 

For proving the isometric character, we turn to corollary B.l 8 with g' = g e X2, which yields, 
using that T{J) is an *-homomorphism, 
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(w 2 ( g )a. 2 1 [r(y)m (/)]*• [r(J)Wi(/)] | w 2 (g)n 2 

x (w, | iv,* (/) • (/) | IV! {j* g )ai ) (/,/' e % ) 

Since T(J) (and the inner product) is linear, we derive from the above that 
(w 2 (g)£l 2 | [r(J)F]*-[T(J)F] | W 2 (s)£> 2 ) 
cofe) 2 

-^i>.cx V (( g \(i-jr) g )) 



o>(J*g) 2 

x I f*-f) I ) (f e r (% 



and by continuity this extends equality extends to F £ r(3Q). In other words, 

||[r(7)F] (V 2 (g)n 2 )| 2 = 



(fl(7*g) ; 



■exp((g|(/-77*)g» 1^(^(7*^!) || 2 (Fer(3ti),«6 3!2). 



>0 

Therefore, 



[r(7)F] (w 2 ( g )n 2 ) =o f(wi(7^)£2i) =o (Fer(sQ),*es6). 

Since 7* (sq) = 3Q (7*7 = 1^, ) we see that if r(7)F = on 9i 2 (so on the dense linear subspace 
ro(^2)£22 of H 2 \ then F = on the dense linear subspace Tq{%_\)0.\ of so F = (by con- 
tinuity). In other words, the linear map T(7) : r(aj) — > T(ac 2 ) has kernel {0} i.e. is injective. 

> 



B.8.3 r(C) if C = J* is a projection 

Let j be a injective *-homomorphism from (-#,<])) into (-#,<|>). Recall that a map P : 
(j? , <|)) — ► (A , (|)) is called a conditional expectation with respect to /' if it is the (necessarily) 
unique quantum operation satisfying Po j = I A . 

On the level of Hilbert spaces, we have that if 7 is an isometry from H\ to J(^i) C #2, 
then it adjoint 7* = P is a projection from #2 to H\ if we identify #1 ~ J(tii ). 

Theorem B.24 

7f J* is a projection (J an isometry) % to 3C then P := T(J*) is a conditional expectation 
(by definition ). 
In particular: 

^ Bl ,B 2 er(^Aer(K) B i p (A)B2=P(j(Bi)Aj(B 2 )) 

< We check the claim for B, = W(bt) € r(%) and A = e r(ac). By linearity and continuity 
the claim then extends to arbitrary fl; e T(aCo) and A e r(ac). 
Let ii , b 2 6 3Q) and a e 3C . Since for f,g,hm (Ki respectively): 

W(/)W(g)W(A) = e -''Im(/l8> e -'1m(/l''> e -''I™<«|A>lV(/ + g + ft) . 

we have 

W{b\)-T{J*)W(a)-W{b 2 ) 

= e -i(ll«ll 2 -ll^ll 2 ) e -^il^) e -^|62) e -^IM W(fol +fl + fc2 ) . 

On the other hand, T(J)W{bi) = W(Jb,) since 7 is an isometry, and it follows that 
T(J)W(biyw{a)-T{J)W{b 2 ) 

= W(Jb l )W(a)W{Jb 2 ) 

e -ilm(Jb l \a) e -ilm(jb 1 \jb 2 ) e -am(a\Jb 2 ) w ( Jhl +a + Jb 2 ) , 
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and therefore, 

r(j*)[r(j)w(bi)-w(a)-r(j)w(b 2 )\ 

= e -\(\\Jb l +a+Jb 2 1| 2 - \\r (Jbt +a+Jb 2 ) || 2 ) 

(* ) 

x g-ilm {bi \J'a) g-ilm {hi \b 2 ) g-ilm (-/*a|*2> 

Using that 7 isisometry and that (Jf\g) = (f\J*g), the term \\Jb\ +a + Jb 2 \\ 2 - \\J*(Jb\ +a + Jb 2 )\\ 
is seen to be ||a|| 2 — ||/*a|| 2 . Further J* (Jb\ + a + Jb 2 ) = b\ +Ja + b 2 , so that (*) equals (*')•!> 

B.8.4 r(C) if C = U is a unitary map 

Since f/ : 9Q — > as well as U^ 1 = U T : — > ^Ci is an isometry, we have that r(f/) is a 
♦-isomorphism with inverse r(I/) _1 = r(f/^). 

In the particular case that sq = %2 =: r(J/) is called an *-automorphism of 
Example B.25 

The ^-automorphism of £ (R") are of the form B -> f/Bf/ f where t/ : M" — > M" is itself a 
unitary map. 

The above example illustrates that the time evolution of a quantum system in $ (W) im- 
plemented by a unitary map, is in fact a ^-automorphism. Via this identification, time 
evolutions on more general Von Neumann algebras can be defined. 

B.9 The usual approach to second quantization for the 
free (quadratic) scalar field 

The approach taken here largely follows [9, p. 107-1 14]. 
Suppose the Lagrangian density is given by 

with V a differentiable real function. The corresponding classical field equations (of mo- 
tion) read 

For simplicity we now restrict to the case of a quadratic potential (i.e. a harmonic oscillator 
potential), V ((])) = jm 2 (j) 2 , which has the Klein-Gordon equation 

□ (j) + m 2 (|) = , 
as classical field equation, and 

as conjugate field (i.e. field impulse). 
From the above, the Hamiltonian is given by 

H = jd\{\ [n(xf + (V^)(x)) 2 ] + m>(x) 2 } , 

which can be recognized as a continuum of coupled oscillators. To elucidate this, assume 
that space has one dimension (coordinate x) and that the coordinate x can only take discrete 
values which are integral multiples of a an elementary length taken as unity. In that case 
the Hamiltonian 



2 
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H = 



(B.2) 



is replaced by 




2 L {^n + ^n+l-^nf+m 2 ^ 



(B.3) 



with commutators 

[7C„,7t„/] = [<t>n,<M =°: and M> B ,JC„/] = • 

A physical model having this Hamiltonian is a chain of mass points with the following 
characteristics: 

1 . the displacement (along the x-axis) of the n th mass point from its equilibrium is given 
by (j)„, its impulse by jc„; 

2. mass point n is attached with a spring to its equilibrium position, the restoring force 
being provided by the m 2 ^ 2 term in the Hamiltonian; 

3. mass point n is coupled to its two neighbors n — 1 and n+ 1 which induces a additional 
force incorporated by the ((j)„ — (j)„±i) 2 contributions to the Hamiltonian. 

Such a chain could model the vibrations of a one-dimensional crystal where (j)„ is the dis- 
placement of the n th atom from its equilibrium, and each atom is coupled to his equilibrium 
position and to its neighbors. Like in the case of a string (which can also seen as a con- 
tinuum limit of a chain of oscillators), it is fruitful to look for normal eigenmodes for this 
model. 

Instead of taking a Fourier series expansion of the discrete version of the system (B.3) 
and then taking the continuum limit, we shall directly apply a Fourier transform to the 
continuous system (B.2). 



Using that 




and that formally 




— — e' 



.ik+k'x 




we see that 




where the last equality uses that %{x) is real. 
Likewise fdx§(x) 2 = /|| tf(k)§{k), whereas using 




we see that 




Putting pieces together a Fourier transformed Hamiltonian obtains: 
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H=^jdki^\k)%(k) + [m 2 + k 2 ]^\k)^{k)\ . 



(B.4) 

o»(*) 

The resemblance with a continuum of harmonic oscillators 

H(k) = ^?^+ma(k) 2 q(k) 2 ^ withm=l 

begins to arise, except for the fact that §(k) and %(k) are no real quantities, hence cannot be 
quantized to self-adjoint operators. However, %^(k)%(k) and § T (k)§(k) are, and one could 
quantize the (positive) square root of those. 

Here we will adopt another approach by quantization of the corresponding creation and 

annihilation operators (which are not required to be self-adjoint). 

Let 



a(k) :-- 
a\k) = 



1 



y/4KK)(k) 
1 



(o(k)^(k)-in\k) 



Then 



Jdko)(k) a t (k)a(k)+a(k)a f (k) 

= -L Jdk cd(*) 2 ($*(*)?(*) +?(*)?*(*)) 

+ (n^k)n(k)+n(k)n^(k^ 
= T- f dk 



(B.5) 
(B.6) 

(B.7) 
(B.8) 



C0(fe) 2 $ t (fe)$(fc)+7t t (fe)7t(fe) 



where we have used in the first equality that the "mixed" terms compensate each other 

srJHk)K(k)dk = s:j(-k)n(k)dk = srj(k)n(-k)dk 

=s:J(k)%Hk)dk, 
r^(k)mdk =...= s?ji(k)$(k)dk, 

and in the second equality that the quadratic terms in and % respectively double 

srj(k)^(k)dk = s:Jm(-k)dk = f?j{-k$(k)dk 
=rJ*m(k)dk, 

j:„n(k)nHk)dk = j:«,n(k)tn(k)dk. 
Summarizing 



1 f°° r 
H=- dk co(fe) a\k)a(k)+a(k)a\k) 



Quantization now proceeds by postulating 

[a{k)^{k')] =8(k-k'), [a{k),a{k!)] = [a\k),a\k')]=Q. 

Using the commutator product rule [01,0203] = [01,02)03 + 02(01,03], 

[H,a(k )] =-(o(k )a(k ), 

[H^iko)] = +(S)(k )a(k ), 
so that, if \E) is an eigenstate of H with energy E, then 



(B.9) 



(B.10) 
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H 



(a(ko)\E)) 



a{k )H \E) + [H,a(ko)} \E) = Ea(k ) \E) + ~(o(k )a(k ) \E) 
[E -©(*&)] (a(*o)|E>) 

a\ko)H\E) + [H,a\ko)] \E)=Ea\k ) \E) + 0){ko)a(koY \E) 
[£ + co(fc )] («%)!£)) 



H 



(a\ko)\E)) 



i.e. a(k) and cfi(k) lower respectively raise the energy of the system by an amount of (0(k), 
which means that with respect to energy quanta, they play the role of annihilation and 
creation operators respectively. 

Here we see the relation between particles and fields: The field (B.2) consists of a contin- 
uum of (uncoupled) harmonic oscillators, each with its own Hamiltonian 



and corresponding creation (annihilation) operator a'(k) (a(k)); each oscillator contributes 
an integral multiple of its eigen-energy (o(k) to the total energy of the field. Therefore the 
energy of the system can be considered as made up of energy quanta of different kinds 
(labeled by (i)(k)), each with its own occupation number n(k). 

Eigenstates of the field are normally not taken against the field operator (j) (position), or for 
that matter, the field impulse operator Jt, but against the energy operator H. The ground 
state |0) of H is the state with all oscillators H(k) = 0)(k){a^ (k)a(k) + a(k)cft(k)] in their 
ground state i.e. with energy \(S)(k) and occupation number n(k) = 0. Working upwards, 
all other eigenstates (excitations) are obtained by creating energy quanta of different (fre- 
quency) modes. 

However, since we are dealing with a continuum of oscillators (each with non-zero ground- 
level energy 0)(k)/2), the ground-state of H has an infinite energy. Indeed, since 



(note the difference between 0, no state i.e. there is no field nor oscillators, and |0), the 
field exists and is in its ground state with all oscillators in their ground-level energy) we 
can evaluate 



although it would be more appropriate to say that the above is undefined by the appearance 
of 8(0). Since the ground-level energy of the field is not observable (unless the field is 
destroyed), and only differences of energy levels are physically observable, we can formally 
declare the ground-level energy of the field to be 0, and furthermore, take the ground-state 
to be normalized, (0|0) = 1. 
A manner to accomplish this is to split up H as 



H(k) = ico(fe) a*(k)a(k)+a(k)aT(k) 



a 



l(k)a(k)\0) =0, and a{k)a\k)\0) = |0) 
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where we have used in the third equality that [a(k),cft(k)] is a scalar ("c-number" in physi- 
cal terminology). 

The above way of handling infinite expected values in the ground-state, i.e. rewriting occur- 
rences of a(k)a f (k) to a* (k)a(k) + [a(k),cft(k)] (the latter term of which is formally 8(0)) 
and then declaring the [a(k),a^{k)\ to vanish in the ground-state expectation, is referred to 
as putting the operator in "normal ordering" (also called Wick's ordering). 
Having declared away the infinite ground state energy, problems arise again when forming 
excited states such as a r (k) |0) (one energy quantum in mode 0)(k)), since these states are 
not normalizable: 

(0|a(Jfc)flt(Jfc)|0) = (0|[fl(/t),a + (/t)]|0) = 8(0) (0|0> 

where we have again used that a(k) |0) =0. 

Wave packets are therefore build using linear superposition: a state of the form 



Zoo 
dkf{k)a\k)\0) 



a f 

is normalizable provided that 

oo > ^0|fl/a}|0^ = jdk jdk' JW)f(k) (0K£')a + (kM 

at(jfc)a(f)+8(t-f) 

= fdk\f(k)\ 2 , 

J — oo 

i.e. the "density kernel" / has to be square integrable. Assuming |0) to be normalized 
(i.e. Jdk \ f(k)\ 2 = 1), we see that if f(k) has a peak around ko with small spread, i.e. 

<0|jfca/a^|0) = Jdkk\f{k)\ 2 = k , 

and 

(0| (k - h) 2 ■ a f a}\0) = jdk (k - k Q ) 2 \f(k)\ 2 

is small, then a|(fc) |0) can be interpreted as a wave packet around wave number k = ko. 
Using that 

[a{k),a\] - jdk' f(k') [a(k),a\k')} = f(k) (B.ll) 

8(*-f) 

we see that 

fl(Jfc)flJ|0)=fl}o(^|0) + [fl(Jfc),flJ]|0) =/(Jfc)|0) . 

~wT 

If we use the latter in the first and third equality of 



(B.12) 



0|fl / fl t (/t)fl(/t)a}|0^ = f{k)(Q\a f a\k)\<yj = f{k)(Q\a(k)a\\0 
= f(k)f(k)(0\0) = \f(k)\ 2 , 
we derive that the energy in state a J- 10) is 

/•oo 
dk (i>{k)(0\a f a\k)a{k)a\\0) 

/•oo 
dk o)(k)\f(k)\ 2 k 0)(k ) 
-oo 



(B.13) 



as one would expect from a state that is concentrated around k = ko- 
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In conclusion, to obtain meaningful expressions in quantum field theory the fundamental 
creation and annihilation operators a\k) and a(k) have to be smeared (smoothed) with 
square integrable functions in fe-space. The resulting operators a J- and ay are so-called 
operator-valued distributions. 

Being a linear superposition of one-particle states, a/|0) is to be interpreted as a one- 
particle state. Multiple particle states are created by successively applying creation operator- 
valued distributions a J. to the ground state |0). For example a (raw) r-particle state is given 
by 

\r)=a}a} 2 ...a} r \0). 

n-particle states are orthogonal to m-particle states whenever n^m. Indeed, from (B.ll) 
we have 

/■co />co _ _ 

dk g(k) [a{k),a\\ = Jdk g(k)f{k) = (g\f) L 2 , (B.14) 

/(*) 

and from there we derive successively that the vacuum is orthogonal to n > 1 -particle states: 
(0|4 i 4...4J0) = (0|44 2 ...4J0)=0, 



that 1 -particle states are orthogonal to n > 2-particle states 
(0\a gl a\a\ 2 . . .aJjO) = (0| a gl a T fi a\ . ..a\ |0) 

= <0|4 a gi a\ . . . 4 10) + (gi|/i ) (0|4 2 . . . fl }jo) = 
o 

(the last term vanishes because the vacuum is orthogonal to a (n — l)-particle state), and so 
on using that 

= [a gl , aji ] fl «2 • • • a sm + a gi [a g2 , a\]a gi . . . a gm + ... + a gl a g2 ... a gn _ , [a gn , a\ ] 
= {g\\f\)a g2 ..- a gm + a gl (g 2 |/i ) a g3 . . . a gm + a gi a g2 (g 3 |/i ) a g4 . . . a gm + ... 
+ a gl a g2 ...a gm _ l (g m \fi) . 

The number of quanta in a multiple particle state is counted by the number operator 

/■co 
dk a\k)a{k), (B.15) 
-co 

since 

[a(k),a J( fi a f2 ...a fn ] 

= [a{k),a\y h ...a fn +a} ] [a(fe),a} 2 ]a} 3 . . .a f + 

...+a\...a\ n Ja{k),a\ n ] 
= Mk)a} 2 ...al+a} i Mk)a} 3 ...a} + 

...+a\..J fn J n (k), 
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so that 

Na\a\ 2 . . . fl J |0> = fdk a\k) a(k) a\a\ . ..a\ |0) 

= fdk a 1 " (k) a\ a\ 2 . . . a) n a (k) |0) + f dk a f (k) [a (k) , a\ a\ . . . ajj |0) 



J dk h M« f (k)a} 2 •••«/„+/ dk a\ f 2 (*) fl t {k)a\ . . . a\ + 
... + fdk a\...a) n J n {ky{k)^j |0) , = n-a\...a) n _f n \0) 

Using the a\ coherent states of the field can be introduced. 

The field Hamiltonian is an (infinite) sum of independent (uncoupled) harmonic oscillators. 
Each harmonic oscillator has its own coherent states Coh(A-) := exp(A,a T (&)) |0), X £ C, 
(where k is the wave number identifying the oscillator). Since the different modes (oscil- 
lators) are independent, simultaneous coherent states of the system (coherent states of the 
field) can be formed by taking infinite products of individual harmonic oscillator coherent 
states: their form is 

Coh(/) :=K(f)exp(fdkf(k)a^k))\0) , 
where K(f) is a normalization factor and f(k) is the weight. 
We now calculate (Coh(g)|Coh(/)). 

SetA:= [fdk g(k)a f (k)] f = fdk g(k)a(k) and B := fdk' f(k')a f {k'). Then 
[A,B] = fdk fdk' gjk)f(k') [a{k),a\k')] = fdkgjk)f{k) = (g\f) L 2 , 

S v ' 

S(k-k') 

whence A and B commute with [A,B], so that by the Baker-Hausdorff-Campbell formula 

e A £ B = e A+B e ^[A,B]^ L i k ewise, e B e A = e B+A e^ B ^ = eVe"^, so that e A e B = 
As a result, 

[Z(i)^(/)]- 1 (Coh(g)|Coh(/)) 
= (0| (exp( fdkg{k)a\k))y (exp(jdk' f(k')a\k'))) |0) 

= (o| e A e B |o)- e ^(o| e V|o) = e ^](o|o) <B '' 6) 



= exp(Jdkg(k)f(k)) , 



where we have used in the before last equality that A = fdk g(k)a(k) annihilates |0) so that 
<^|0) = (I+A+A 2 /2\+A 3 /3l + ...)\0) = |0), and by the same token (0|e B = (0|. 

In particular, K(f) = exp(-± /dA; /(£)/(£)) = e ~^ m l 2 , and therefore 

(Coh( g )|coh(/)) =e -lii/n 2 -iii*ii 2 +a'i/> =e -.im( g \f)Au-f\\ 2 . 

From (B.16), we see that the coherent vectors are not orthogonal. Nevertheless, they form 
a complete set, since 

|Coh(;/) |,=o=4l°>. 
and more generally, 

&Coh(f/) | f =o=(4)"|0) , 
from which via a polarization formula, 

aia 2 ---an^ (2 n n!) _1 £ ei) ... iEj!=±1 (eiai + . . .+£„«„)" , 
the state 



4i4 

can be obtained. 
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